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ABSTRACT 


This  annual  report  presents  a  summary  of  the  scientific 
progress  and  accomplishments  on  research  projects  funded  by 
the  Joint  Services  Electronics  Program  (JSEP)  for  the  contract 
period  from  1  April  1980  through  31  March  1981.  It  does  not 
contain  information  regarding  accomplishments  on  research 
projects  funded  in  other  ways. 

The  Joint  Services  Electronics  Program  at  the  Polytechnic 
is  the  core  of  interdisciplinary  research  in  electronics  encom¬ 
passing  programs  in  the  Departments  of  Electrical  Engineering, 
Physics  and  Chemistry  under  the  aegis  of  the  Microwave 
Research  Institute.  The  research  encompassed  by  this  pro¬ 
gram  is  grouped  under  three  broad  categories :  Electro¬ 
magnetics,  Solid  State  and  Information  Electronics.  The 
detailed  projects  (work  units)  comprising  the  complete  pro¬ 
gram  are  listed  in  the  Table  of  Contents. 

In  addition  to  the  progress,  we  include  for  each  of  the 
work  units  the  publications  during  the  past  year  resulting 
from  the  research,  and  some  of  the  recent  interactions  with 
the  DoD  or  industry  or  other  academic  institutions.  At  the 
front  of  this  Annual  Report  we  present  our  Report  on  Signifi¬ 
cant  Scientific  Accomplishments,  and  at  the  end  we  append  a 
comprehensive  list  of  all  the  journal  publications  and  confer¬ 
ence  proceedings  during  the  period  from  August  1,  1980  to 
the  present,  arranged  alphabetically  according  to  the  name  of 
the  first  author. 
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1.  Rigorous  Solution  for  the  Oblique  Guidance  of  Electromagnetic  Sur¬ 
face  Waves  by  Periodically -Grooved  Dielectric  Waveguides,  and  Its 

Many  Implications 

Professors  S.T.  Peng  and  A. A.  Oliner 

When  electromagnetic  surface  waves  are  guided  by  a  periodically- 
grooved  dielectric  layer  waveguide,  and  the  waves  propagate  in  a 
direction  perpendicular  to  the  periodic  grooves,  the  TE  and  TM  surface 
waves  are  not  coupled  and  the  boundary -value  problem  is  a  scalar 
two-dimensional  one.  When  the  surface  waves  are  guided  at  an  oblique 
angle  to  the  grooves,  however,  TE-TM  mode  coupling  occurs,  and  the 
problem  becomes  a  formidable  vector  three-dimensional  one.  We  have 
recently  obtained  a  rigorous  solution  to  the  oblique  guidance  problem 
for  a  dielectric  periodic  waveguide  with  rectangular  grooves.  To  our 
knowledge,  this  is  first  rigorous  solution  to  a  vector  problem  involving 
any  kind  of  periodic  structure. 

The  problem  which  was  solved  has  great  practical  value  and  poten¬ 
tial  importance.  First  of  all,  it  represents  a  structure  which  is  cur¬ 
rently  being  explored  in  integrated  optics  for  applications  such  as 
multiplexing.  Because  of  its  potential  for  many  other  applications, 
careful  measurements  were  made  recently  in  Germany  but  there  was  no 
theory  against  which  these  measurements  could  be  compared.  Our 
theory  provides  excellent  agreement  with  those  measurements ,  which 
correspond  to  the  Bragg  region  of  the  wavenumber  plane,  and  which 
exhibited  four  stop  bands  instead  of  the  ordinary  two  because  of  the 
added  TE-TM  mode  coupling. 

We  have  recognized  that  the  solution  we  derived  is  also  applicable 
to  a  variety  of  problems  in  the  millimeter-wave  integrated  circuit  field, 
in  particular  to  periodic  gratings  on  dielectric  waveguides  of  finite 
width,  for  which  no  solutions  are  yet  available.  We  are  currently 
utilizing  this  solution  in  the  analysis  of  two  basic  periodic  grating 
structures  of  this  type:  (1)  a  grating  whose  dimensions  place  it  in  the 
Bragg  reflection  region  so  that  it  finds  application  as  a  resonator  or 
filter;  there,  we  find  that  an  extra  (unwanted)  reflection  peak  or 
transmission  dip  is  obtained,  and  (2)  a  grating  whose  dimensions  permit 
leakage  of  radiation  so  that  the  structure  becomes  a  leaky  wave  an¬ 
tenna.  In  the  first  example,  we  find  a  new  effect  which  can  seriously 
impair  device  performance ;  in  the  second  example ,  our  new  solution 
permits  the  development  of  a  theory  for  this  class  of  antennas. 

Our  significant  contribution  here  is  that  we  have  derived  a  rigor¬ 
ous  solution  which  not  only  is  the  first  of  its  type  but  also  is  a  key  to 
the  opening  of  many  new  doors. 

2.  Extension  of  Dynamical  Theory  to  X-Ray  Phonon  Interaction 

Professor  H.J.  Juretschke  and  Mr.  F.  Robbins 

The  interaction  of  x-rays  with  many  phonons  has  been  formulated 
under  the  conditions  for  which  both  the  primary  and  the  scattered 
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x-ray  modes  consists  of  two  or  more  coherently  coupled  plane  waves. 
In  addition,  the  scattered  modes  are  treated  as  subject  to  boundary 
conditions  appropriate  to  entirely  internally  generated  fields.  This 
dynamically  self-consistent  theory  predicts  new  features  of  x-ray 
phonon  interaction  under  diffraction  conditions,  such  as  satellites  in  the 
main  reflection,  line  narrowing  and  peak  reduction  beyond  that  of  the 
Debye-Waller  factor,  and  collective  phonon-induced  changes  in  the  index 
of  refraction  of  all  the  x-ray  modes.  Some  aspects  of  these  features 
have  been  seen  in  already  published  experiments.  Others  are  below  the 
resolution  available  under  most  presently  employed  experimental  condi¬ 
tions,  but  they  will  become  accessible  with  the  appropriate  use  of 
synchrotron  sources.  The  theory  also  contains  the  full  description  of 
the  interaction  of  x-rays  with  photons  in  those  solids  where  the  photon 
propagates  within  the  crystal  in  strong  coupling  to  optical  phonons  - 
the  so-called  polariton.  This  particular  type  of  nonlinear  interaction 
between  x-rays  and  lasers  in  the  optical  range  is  the  strongest  and 
therefore  the  most  likely  to  be  experimentally  observable.  The  obser¬ 
vations  predicted  by  our  calculations  will  constitute  a  test  of  the  extent 
to  which  current  dynamical  theory,  which  is  self-consistent  but  is 
essentially  a  linear  theory,  can  account  for  such  interactions  and  also 
the  extent  to  which  a  more  microscopic  theory  of  x-ray  propagation  in 
crystals  has  to  be  sought,  especially  as  the  intensities  of  such  x-rays 
continue  to  increase. 

3.  A  Robustizing  Approach  to  System  Identification 
Professor  L.  Kurz 

The  advent  of  modern  high-speed  sampling  techniques  has  resulted 
in  additional  stress  on  the  problem  of  signal  processing  in  a  dependent 
corrupting  noise.  Since  in  many  systems  of  practical  interest  the  knowl¬ 
edge  of  the  noise  statistics  is  either  inexact  or  unspecified,  some  adap¬ 
tive  techniques  which  are  frequently  ineffective  or  too  complicated  to 
implement  were  suggested  in  the  past.  If  the  methodology  of  system 
identification  can  be  adjusted  appropriately,  a  process  equivalent  to  the 
prewhitening  filter  may  be  evolved.  To  that  end,  a  robust  identifica¬ 
tion  of  the  autoregressive  (AR)  model  was  introduced. 

The  robustizing  process  is  in  the  form  of  a  robustized  Robbins- 
Monro  stochastic  approximation  (RMSA)  algorithm  and  is  based  on  the 
m-interval  polynomial  approximation  (MIPA)  method  proposed  recently  by 
the  author.  The  resulting  algorithm  represents  a  recursive  version  of 
the  well-known  maximum  entropy  method  (MEM)  for  spectral  estimation 
introduced  by  Burg  or  of  the  popular  Widrow  least-mean-square  (LMS) 
adaptive  filter  adopted  in  many  engineering  disciplines.  Furthermore, 
the  robustized  algorithm  leads  naturally  to  a  robustized  version  of 
Akaike's  information  criterion  (AIC)  to  determine  the  order  of  the 
autoregression.  The  simplicity  of  implementation  and  flexibility  make 
the  application  of  identification  algorithms  based  on  the  MIPA  concept 
particularly  attractive  in  practice.  The  flexibility  and  robustness  of 
the  procedures  were  confirmed  by  Monte  Carlo  simulations. 
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4.  Collective  Approach  to  High  Frequency  Transmission  of  Electromag¬ 
netic  Waves  Through  Dielectric  Layers 

Professor  L.B.  Felsen  and  Mr.  P.D.  Ein2iger 

Transmission  through  layered  media  is  of  interest  for  various 
applications  in  electromagnetics  and  optics.  Examples  are  provided  by 
curved  dielectric  radome  covers  for  microwave  antennas,  and  by  layered 
optical  filters.  A  common  approach  to  this  problem,  when  the  excitation 
comes  from  a  realistic  localized  source,  is  to  trace  rays  successively 
through  the  layered  structure.  Due  to  internal  reflections  at  the  layer 
interfaces,  this  procedure  gives  rise  to  a  multiplicity  of  ray  fields, 
which  renders  computation  inconvenient  and  physical  interpretation 
awkward . 

By  a  novel  approach,  based  on  selective  Poisson  summation  of  the 
above-noted  ray  fields,  we  have  managed  to  treat  higher-order  multiple 
internal  reflection  events  in  a  collective  manner.  One  thereby  obtains  a 
representation  for  the  transmitted  field  comprised  of  a  certain  number  N 
of  ordinary  ray  fields  with  internal  reflections  0  <  n  <  N,  plus  a  "collec¬ 
tive"  ray  field  with  weighted  amplitude  that  accounts  for  all  orders  of 
reflection  n  >  N.  The  collective  weighting  incorporates  composite  plane 
wave  transmission  coefficients  for  the  equivalent  local  layer  structure 
(with  previously  ignored  curvature  corrections  for  nonplanar  layers) 
rather  than  the  transmission  and  reflection  coefficients  for  a  single 
interface,  which  appear  in  the  ordinary  ray  fields.  The  number  N  of 
ordinary  rays  must  be  large  enough  so  that  all  rays  with  n  >  N  exhibit 
the  phase  coherence  descriptive  of  local  plane  wave  transmission.  For 
relatively  thin  layers,  the  collective  ray  concept  is  found  to  be  valid 
even  for  the  direct  ray  path  without  internal  reflections  so  that  the 
entire  source-excited  transmitted  field  is  represented  accurately  by  a 
single  species  of  collective  rays. 

"Collective  ray  optics’’  may  be  regarded  as  a  generalization  of 
conventional  ray  optics  for  layered  media.  The  validity  of  the  field 
formulation  in  terms  of  ordinary  and  collective  rays  has  been  verified 
by  comparison  with  independently  derived  results  for  the  rigorously 
solvable  prototype  of  a  circular  cylindrical  layer  of  constant  thickness. 
The  investigations  so  far  have  been  restricted  to  line-source-excited 
two-dimensional  configurations.  Extensions  to  other  source  types  (dis¬ 
crete  and  continuous  aperture  distributions)  and  to  three-dimensional 
structures,  with  emphasis  also  on  polarization  effects,  are  planned  for 
the  future. 

5.  Optimal  Multivariable  Controllers  Employing  a  Two-Degree-of-Free- 
dom  Design  Methodology 

Professors  D.C.  Youla  and  J.J.  Bongiorno,  Jr. 

By  and  large,  the  design  of  optimal  Wiener- Hopf  controllers  for 
linear  large-scale  time -invariant  plants  has  been  restricted  to  one- 
degree-of-freedom  control  laws.  In  turn,  this  has  practically  precluded 
the  use  of  reduced-order  modeling  and  has  made  it  impossible  to  ade- 
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quately  incorporate  other  "hard"  design  specifications  imposed  by  sen¬ 
sitivity  constraints  and  phase-gain  stability  margin  requirements  into  a 
bonafide  optimality  theory.  The  conventional  single-loop  configuration 
satisfies  the  required  closed-loop  stability,  but  forces  the  rejection  of 
other  solutions  which,  although  they  could  provide  performance  better 
than  an  order  of  magnitude,  are  inconsistent  with  single-loop  stability. 

Guided  and  motivated  by  the  desire  to  control  large-scale  systems 
simply  and  effectively,  we  have  reexamined  these  previously-rejected 
solutions  and  have  discovered  that  the  potential  of  some  of  those  solu¬ 
tions  can  be  realized,  consistent  with  system  stability,  if  the  topology 
of  the  control  structure  is  appropriately  altered  to  include  suitable 
minor-loop  feedback.  Using  this  approach,  we  have  developed  a  two- 
degree-of- freedom  design  methodology.  By  exploiting  the  additional 
degree  of  freedom,  it  has  been  possible  to  combine  transient  and  steady - 
state  performance  requirements  with  "hard"  constraints  imposed  by 
sensitivity  and  stability  margins.  In  some  cases  we  have  obtained  an 
improvement  of  two  orders  of  magnitude  over  the  more  conventional 
one-degree-of-freedom  design  techniques!  The  possibilities  appear 
exciting  but  further  research  is  necessary  before  a  complete  theory 
becomes  available. 

6.  New  Interpretation  of  Si2e  Effects  in  the  Conduction  Properties  of 

Thin  Metallic  Films" 

Professor  H.J.  Juretschke  and  Mr.  R.  Pimpinella 

Size  effects  describe  the  influence  of  surface  scattering  of  con¬ 
duction  electrons  on  the  conductance  of  thin  metallic  films.  If  properly 
understood ,  they  contain  information  about  the  environment  in  the 
surface  region  sensed  by  these  electrons.  We  have  proposed  a  new 
interpretation  of  size  effects  that  is  based  on  the  proposition  that  the 
so-called  surface  scattering  actually  takes  place  throughout  a  small  but 
finite  layer  of  the  solid  just  beneath  the  geometrical  surface,  rather 
than  being  due  entirely  to  the  geometrical  features  of  the  surface 
boundary  itself,  as  has  been  commonly  assumed.  The  theory  has  been 
formulated  in  terms  of  two  principal  parameters  that  define  the  prop¬ 
erties  of  the  scattering  layer:  its  thickness  and  the  surface  region 
mean- free  path.  It  incorporates  contributions  to  the  surface  mean-free 
path  arising  both  from  a  modified  surface  phonon  spectrum  and  from 
static  structural  disorder.  Most  of  the  results  of  current  models  for 
size  effects  are  either  contained  in  our  formulation,  or  are  mimicked  by 
it,  though  for  different  underlying  reasons.  The  theory  has  so  far 
been  very  successful  in  three  different  applications:  (a)  the  surface 
scattering  by  clean  and  smooth  metal  surfaces  over  a  large  range  of 
temperatures,  (b)  the  influence  of  surface  charge  on  this  surface  scat¬ 
tering,  and  (c)  the  scattering  by  controllably  produced  very  thin  rough 
overlays  on  the  same  originally  very  smooth  surface.  The  results  of 
this  interpretation  imply  that  well-annealed  metal  surfaces  have  a  char¬ 
acteristic  domain  structure  of  atomically  smooth  regions  separated  by 
extended  faults,  rather  t*  i  disorder  on  the  scale  of  individual  atoms. 
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A.  MILLIMETER  WAVE  WAVEGUIDES  AND  ANTENNA  ARRAYS 


Professors  S.T.  Peng,  A.  Hessel,  and  A. A.  Oliner 

Unit  EMI-1 

1.  OBJECTIVE(S) 

To  systematically  examine  the  properties  of  a  class  of  periodically- 
modulated  dielectric  waveguides  to  assess  their  suitability  as  scanning 
antennas  for  the  millimeter  wave  range.  The  study  would  involve  first 
the  analysis  and  measurement  of  several  types  of  dielectric  waveguides 
suitable  for  mm  wave  applications.  The  results  would  then  be  employed 
in  an  investigation  of  periodically-modulated  leaky  wave  structures  for 
use  as  linear  or  two-dimensional  antenna  arrays.  These  antennas  can 
have  a  low  profile  or  be  flush-mounted,  are  inherently  rugged,  are 
simple  to  fabricate,  and  can  be  made  conformal  to  the  surface  of  a 
moving  vehicle. 

This  comprehensive  investigation  involves  the  solution  of  a  number 
of  constituent  unsolved  theoretical  problems  which  are  basic  in  nature 
and  of  substantial  interest  in  their  own  right. 

2.  APPROACH 

The  antennas  are  viewed  as  leaky  wave  structures  and  are  ana¬ 
lyzed  by  considering  them  as  waveguides  with  a  complex  propagation 
constant;  the  propagation  behaviors  are  determined  by  microwave  net¬ 
work  methods.  For  a  systematic  and  thorough  understanding  of  the 
wave  processes  involved,  the  complex  waveguiding  problem  is  broken 
into  simpler  constituent  ones  which  are  of  great  theoretical  importance 
in  themselves  and  can  be  handled  independently.  The  overall  wave- 
guiding  characteristics  are  then  analyzed  in  terms  of  the  wave  pro¬ 
cesses  associated  with  each  constituent  problem.  In  fact,  the  computer 
program  for  the  grating  antenna  problem  will  consist  of  subprograms 
separately  developed  for  each  constituent  problem. 

3.  PROGRESS 

This  research  program  is  divided  into  two  basic  parts:  the  prop¬ 
erties  of  open  dielectric  waveguides  for  millimeter  waves,  and  new  types 
of  scanning  antennas  based  on  these  waveguides .  Since  a  proper 
understanding  of  the  waveguides  must  underlie  any  careful  investigation 
into  the  antennas,  the  major  part  of  the  waveguide  investigation  had  to 
be  completed  first.  As  pointed  out  in  last  year's  progress  section,  it 
was  anticipated  that  the  study  of  the  waveguide  properties  would  yield 
interesting  and  valuable  results,  but  that  the  calculations  would  be 
relatively  straightforward.  For  that  reason,  an  ambitious  program  was 
outlined  for  the  antenna  portions. 

It  was  discovered,  however,  that  serious  and  unexpected  conver¬ 
gence  difficulties  were  found  for  the  scattering  of  a  surface  wave 
incident  obliquely  on  a  dielectric  step.  Such  a  step  corresponds  to  the 
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side  of  a  dielectric  waveguide,  and  this  scattering  problem  is  a  basic 
element  of  the  dielectric  waveguide  investigation  when  accurate  and 
reliable  values  of  leakage  are  sought.  The  guiding  of  waves  by  a 
dielectric  strip  guide  may  be  viewed  as  the  multiple  scattering  of  a 
uniform  surface  wave  by  the  two  step  discontinuities  forming  the  side 
walls  of  the  strip  guide.  Thus,  the  scattering  of  a  surface  wave  by  a 
dielectric  step  discontinuity  is  a  basic  problem  in  the  analysis  of  such 
dielectric  waveguide  structures. 

The  previously-published  information  on  surface  waves  scattered 
by  a  dielectric  step  corresponds  to  TE  surface  waves  normally  incident 
on  the  step.  For  that  case,  convergence  is  rapid,  and  no  problems 
arise.  Our  needs  correspond  to  oblique  incident  of  either  TE  or  TM 
waves,  where  the  two  mode  types  couple  to  each  other,  and  are  no 
longer  independent,  as  occurs  for  normal  incidence.  Due  to  the  poor 
convergence  that  results,  a  very  large  number  of  modes  is  required  in 
the  mode-matching  procedure  that  one  must  use. 

Because  of  this  need  to  employ  very  many  modes  in  the  mode¬ 
matching  method  when  oblique  incidence  is  involved,  we  developed  an 
improved  calculation  procedure1  for  handling  a  very  large  number  of 
modes .  This  procedure  is  based  on  a  recognition  of  the  properties  of 
the  higher-order  modes.  The  new  method  permits  one  to  include  many 
higher-order  modes  in  the  calculation  without  requiring  a  large  com¬ 
puter  memory  core  and  a  long  computing  time.  In  addition,  the  scat¬ 
tered  mode  amplitudes  for  the  lower -order  modes  can  be  easily  deter¬ 
mined  while  including  the  effects  of  the  higher-order  modes. 

The  physics  of  the  scattering  of  surface  waves  obliquely  incident 
on  a  dielectric  step  is  of  interest  in  itself,  and  such  a  study  was  pre¬ 
viously  reported  only  briefly  by  us.2  Our  solution  at  the  time  was  only 
approximate,  but  it  contained  all  the  physics.  We  found  that  the  scat¬ 
tering  is  accompanied  by  mode  conversion  of  TE  waves  into  TM  waves, 
and  vice  versa,  so  that  an  incident  wave  gives  rise  to  reflected  and 
transmitted  waves  of  the  opposite  polarization;  by  a  Brewster  angle 
effect  for  TE  wave  incidence  (it  is  for  TM  incidence  when  uniform 
dielectric  interfaces  are  involved);  and  by  electric  field  singularities 
that  occur  at  the  corners  of  the  dielectric  step  due  to  the  TM  com¬ 
ponents  of  the  field.  An  updated  investigation  using  all  modes  was 
conducted  during  the  past  year  in  parallel  with  the  waveguide  studies, 
and  is  nearing  completion.  A  comprehensive  paper  on  this  work  is 
planned . 3 

The  first  phase  of  the  dielectric  waveguides  investigation  during 
this  past  year,  while  the  step  convergence  inquiries  were  proceeding, 
involved  examining  when  leakage  would  or  would  not  occur,  and  for 
which  modes,  on  the  insular  waveguide  and  on  the  inverted  strip  guide. 
The  criterion  used  was  developed  first  in  the  context  of  dielectric 
waveguides  for  integrated  optics  (reported  in  Work  Unit  EMl-4),and  it 
was  first  applied  there  to  rib  waveguides.  When  adapted  to  the  dielec¬ 
tric  structure  for  millimeter  waves,  it  was  necessary  to  generalize  the 
procedure.  As  a  result,  it  has  become  a  very  simple,  general  method, 
essentially  graphical,  and  based  upon  the  easily-found  dispersion  curves 
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for  the  separate  TE  and  TM  surface  waves  that  comprise  the  elements  of 
the  final  guided  wave.  A  paper  on  this  simple  and  systematic  method  is 
in  the  process  of  being  written.4 

When  reliable  recommendations  became  available  during  the  past 
year  regarding  the  convergence  properties  of  the  new  calculational 
method  for  dielectric  steps  mentioned  above,  this  procedure  was  incor¬ 
porated  into  the  computation  procedure  for  the  properties  of  dielectric 
waveguides.  Accurate  numerical  results  were  obtained  for  several 
different  cases.  It  was  found  that  the  insular  waveguide  is  usually 
more  strongly  leaky  than  the  inverted  strip  guide,  but  that  for  the 
insular  guide  the  lowest  mode  never  leaks.  In  contrast,  the  lowest 
mode  (TM)  on  the  inverted  strip  guide  can  indeed  leak,  depending  on 
parameter  values. 

In  recognition  of  the  fundamental  accomplishments  associated  with 
our  program  on  leakage  effects  on  these  open  dielectric  waveguides,  we 
were  invited  to  present  a  pair  of  papers5'6  on  this  work. 

The  measurement  program  on  leakage  from  dielectric  waveguides 
has  progressed  during  the  past  year,  but  much  difficulty  has  been 
encountered  with  spurious  radiation  effects.  The  building  of  the  set-up 
was  completed  some  time  ago,  but  unexpected  radiation  contributions 
have  distorted  the  measured  results.  It  is  clear  that  leakage  does 
occur,  and  in  the  proper  polarization,  and,  in  fact,  part  of  the  problem 
is  that  the  leakage  is  so  strong  that  multiple  reflections  of  that  energy 
are  found.  Theory  predicts  that  the  energy  should  leak  at  a  specified 
angle,  and  that,  as  one  probes  in  a  direction  perpendicular  to  the 
guide  axis,  the  leakage  field  should  increase  up  to  a  critical  distance 
and  then  drop  off  rapidly.  We  find  that  the  angle  agrees  well  quanti¬ 
tatively  with  theory,  and  that  the  transverse  variation  follows  the 
theory  qualitatively,  but  the  spurious  fields  due  to  the  feed  horn  and 
reflections  introduce  additional  puzzling  variations.  We  are  in  the 
process  of  trying  to  improve  the  set-up. 

Significant  progress  has  been  made  in  the  program  relating  to 
grating  antennas  suitable  for  millimeter  waves.  The  antenna  structure, 
which  we  are  analyzing  first,  comprises  a  rectangular  image  line  with  a 
periodic  array  of  grooves  in  its  top  face.  Due  to  the  periodicity, 
radiation  of  energy  into  space  at  some  angle  can  occur.  The  antenna  is 
analyzed  by  viewing  it  as  a  leaky  waveguide  (different  from  the  pre¬ 
vious  leaky  waveguides),  where  the  leakage  occurs  from  the  top  surface 
into  space.  We  thus  view  the  guidance  of  the  mode  in  terms  of  surface 
waves  guided  by  the  corrugated,  or  grooved,  surface,  bouncing  back 
and  forth  at  an  angle  between  the  side  walls  of  the  guide. 

Thus,  one  must  first  know  the  propagation  or  guidance  behavior  of 
a  surface  wave  traveling  along  a  corrugated  dielectric  surface  at  an 
oblique  angle  with  respect  to  the  grooves.  The  case  of  a  surface  wave 
normally  incident  on  such  a  grooved  surface  has  been  solved  previous¬ 
ly,  but  when  the  incidence  is  oblique  to  the  grooves  the  problem  is 
new,  since  now  TE-TM  mode  coupling  occurs,  and  the  problem  becomes 
a  three-dimensional  vector  one.  This  vector  problem  was  solved  by  us, 
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and  a  computer  program  based  on  that  exact  solution  was  developed, 
and  phrased  in  general  terms.  This  solution7  represents  the  first 
successful  exact  analysis  of  a  three-dimensional  vector  boundary-value 
problem  involving  a  periodic  structure. 

During  this  past  year,  we  have  attempted  to  verify  the  accuracy 
of  this  solution  by  comparing  its  predictions  with  available  experimental 
data,  in  two  different  ways.  First,  the  computer  solution  was  adapted 
to  apply  to  such  corrugated  dielectric  surfaces  for  integrated  optics, 
and  then  comparison  was  made  with  accurate  measurements  taken  in 
Germany.  The  agreement  was  excellent;  further  details  are  presented 
in  Work  Unit  EMI -4,  on  optical  waveguides. 

The  second  way  in  which  experimental  verification  is  proceeding 
involves  our  own  measurements  at  millimeter  wave  frequencies.  These 
measurements  are  not  part  of  the  original  three-year  plan,  but  they  are 
felt  to  be  worthwhile  in  view  of  the  basic  nature  of  the  oblique  guid¬ 
ance  solution  mentioned  above.  We  have  constructed  a  waveguide 
comprised  of  a  dielectric  strip  whose  top  surface  is  grooved  periodi¬ 
cally;  instead  of  the  sides  left  open,  however,  the  grooved  dielectric 
strip  is  placed  between  vertical  metallic  plates.  The  resulting  arrange¬ 
ment  is  a  trough  waveguide,  which  is  an  approximation  to  the  antenna 
itself,  but  for  which  we  can  calculate  the  propagation  and  radiation 
behavior  exactly.  In  addition,  accurate  measurements  are  possible  and 
spurious  radiation  may  be  kept  low.  We  have  already  constructed  the 
experimental  set-up,  and  we  have  made  measurements  of  the  phase 
constant  in  the  non-radiating  passband.  These  preliminary  results 
agree  well  with  theory.  Near  the  stop  band,  however,  our  early  mea¬ 
surements  show  unsatisfactory  pattern  deterioration,  but  that  region  is 
always  the  most  difficult  one. 

Many  measurements  have  been  taken  on  leaky  wave  antennas  com¬ 
prised  of  gratings  on  dielectric  image  guide,  but  no  theory  for  the 
leaky  wave  behavior  has  appeared  which  takes  the  finite  width  of  the 
guide  into  account.  We  have  recently  developed  a  theory8 ' a  which 
adapts  the  rigorous  result  mentioned  above  for  oblique  guidance  by  a 
grooved  dielectric  surface;  it  presents  a  straightforward  procedure  for 
determining  the  a  and  the  p  of  the  leaky  wave,  thus  permitting  one  to 
design  such  an  antenna  in  accordance  with  desired  radiation  charac¬ 
teristics  . 

Another  important  class  of  antennas  for  millimeter  waves  is  the 
phase-frequency  two-dimensional  steerable  arrays  described  in  our 
proposal.  Such  arrays  would  be  low-cost  substitutes  for  the  more¬ 
customary  and  more-expensive  phase-phase  scanning  planar  arrays. 
The  phase-frequency  structure  is  comprised  of  a  periodic  array  of 
leaky-wave  line  sources,  and  is  particularly  suitable  for  the  millimeter 
wave  range.  The  analysis  is  made  difficult  because  of  the  strong 
mutual  coupling  which  is  present  between  neighboring  line  sources,  and 
it  calls  for  a  new  approach  along  the  lines  of  our  proposal. 

The  overall  structure  is  decomposed  into  unit  cells  with  phase-shift 
walls,  and  then  each  unit  cell  is  viewed  transversely  and  seen  to  con- 
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sist  of  sections.  A  fundamental  section  is  that  composed  of  a  unit  cell 
in  an  infinite  medium  made  up  of  parallel  dielectric  rods,  each  rod  of 
rectangular  cross  section.  One  must  therefore  find,  as  an  intermediate 
problem,  the  dispersion  characteristics  and  mode  structure  of  the  modes 
which  propagate  in  such  a  rodded  medium. 

The  problem  of  mode  propagation  in  a  rodded  medium  of  that  type 
is  actually  a  basic,  canonical  problem  in  its  own  right,  of  interest  in  a 
number  of  antenna  applications.  In  addition  to  its  need  in  the  phase- 
frequency  class  of  dielectric  antennas,  it  is  directly  applicable,  for 
example,  to  two-dimensional  phased  arrays  which  have  dielectric  rods 
protruding  beyond  the  ground  plane  for  matching  and  bandwidth  pur¬ 
poses  . 

The  various  steps  required  for  the  solution  of  this  intermediate 
problem  were  described  in  last  year's  progress  section.  Following  those 
steps  during  this  past  year,  we  have  determined  the  properties  of  all 
the  modes  of  parallel  dielectric  slabs,  including  their  orthogonality 
properties,  and  we  developed  suitable  computer  programs  for  them.  We 
next  obtained  a  transfer  matrix  for  a  dielectric  rod  of  finite  cross 
section,  but  of  infinite  length,  and  located  in  a  unit  cell  with  phase- 
shift  walls,  thereby  accounting  completely  for  all  mutual  coupling  effects 
between  neighboring  rods  in  the  medium.  In  establishing  a  computer 
program  for  this  transfer  matrix,  however,  we  ran  into  a  difficulty 
which  we  are  in  the  process  of  overcoming. 

4.  PUBLICATIONS 

1.  S.T.  Peng  and  A. A.  Oliner,  "An  Improved  Calculation  Pro¬ 
cedure  for  Scattering  by  Dielectric  Steps,"  Proc.  Internat. 
URSI  Sympos.  on  Electromagnetic  Waves,  Paper  No.  321B, 
Munich,  Germany  (August  26-29,  1981). 

2.  J.P.  Hsu,  S.T.  Peng  and  A. A.  Oliner,  "Scattering  by  Di¬ 
electric  Step  Discontinuities  for  Obliquely  Incident  Surface 
Waves,"  Digest  of  URSI  Meeting,  p.  46,  College  Park, 
Maryland  (May  1978). 

3.  A.  Sanchez,  T.I.  Hsu,  S.T.  Peng  and  A. A.  Oliner,  "Oblique- 
Incidence  Scattering  and  Mode  Conversion  of  Surface  Waves 
by  Dielectric  Steps,"  to  be  published. 

4.  A. A.  Oliner  and  S.T.  Peng,  "A  Simple  Graphical  Method  for 
Predicting  When  Modes  Will  Leak  on  Dielectric  Waveguides  for 
Millimeter  Waves,"  to  be  published. 

5.  S.T.  Peng  and  A. A.  Oliner,  "Guidance  and  Leakage  Prop¬ 
erties  of  a  Class  of  Open  Dielectric  Waveguides,  Part  I: 
Mathematical  Formulations,"  IEEE  Trans.  Microwave  Theory 
Tech.,  Vol.  MTT-29  (Special  Issue  on  Open  Guided  Wave 
Structures),  to  be  published  in  September,  1981.  Invited 
Paper. 
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6.  A. A.  Oliner,  S.T.  Peng,  T.I.  Hsu  and  A.  Sanchez,  "Guid¬ 
ance  and  Leakage  Properties  of  a  Class  of  Open  Dielectric 

Waveguides,  Part  II:  New  Physical  Effects,"  same  as  refer¬ 
ence  5.  Invited  Paper. 

7.  S.T.  Peng,  "Oblique  Guidance  of  Surface  Waves  on  Corru¬ 
gated  Dielectric  Layers,"  Proc.  Internat.  URSI  Sympos.  on 

Electromagnetic  Waves,  Paper  No.  341B,  Munich,  Germany 
(August  26-29,  1980). 

8.  S.T.  Peng,  A. A.  Oliner  and  F.  Schwering,  "Theory  of  Di¬ 
electric  Grating  Antennas  of  Finite  Width,"  Proc.  IEEE  AP-S 
Internat.  Sympos.,  pp.  529-532,  Los  Angeles,  California 
(June  16-19,  1981). 

9.  S.T.  Peng  and  A. A.  Oliner,  "Radiation  from  Grating  An¬ 

tennas  on  Dielectric  Waveguides  of  Finite  Width,"  Proc.  11th 
European  Microwave  Conf.,  Paper  No.  B8.5,  Amsterdam,  The 
Netherlands  (September  7-11,  1981). 

5.  DoD  AND  OTHER  INTERACTIONS 

(a)  Professor  Peng  spent  two  summers  (1977,  1978)  at  the  US 
Army  CORADCOM  at  Fort  Monmouth,  New  Jersey,  working  with  Dr. 
Felix  Schwering  on  dielectric  grating  and  taper  antennas,  under  an 
LRCP  arrangement  with  the  US  Army  Research  Office.  These  inter¬ 
actions  have  continued  in  the  form  of  a  Post-LRCP  contract  and  suc¬ 
cessor  contracts. 

(b)  Both  Professors  Oliner  and  Peng,  but  particularly  Professor 
Peng,  have  held  discussions  (1979-1981)  on  electronically-scanned 
millimeter-wave  dielectric  grating  antennas  with  Dr.  Harold  Jacobs  of 
the  US  Army  ERADCOM,  Fort  Monmouth,  New  Jersey,  and  have  sup¬ 
plied  information  or  the  results  of  analyses  on  several  occasions. 

(c)  Professors  Oliner  and  Peng  participated  in  the  millimeter- wave 
workshop  sponsored  by  the  US  Army  in  October  1980  in  Colorado; 
Professor  Oliner  presented  a  short  talk  on  his  JSEP-supported  studies 
and  was  invited  to  be  a  member  of  the  panel  that  evaluated  the  state  of 
the  art. 

(d)  Professor  Oliner  has  a  contract  with  the  Air  Force's  RADC/ET 
(Hanscom  Field)  on  novel  millimeter-wave  antennas,  but  other  than  the 
grating  structures  studied  under  JSEP  sponsorship. 
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B.  HYBRID  RAY-MODE  FORMULATION  OF  DUCTED  PROPAGATION 


Professor  L.B.  Felsen 

Unit  EMI-2 

1.  OBJECTIVE(S) 

Ray  optical  and  guided  mode  expansions  have  been  used  as  alter¬ 
native  methods  for  analyzing  high-frequency  propagation  in  layers  and 
ducts.  When  the  waveguide  boundaries  are  penetrable,  the  discrete 
spectrum  of  guided  modes  is  generally  augmented  by  a  continuous 
spectrum.  In  ducts  of  large  cross  section  compared  to  the  operating 
wavelength,  the  number  of  trapped  modes  may  be  exceedingly  large  so 
that  truncation  of  the  mode  series  has  frequently  been  implemented,  and 
the  continuous  spectrum  ignored.  Alternatively,  ray  models  may  require 
a  very  large  number  of  reflected  or  refracted  rays  so  that  the  "ray 
series"  is  terminated  at  some  finite  limit.  Criteria  justifying  the  trun¬ 
cation  and  assessing  the  resulting  error  have  generally  not  been  avail¬ 
able.  Numerical  results  for  fields  generated  in  this  fashion  may  there¬ 
fore  be  open  to  question. 

In  view  of  these  observations,  it  is  desirable  to  develop  a  sys¬ 
tematic  procedure  that  can  account  for  the  remainder  of  a  truncated  ray 
or  mode  series.  A  properly-chosen,  range-dependent  combination  of 
rays  and  modes  would  not  only  furnish  an  accurate  and  convenient  field 
representation,  but  would  also  provide  new  basic  insight  into  the  propa¬ 
gation  process. 

Our  objective  is  to  develop  a  new  and  systematic  procedure  that 
can  account  for  the  remainder  of  a  truncated  ray  or  mode  series,  and 
thus  improve  the  accuracy  of  the  field  representation  for  high  fre¬ 
quency  propagation  in  layers  and  ducts.  This  is  to  be  accomplished 
with  a  properly-chosen  range-dependent  combination  of  rays  and  modes, 
which  also  should  provide  new  basic  insight  into  the  propagation  pro¬ 
cess. 

2.  APPROACH 

Field  representation  for  propagation  in  layers  and  ducts  in  the 
high  frequency  range,  where  the  field  is  generally  characterized  by 
very  many  rays  or  very  many  modes,  will  be  determined  by  means  of  a 
suitably  chosen  mixture  of  both  rays  and  modes  in  order  to  produce  a 
more  convergent  and  more  accurate  formulation. 

3.  PROGRESS 

Substantial  progress  has  been  made  on  the  numerical  calculations 
for  source-excited  waveguides  and  ducts  formed  by  two  boundaries 
confining  a  homogeneous  or  transversely  inhomogeneous  medium.  The 
hybrid  ray-mode  theory  has  been  compared  with  reference  solutions 
based  on  a  normal  mode  (discrete  and  continuous  spectrum)  represen¬ 
tation.  Calculations  have  been  performed  for  the  simplest  case  of  a 
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homogeneously  filled  parallel  plane  waveguide  with  impedance  walls.1 
They  have  also  been  done  for  open  waveguides  comprised,  respectively, 
of  a  transversely  homogeneous  and  of  a  transversely  inhomogeneous 
layer  above  a  semi-infinite  homogeneous  bottom.  The  former  configura¬ 
tion  was  explored  in  connection  with  seismic  propagation,  and  required 
extension  of  the  hybrid  method  to  transient  fields.2  The  latter  con¬ 
figuration  is  of  interest  for  sound  propagation  in  an  ocean  channel  with 
penetrable  bottom.3  Here,  the  inhomogeneous  refractive  index  was 
assumed  to  vary  exponentially  with  depth;  as  noted  previously,  this 
model  (by  a  simple  coordinate  mapping)  also  describes  propagation  along 
a  cylindrically  curved  concave  boundary  in  a  uniform  medium.4 

In  all  of  the  cases  studied,  the  hybrid  ray-mode  formulation  has 
exhibited  the  versatility  and  accuracy  found  in  previous  investigations 
of  simpler  problems.  All  of  the  problems  are  relevant  to  electromagnetic 
wave  propagation  in  various  environments  (ground  wave,  tropospheric, 
optical  graded  index,  circular  waveguide)  although  some  of  them  have 
been  phrased  in  the  context  of  wave  propagation  in  other  (acoustic  and 
elastic)  media.  Thus,  the  hybrid  method  is  now  established  as  an 
effective  tool  for  analyzing  in  a  novel,  numerically  efficient  and  physi¬ 
cally  appealing  manner  source-excited  propagation  in  multimode  guiding 
regions.  This  assertion  is  supported  by  the  various  publications  that 
have  been  prepared  on  this  subject  area.1-11 

The  hybrid  method  has  far  reaching  implications  for  relatively 
arbitrary  propagation  in  guiding  environments,  that  can  now  be  ex¬ 
plored  within  a  general  framework.  This  important  problem  area  is 
addressed  in  the  new  proposal. 

4.  PUBLICATIONS 

1.  L.B.  Felsen  and  A.  Kamel,  "Hybrid  Ray-Mode  Formulation  of 
Parallel  Plane  Waveguide  Green's  Functions,"  IEEE  Trans,  on 
Antennas  &  Propag.,  AP-29,  pp.  637-649  (1981). 

2.  A.  Kamel  and  L.B.  Felsen,  "Hybrid  Ray-Mode  Formulation  of 
SH  Motion  in  a  Two-Layer  Half  Space,"  accepted  for  publica¬ 
tion  in  Bull.  Seismol.  Soc.  Am. 

3.  E.  Niver,  S.H.  Cho  and  L.B.  Felsen,  "Rays  and  Modes  in  an 
Acoustic  Channel  with  Exponential  Velocity  Profile,"  accepted 
for  publication  in  Radio  Science. 

4.  L.B.  Felsen  and  T.  Ishihara,  "Hybrid  Ray-Mode  Formulation 
of  Ducted  Propagation,"  J.  Acoust.  Soc.  Am.,  (March  1979). 

5.  T.  Ishihara,  L.B.  Felsen  and  A.  Green,  "High  Frequency 
Fields  Excited  by  a  Line  Source  Located  on  a  Perfectly  Con¬ 
ducting  Concave  Cylindrical  Surface,"  IEEE  Trans,  on  An¬ 
tennas  and  Propagation,  AP-26,  No.  6,  (November  1978). 

6.  L.B.  Felsen  and  T.  Ishihara,  "High-Frequency  Surface  Fields 
Excited  by  a  Point  Source  on  a  Concave  Perfectly  Conducting 
Cylindrical  Boundary,"  Radio  Science,  24,  No.  2,  pp.  205- 
216,  (March- April  1979). 
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7.  T.  Ishihara  and  L.B.  Felsen,  "High  Frequency  Fields  Excited 
by  a  Line  Source  Located  on  a  Concave  Cylindrical  Impedance 
Surface,"  IEEE  Trans,  on  Antennas  and  Propag.,  AP-27,  No. 
2,  (March  1979). 

8.  S.H.  Cho,  C.G.  Migliora  and  L.B.  Felsen,  "Hybrid  Ray-Mode 
Formulation  of  Tropospheric  Propagation,"  AGARD  Symposium 
on  Special  Topics  in  HF  Propagation,  Lisbon,  Portugal,  p.  11, 
1-8,  15  (May  28- June  1,  1979). 

9.  L.B.  Felsen,  "Hybrid  Ray-Mode  Fields  in  Inhomogeneous 
Waveguides  and  Ducts,"  J.  Acoust.  Soc.  Am.,  69(2),  pp. 
352-361  (February  1981). 

10.  E.  Topuz  and  L.B.  Felsen,  "High-Frequency  Electromagnetic 
Fields  on  Perfectly  Conducting  Concave  Cylindrical  Surfaces," 
IEEE  Trans,  on  Antennas  and  Propag.,  AP-28,  No.  6, 
(November  1980). 

11.  E.  Topuz,  E.  Niver  and  L.B.  Felsen,  "Electromagnetic  Fields 
Near  a  Concave  Perfectly  Conducting  Cylindrical  Surface," 
accepted  for  publication  in  IEEE  Trans,  on  Antennas  and 
Propag . 

5.  DoD  AND  OTHER  INTERACTIONS 

(a)  Two  seminars  on  hybrid  ray-mode  methods  for  concave  sur¬ 
faces  presented  at  ARO  (Durham)  1980. 

(b)  Three  seminars  presented  during  1980-81  at  U.S.  Geological 
Survey,  Menlo  Park,  California,  on  hybrid  methods  applied  to  seis¬ 
mology  . 

(c)  Collaboration  with  Dr.  G.  Jacobsen,  Technical  Univ.  of  Den¬ 
mark,  Lyngly,  Denmark,  on  hybrid  methods  applied  to  beam  propagation 
in  optical  fibers. 

(d)  Collaboration  with  Dr.  J.M.  Arnold,  Univ.  of  Nottingham, 
Nottingham,  England,  on  hybrid  methods  applied  to  tapered  dielectric 
waveguide . 

(e)  Collaboration  with  Professor  A.M.  Scheggi,  Institute  for 
Electromagnetic  Wave  Propagation,  Florence,  Italy,  on  hybrid  methods 
applied  to  optical  fibers. 

(f)  Collaboration  with  NATO  Underwater  Sound  Laboratory,  La 
Spezia,  Italy,  on  hybrid  methods  for  propagation  in  ocean  channels. 

(g)  Collaboration  with  Professor  E.  Topuz,  Department  of  Elec¬ 
trical  Engineering,  Technical  University  of  Instanbul,  Turkey,  on 
application  of  hybrid  methods  to  imaging  in  multi-mode  optical  wave¬ 
guides  . 
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C.  WAVE-MATTER  INTERACTIONS 


Professors  S.  Barone,  S.  Gross  and  N.  Marcuvitz 
Unit  EMI -3 

1.  OBJECTIVE(S) 

Our  objective  is  to  extend  tc  nonlinear  and/or  higher  power  levels 
and  to  develop ,  both  theoretically  and  experimentally ,  engineering 
expressions  for  "collision"  (interaction)  terms  necessary  in  the  descrip¬ 
tion  of  particle-particle,  particle-wave  and  wave-wave  interaction  prob¬ 
lems  arising  in  electromagnetic  v/ave  propagation  through  media  in  which 
ionization  and  related  phenomena  arise .  The  latter  include  microwave 
breakdown  processes  associate,  with  high  power  wave  propagation 
through  the  atmosphere,  -rucrowave  propagation  through  and  modifica¬ 
tion  of  the  ionosphere,  plasma  screening  effects  in  laser  wave  absorp¬ 
tion  and  reflection  by  materials,  etc. 

2.  APPROACH 


An  overall  analytical  theme  is  to  obtain  a  kinetic  basis  for  models 
of  "collisional"  interactions.  Kinetic  models  must  be  consistent  with  the 
rate  equations  and  the  experimental  rate  coefficients  and  must  provide 
reasonably  correct  representations  of  physical  phenomena  both  at  the 
kinetic  and  fluid  levels.  For  example,  a  possible  general  form  for  the 
kinetic  equation  describing  the  density  distribution  function  f,(v,r,t)  of 

a 

particles  of  type  "a"  from  which  fluid  equations  may  be  derived  is: 


d 

3t 


D 


Va 


+  v. 


A  fa  +  Ca(f1,  f2,...) 


where  D  and  A  are  diffusion  and  friction  coefficients  indicative  of  long 
range  (particle-wave)  collisions  and  C  is  indicative  of  short  range 

a 

collisions.  While  it  is  also  possible  to  include  some  short  range  (particle- 
particle)  collisional  contributions  in  D  and  A,  others  involving  transfer 
of  non-infinitesimal  amounts  of  energy  and  non- particle  conserving 
collisions,  such  as  ionization  and  recombination,  must  be  included  in  the 
Ca  term.  Our  initial  collisional  models  are  suitable  only  for  kinetic 

equations  that  are  isotropic  in  velocity  or  energy.  However,  in  the 
presence  of  large  applied  electric  field  the  isotropy  assumption  has 
limited  validity  particularly  when  electron  velocities  are  larger  than  the 
thermal  speed .  The  use  of  the  microscopic  kinetic  approach  of 
Klimontovitch ,  together  with  the  renormalization  and  quasiparticle  tech¬ 
niques,  are  being  explored  to  remove  the  isotropy  limitation.  Moment 
or  fluid  equations  derived  from  such  kinetic  models  will  be  checked  for 
consistency  with  experimentally  determined  momentum,  energy  transfer, 
ionization,  etc.,  rates.  The  sensitivity  of  kinetic  and/or  fluid  descrip¬ 
tions  to  various  collisional  models  will  be  studied  by  developing  and 
comparing  analytical  and  numerical  models  of  overall  physical  processes 
on  our  interactive  graphic  computer  facility. 
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3.  PROGRESS 

A .  Introductory  Remarks 

General  approaches  have  been  developed  for  "collision"  (interac¬ 
tion)  terms  required  in  the  description  of  particle- particle,  particle- 
wave  and  wave-wave  interaction  processes  which  arise  at  those  high- 
power  levels  at  which  ionization  and  related  phenomena  take  place 
during  the  propagation  of  an  electromagnetic  wave  through  various 
media.  The  common  underlying  analytical  theme  has  been  to  obtain  a 
kinetic  base  for  models  of  "collisional"  interactions.  Although  the 
research  has  been  applied  mostly  at  laser  frequencies,  a  general  ap¬ 
proach  in  terms  of  a  quasiparticle  treatment  of  wave  propagation 
through  nonlinear  media1-3  is  also  applicable  at  microwave  frequen¬ 
cies.  We  have  applied  this  approach  to  specific  electromagnetic  wave- 
media  interactions  such  as  ionospheric  scintillation  and  turbulence,4-8 
laser  generation  and  propagation,9-12  laser  interaction  with  metallic 
surfaces13-19  and,  to  a  lesser  extent,  high-power  microwave  propaga¬ 
tion  and  breakdown.20'21  Investigation  of  the  basic  interactions  in 
different  electromagnetic  frequency  regimes  has  been  helpful  in  the 
development  of  both  analytical  treatments  and  computer  models  for 
comparison  with  experiment. 

B .  Ionospheric  Turbulence 


A  three-dimensional  extension  of  our  previously  reported  two-fluid 
theory  of  ionospheric  plasma  bubbles  in  the  collision-dominated  regime 
was  completed  in  the  past  period.  Effects  of  magnetic  dip  and  the  finite 
extent  of  the  bubble  along  the  magnetic  field  have  been  included.  The 
theory  clarifies  the  nature  of  the  lifting  force  on  a  bubble.  For  closed 
bubbles  which  do  not  deplete  an  entire  flux  tube,  two-dimensional 
theories  are  valid  for  bubbles  having  the  smaller  of  observed  transverse 
dimensions.  For  closed  bubbles  with  the  larger  of  observed  transverse 
dimensions,  or  for  smaller  bubbles  that  do  not  extend  most  of  the 
length  of  a  flux  tube,  three-dimensional  effects  become  significant.  We 
find  (1)  bubbles  rise  more  slowly  than  predicted  by  two-dimensional 
theories,  (2)  bubbles  with  different  ratios  of  longitudinal  to  transverse 
dimensions  rise  at  different  rates,  (3)  bubbles  which  are  not  sym¬ 
metrical  about  a  meridian  plane  drift  east- west.  This  is  also  a  feature 
of  two-dimensional  theory  which  to  our  knowledge  has  not  been  pre¬ 
viously  noted. 

Efforts  to  develop  renormalization  theory  in  a  form  suitable  for  the 
practical  solution  of  ionospheric  turbulence  problems  were  continued. 
All  known  techniques  can  be  viewed  as  approximations  that  strive  to 
evolve  a  precise  renormalization  theory  valid  for  strong  (nonlinear) 
turbulence.  Our  work  developed  along  relatively  conventional  lines, 
and  also  in  directions  indicated  by  renormalization  group  and  divergent 
series  techniques.  During  the  past  year,  previous  general  results5 
have  been  related  to  the  direct  interaction  approximation  in  such  a  way 
that  systematic  improvements  may  be  possible.  This  is  one  approach 
towards  correcting  some  practical  incompleteness  in  the  direct  inter¬ 
action  approximation  for  calculating  the  spectral  index.  We  have  made 
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some  progress  in  carrying  out  this  improved  approximation  in  simple 
turbulent  systems  and  in  systems  driven  by  stochastic  fields. 

One  of  the  uncertainties  in  turbulence  calculations  revolves  around 
the  accuracy  of  the  first  term  in  renormalized  perturbation  expansions 
that  are  used  to  represent  the  statistical  properties  of  waves  propagat¬ 
ing  through  the  ionosphere.  The  accuracy  is  certainly  in  question  for 
strongly  turbulent  ionospheres  wherein  the  medium  statistics  impose 
radically  different  statistics  on  waves  propagating  through  the  medium. 
We  have  investigated  the  accuracy  of  a  number  of  renormalization  and 
summable  divergent  series  procedures  for  expansions  arising  from 
different  medium  statistics  and  levels  of  turbulence.  No  definitive 
results  have  as  yet  been  achieved.  This  difficult  problem,  which  ap¬ 
pears  in  a  different  context  in  the  ionospheric  turbulence  work  above, 
is  generic  to  strong  turbulence  theory  and  has  only  been  partially 
solved . 

C.  Laser-Metal  Interactions 


Some  of  the  basic  phenomena  which  arise  in  the  interaction  of 
electromagnetic  fields  with  metal  surfaces  have  been  examined.  The 
major  emphasis  has  been  on  the  phenomena  which  occur  during  the 
formation  of  the  hot  layer  of  ionized  gas  near  the  metal  surface.16  The 
relationship  between  the  growth  of  the  laser-generated  discharge  and 
the  time  dependence  of  the  plasma  absorptivity  was  of  particular  inter¬ 
est  because  of  the  possibility  of  experimental  verification.  During  the 
time  period  that  the  pulse  power  is  near  its  maximum  value,  a  large 
variation  takes  place  in  the  density  of  the  ablated  material;  thus  the 
effect  of  thermal  diffusion  (via  electrons)  is  important  and  is  assumed 
to  dominate  the  effect  of  hydrodynamic  expansion  in  determining  cooling 
rates.  The  interaction  process  proceeds  in  a  sequence  of  stages  in 
each  of  which  different  phenomena  take  on  primary  importance.  The 
various  stages  were  isolated  in  our  formulation  for  the  purpose  of 
analysis.  These  were  (1)  heating  of  the  metal  and  formation  of  blow-off 
vapor,  (2)  heating  of  the  vapor  leading  to  breakdown,  (3)  propagation 
of  the  breakdown  front  and  (4)  build-up  of  higher  stages  of  ionization 
in  the  already  formed  plasma. 

In  the  analysis  of  the  metal  heating  stage  we  presented  a  method 
of  extrapolating  existing  data  on  material  properties  into  the  tempera¬ 
ture  range  of  interest.  Because  of  the  sensitive  nature  of  the  depen¬ 
dence  of  vapor  pressure  on  temperature  the  results  based  on  the  ex¬ 
trapolated  values  differed  significantly  from  those  obtained  by  Harrach 
using  a  more  approximate  expression.22  A  numerical  integration  of  the 
nonlinear  heat  transfer  equation  led  to  surface  temperature  values 
substantially  different  from  the  predictions  of  the  steady  state  model  of 
surface  evaporation. 

The  initial  stages  of  the  growth  toward  avalanche  breakdown  in  the 
vapor  proceeded  in  a  way  predicted  previously  when  the  spatial  varia¬ 
tion  of  the  blow-off  vapor  density  and  pressure  were  ignored.22  How¬ 
ever,  the  time  required  for  breakdown  was  substantially  increased 
(doubled)  over  what  was  observed  in  the  absence  of  the  thermal  dif- 
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fusion  terms.  This  term  played  a  critical  role  during  the  avalanche 
process  and  during  the  subsequent  propagation  of  the  discharge  front. 

The  propagation  speed  of  the  discharge  front,  in  our  model,  was 
determined  by  the  effectiveness  of  the  thermal  diffusion  mechanism  in 
equilibrating  the  electron  temperature  within,  and  in  layers  adjacent  to, 
the  discharge  region.  Thus  the  energy  absorbed  within  most  of  the 
plasma  goes  toward  the  production  of  ionization  in  a  spatial  region  in 
advance  of  the  discharge  front.  During  the  propagation  period,  lasting 
several  nanoseconds ,  the  net  transmission  of  the  plasma  decreases 
sharply  as  the  absorption  optical  path  length  increases.  The  absorbed 
energy  is  stored  in  the  ionized  gas  and  is  available  for  subsequent 
heating  of  the  metal,  possibly  through  emission  of  uv  radiation. 

Experiments  by  W.T.  Walter,  et  al18  have  been  performed  under 
conditions  which  correspond  well  to  the  assumptions  of  our  model. 
Direct  real-time  reflectance  measurements  were  made  on  copper  targets 
irradiated  with  30  ns  pulses  from  a  ruby  laser.  When  the  pulse  inten¬ 
sity  was  around  5  x  108  watts/cm2  the  observed  reflectivity  variation 
was  in  quantitative  agreement  with  our  results.  Our  studies  also  led  to 
predictions  on  plasma  temperature  and  state  of  excitation.  These  will 
be  subject  to  experimental  test  when  more  complete  measurements  of 
neutral  and  ionized  Cu  lines  are  made  in.  the  "light  flash"  radiated  by 
the  plasma  in  front  of  the  metal  surface. 

D.  High-Power  Microwave- Atmosphere  Interactions 

The  original  model  developed  for  the  propagation  of  a  high-power 
pulse  has  been  expanded  to  include  the  temperature  of  the  atmospheric 
gases  and  to  allow  for  their  heating  by  the  propagating  pulse.  As  we 
move  from  the  shorter  nanosecond  pulses  of  interest  for  sneak-through 
propagation  to  longer  microsecond  pulses  of  interest  for  a  dump  of 
energy  into  a  breakdown  plasma,  gas  heating  and  hydrodynamic  effects 
will  play  important  roles. 

The  propagation  of  a  10-ns,  3-GHz  pulse  through  5000  feet  at  sea 
level,  as  obtained  from  our  computer  model,  is  illustrated  in  Figure  1. 
The  peak  power  of  7.3  MW/cm2  is  5.3  times  the  power  density  required 
for  dc  or  static  breakdown  (i.e.,  POW  =  5.3).  Pulse  shapes  at  the 
antenna,  one-quarter,  one- half,  three-quarters  and  full  range  are 
displayed  by  the  solid  curves  along  the  z-axis  in  the  upper  portion  of 
Figure  1.  Dotted  curves  show  the  time  dependence  of  the  electron 
temperature  at  each  range  and  the  additional  dotted  curve  at  full  range 
displays  the  time  dependence  of  the  electron  density.  The  erosion  of 
the  tail  during  propagation  of  the  pulse  is  more  clearly  illustrated  in 
the  lower  portion  of  Fig.  1  where  the  five  solid  curves  in  the  upper 
portion  are  displayed  on  a  single  xy  plot  (i.e.,  the  outer  solid  curve  is 
the  pulse  shape  at  z  =  0  and  the  successively  narrower  solid  curves  are 
pulse  shapes  at  quarter  range  intervals).  At  full  range  of  5000  feet 
only  23.8%  of  the  initial  pulse  is  transmitted  and  an  electron  density  of 
10*  cm"3  is  produced. 
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Fig.  1  Propagation  of  a  3-GHz  microwave  pulse  for  5000  feet  through 
the  atmosphere . 
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To  validate  the  microwave  propagation  computer  model,  we  have 
assembled  a  conventional  S-band  rectangular  waveguide  system  at  the 
Polytechnic.  A  section  of  waveguide  between  two  teflon  windows  is 
connected  to  a  vacuum  pump  so  that  gas  composition  and  pressure  can 
be  varied  as  desired.  Air  breakdown  and  pulse  shortening  are  ob¬ 
served  when  the  pressure  within  this  waveguide  section  is  reduced  from 
atmospheric  to  a  region  in  the  vicinity  of  the  Paschen  minimum.  In 
Fig.  3,  for  example,  pulse  shortening  of  a  30-kW,  80-ns,  3-GHz  pulse 
is  shown  to  take  place  when  the  pressure  is  between  1  and  4  torr. 

Each  experimental  picture  in  Fig.  2  is  a  double  exposure.  One 
exposure  shows  the  pulse  transmitted  in  the  partially  evacuated  wave¬ 
guide  section  when  the  waveguide  pressure  is  reduced  to  the  values 
specified  under  the  picture.  The  second  exposure  shows  the  pulses 
transmitted  when  the  waveguide  contains  atmospheric  pressure.  At  760 
torr  the  field  strength  is  too  low  to  increase  the  electron  density  to  a 
level  where  any  modification  of  the  pulse  shape  occurs.  Also  shown  in 
Fig.  2  are  the  computer  model  results  for  propagation  of  a  pulse  at  the 
specified  pressures.  Similar  to  the  displays  in  Fig.  1,  the  solid  curves 
show  the  pulse  shape  incident  at  the  teflon  window  as  well  as  after 
propagating  10,  20,  30  and  40  feet  in  the  waveguide  section  with  re¬ 
duced  pressure.  In  the  results  shown  in  Fig.  2,  the  actual  experi¬ 
mental  pulse  shape  was  digitized  into  the  computer  and  used  in  the 
model . 

In  Fig.  3  a  similar  comparison  is  shown  for  a  50  kW  400-ns,  3-GHz 
pulse.  Pulse  shortening  is  observed  within  a  pressure  region  of  2  to 
21  torr  for  this  longer  and  higher  peak  pulse.  Above  2  torr  the  agree¬ 
ment  between  experiment  and  computer  model  is  good;  however,  below  2 
torr  the  model  results  indicate  erosion  of  the  pulse  tail,  which  is  not 
observed  experimentally.  The  comparison  in  both  Figs.  2  and  3  has 
been  improved  by  modifying  the  empirical  coefficients  used  in  the  ioni¬ 
zation  coefficient  Hj/  the  drive  velocity  v^  and  the  momentum  transfer 

frequency  u.  In  effect  has  been  divided  by  1.3,  u  has  been  divided 

by  1.2  and  v^  has  been  multiplied  by  1.2.  Further  changes  in  the 

empirical  coefficients  are  anticipated.  Two  additional  processes  must  be 
considered  in  completing  the  validation:  (1)  the  initial  electron  density 
to  be  used  in  the  computer  model  and  any  residual-charged  or  easily- 
ionizable  species  remaining  from  a  previous  pulse,  and  (2)  heating  of 
the  waveguide  window  producing  initiators  of  breakdown.  The  experi¬ 
mental  data  in  Fig.  2  were  taken  at  a  pulse  repetition  rate  of  357;  in 
Fig.  3  the  rate  was  reduced  to  8  pps.  The  initial  electron  density  was 
100  cm-3  in  the  computer  model  results  displayed  in  Figures  1-3. 
Substantial  reductions  (e.g.,  to  10,  1  or  0.1  initial  electrons  cm-3)  do 
produce  significant  changes  in  the  shape  of  the  transmitted  pulse.  It 
is  evident  that  additional  work  both  theoretically  and  experimentally  is 
required  to  complete  the  validation  of  the  computer  model. 
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Fig.  2b  Comparison  of  experiment  with  computer  model  results  for 
propagation  of  a  3-GFte,  30- kW  80-ns  microwave  pulse  in  a 
partially-evacuated  waveguide . 
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3a  Comparison  of  experiment  with  computer  model  results  for 
propagation  of  a  3-GHz,  50- kW  400-ns  microwave  pulse  in  a 
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Fig.  3b  Comparison  of  experiment  with  computer  model  results  for 
propagation  of  a  3-GHz,  50-kW  400-ns  microwave  pulse  in  a 
partially  evacuated  waveguide. 
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5.  DoD  AND  OTHER  INTERACTIONS 

About  three  years  ago.  Professor  Marcuvitz  was  contacted  by  R. 
Wasneski  of  NAV  AIR  about  our  willingness  to  investigate  the  possible 
realization  of  atmospheric  channels  created  by  and  capable  of  propagat¬ 
ing  very-high-power  microwaves  for  distances  of  a  few  kilometers. 
Such  a  possibility  had  been  reported  in  a  Russian  publication.  A  quick 
NAV  AIR  sponsored  study  indicated  that  channel  creation  was  not 
feasible  with  currently  available  microwave  sources.  Other  possibilities, 
however,  for  the  transmission  of  high-power  microwave  pulses  were 
suggested  during  the  study  and  have  led  to  a  continuing  relationship 
with  NAV  AIR  and  NRL. 
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As  a  related  consequence  of  this  program,  Professor  Marcuvltz  was 
invited  to  become  a  member  of  the  recent  special  DoD  committee  to 
undertake  the  development  of  a  preliminary  plan  for  a  U.S.  High-Power 
Microwave  Technology  (HPMT)  program  with  potential  applications  to 
DoD  missions. 

The  background  for  our  quick  response  stemmed  from  some  of  the 
basic  nonlinear  and  turbulent  wave  propagation  studies  that  originated 
from  our  JSEP  research.  More  specifically,  our  program  on  this  work 
unit  provided  basic  technical  contributions  for  our  recent  NRL- 
sponsored  program  on  "High  Power  Microwave  Pulse  Propagation,"  which 
is  both  theoretical  and  experimental .  The  computer  models  resulting 
from  this  analysis  and  their  low  power  experimental  verification  at  MRI 
by  Professor  Walter  are  going  to  be  checked  with  the  gigawatt  magne¬ 
tron  facility  under  construction  at  NRL.  Modification  of  transmitted 
microwave  pulse  shapes,  which  have  been  examined  at  kW/cm2  power 
densities  at  Polytechnic,  can  be  tested  at  MW/cm2  peak  power  densities 
at  NRL  and  used  to  validate  the  computer  model  for  microwave  pulse 
propagation  being  developed  at  the  Polytechnic.  This  is  a  cooperative 
project  in  which  Polytechnic  personnel  and  equipment  will  be  moved 
temporarily  to  Washington,  DC,  to  carry  out  the  measurements  together 
with  the  NRL  personnel. 
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D.  OPTICAL  WAVEGUIDING  STRUCTURES 


Professors  T.  Tamir,  A. A.  Oliner  and  L.B.  Felsen 

Unit  EMI-4 
1.  OBJECTIVE(S) 

(a)  To  investigate  the  waveguiding  and  coupling  properties  of 
dielectric  gratings  having  asymmetric  (blazed)  profiles  for  applications 
to  beam  couplers,  filters,  distributed-feedback  lasers  and  other  inte¬ 
grated-optics  components.  In  particular,  criteria  for  the  design  of 
high-efficiency  gratings  with  prescribed  desirable  characteristics  will  be 
studied . 

(b)  To  analyze  and  characterize  the  behavior  of  basic  discontin¬ 
uity  structures  found  in  integrated-optical  circuits,  to  determine  their 
reflection  and  scattering  characteristics.  No  one  has  yet  performed  such 
analyses  to  any  degree  of  accuracy.  The  simple  dielectric  step  junction 
will  be  examined  first.  It  will  be  analyzed  rigorously,  and  then  to 
varying  degrees  of  approximation,  determining  the  simplest  form  suit¬ 
able  for  recommendation. 

(c)  To  analyze  the  propagation  characteristics  of  various  linear 
(strip-type)  waveguides  of  importance  in  integrated-optical  circuitry. 
Most  such  waveguides  have  so  far  been  analyzed  by  others  only  to  first 
order.  A  better  analysis  has  already  revealed  the  presence  in  some 
waveguides  of  leakage  and  resonance  effects  which  were  heretofore 
unrecognized  and  which  could  have  important  device  implications . 
Several  different  types  of  waveguide  will  be  examined,  and  to  various 
degrees  of  accuracy.  The  results  obtained  under  (b)  above  for  the 
dielectric  step  junction  (since  the  sides  of  strip  waveguides  correspond 
to  such  junctions)  will  be  applied  in  this  connection. 

(d)  To  examine  analytically  the  properties  and  the  implications  of 
a  newly-recognized  total  absorption  effect  associated  with  thin  lossy 
layers  and  with  metallic  periodic  structures  which  are  not  perfectly 
conducting.  Recent  work  has  shown1-3  that,  if  electromagnetic  energy 
is  incident  on  a  planar  structure  consisting  of  a  thin  lossy  layer  and/or 
a  metallic  grating,  all  or  nearly  all  of  the  incident  energy  can  be  ab¬ 
sorbed  by  properly  adjusting  the  incidence  angle  and  the  dimensions  of 
the  planar  structure.  The  interesting  aspect  of  this  effect  is  that  total 
absorption  can  be  realized  without  the  need  for  thick  lossy  absorbers. 
The  aim  of  the  proposed  work  is  to  first  examine  the  total-absorption 
effect  in  detail  for  both  layered  media  and  periodic  (grating)  structures 
and  then  to  apply  the  results  to  a  variety  of  applications,  which  include 
the  design  of  PN  photodiodes  having  both  high  quantum  efficiencies  and 
very  rapid  response  characteristics. 

(e)  To  explore  a  new  asymptotic  theory  for  propagation  in  guid¬ 
ing  regions  with  inhomogeneous  permittivity  properties.  The  results 
are  expected  to  be  relevant  for  integrated-optical  and  optical  fiber 
waveguides,  underwater  sound  wave  propagation  in  the  ocean,  ducting 
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in  the  ionosphere,  and  related  applications.  This  theory  will  be  applied 
to  various  graded-index  waveguides,  including  both  slabs  and  fibers, 
with  a  view  toward  accommodating  analytical  refractive  index  profiles  of 
arbitrary  shape.  Such  structures  would  include  profiles  with  discon¬ 
tinuous  changes  in  the  index  or  its  derivatives;  the  theory  will  later  be 
generalized  to  waveguides  with  asymmetrical  boundaries,  with  longitu¬ 
dinal  variation,  and  with  curved  axes. 

(f)  To  examine  the  mathematical  and  physical  implications  of 
surface  waves  guided  obliquely  on  periodic  planar  dielectric  waveguides. 
When  surface  waves  are  guided  perpendicularly  to  the  periodic  grooves, 
the  mathematical  problem  is  a  scalar  two-dimensional  one  in  which  TE 
and  TM  waves  remain  independent  of  each  other.  For  oblique  guidance, 
TE-TM  mode  coupling  occurs  and  the  problem  becomes  a  vector  three- 
dimensional  one.  This  mode  coupling  results  in  a  variety  of  physical 
effects  that  should  be  explored. 

2.  APPROACH 

The  approaches  used  in  the  various  aspects  of  this  program  employ 
a  number  of  theoretical  techniques  that  have  already  been  developed  to 
a  high  degree  of  sophistication  and  success  in  the  area  of  microwave 
engineering.  In  particular,  the  researchers  on  this  program  have 
personally  contributed  significantly  to  the  development  of  these  tech¬ 
niques,  which  include  transverse  resonance  methods,  various  approxi¬ 
mate  and  rigorous  modal  methods,  and  periodic  structure  theory.  In 
most  aspects  of  the  program,  the  theoretical  techniques  themselves  are 
developed  further  in  the  course  of  the  investigations.  This  is  espec¬ 
ially  true  in  connection  with  the  investigations  relating  to  graded-index 
guiding  structures,  where  a  new  method  embodying  local  evanescent 
fields  with  complex  phase  is  employed,  instead  of  the  conventional 
asymptotic  theories  which  utilize  local  ray  fields  or  mode  fields  with  real 
phase. 

3.  PROGRESS 

A  Dielectric  Gratings  for  Integrated  Optics 

Our  present  study  of  dielectric  gratings  for  integrated-optics 
applications  has  continued  our  prior  work  on  periodic  structures,  which 
had  involved  the  development  of  exact  field  solutions  and  their  accurate 
numerical  evaluation,1  in  general,  and  their  application  to  the  analysis 
and  evaluation  of  beam  couplers,3-5  in  particular.  During  the  earlier 
portion  of  our  current  effort,  we  had  carried  out  the  following  studies: 

(1)  We  developed  approximate6-9  and  rigorous12' 13  methods  to 
deal  with  scattering  by  gratings  and,  in  particular,  to  confirm  and 
assess  the  desirable  directional  effects  that  can  be  achieved  by  suitably 
blazed  gratings. 

(2)  We  derived7'8  simple  criteria  for  the  construction  of  blazed 
gratings  with  nearly  100%  coupling  efficiency. 
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(3)  We  verified10' 11  some  of  the  above  theoretical  results  by  an 
experimental  study  which  used  microwave  models  of  optical  gratings. 

The  past  year  has  been  used  to  consolidate  and  expand  those 
results,  in  general,  and  to  assess  the  effects  of  certain  physical  grat¬ 
ing  parameters,  in  particular.  Specifically,  these  efforts  have  included 
the  following  aspects. 

(4)  We  have  completed  the  analysis  dealing  with  the  scattering 
and  guiding  of  waves  by  gratings  with  arbitrary  profiles. 18,19  In 
particular,  we  have  extended  this  rigorous  analysis  to  include  the  effect 
of  absorption  losses. 

(5)  We  have  extended  the  results  for  preferential  guidance  by 
blazed  gratings  to  profiles  other  than  triangular. 19,20  This  particular 
effort  has  verified  that  the  strong  blazing  effect  can  be  maintained  over 
a  relatively  wide  (±  20%)  frequency  spectrum,  and  for  relatively  large 
changes  in  grating  profile  and/or  in  the  other  physical  parameters. 
Thus,  the  effects  due  to  fabrication  tolerances  or  to  grating  imperfec¬ 
tions  appear  to  be  minor.20 

(6)  We  have  obtained  further  verification  of  our  theoretical  re¬ 
sults  by  means  of  the  experimental  facet  using  microwave  models.  In 
this  context,  we  have  extended  the  modeling  approach  to  a  wider  class 
of  problems.20 

In  summary,  the  effort  dealing  with  dielectric  gratings  for  beam¬ 
coupling  applications  has  achieved  results  that  are  fully  consistent  with 
the  objectives  of  the  originally  proposed  investigation. 

B.  Scattering  Properties  of  the  Basic  Dielectric  Step  Junction 

We  have  been  examining  both  approximate  and  highly-accurate 
analytical  formulations  for  the  dielectric  step  junction.  That  junction  is 
of  importance  because  it  corresponds  to  the  side  wall  of  dielectric  strip 
waveguides,  and  also  because  it  occurs  as  an  integral  element  of  other 
structures  for  integrated  optics  applications,  such  as  multiplexers  and 
modulators . 

With  respect  to  the  accurate  formulation,  we  ran  into  the  conver¬ 
gence  difficulties  encountered  in  the  millimeter  wave  studies,  and  des¬ 
cribed  under  Work  Unit  EMI-1.  The  problems  are  analogous,  since  in 
optics  the  dielectric  layer  is  placed  on  a  dielectric  substrate,  whereas 
at  millimeter  waves  one  uses  a  metal  ground  plane  instead  of  a  sub¬ 
strate.  The  calculational  details  are  different  (as  is  the  computer  pro¬ 
gram),  but  the  basic  features  are  the  same.  We  therefore  used  the  new 
calculational  method  discussed  in  Work  Unit  EMI-1,  and  we  obtained 
results  with  greatly  improved  convergence.  With  the  computer  program 
thus  modified  and  improved,  we  are  obtaining  various  numerical  results 
for  the  scattering  and  mode  conversion  of  surface  waves  occurring  upon 
their  oblique  incidence  on  dielectric  steps.  These  calculations,  anal¬ 
ogous  to  those  described  under  Work  Unit  EMI-1,  are  nearing  com¬ 
pletion  . 
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C.  Leakage  Characteristics  of  Dielectric  Strip  Waveguides 

Our  studies  on  the  propagation  properties  of  three-dimensional 
dielectric  waveguides  for  integrated  optics  were  the  first  to  demonstrate 
that  under  appropriate  conditions  some  modes  can  leak  instead  of  being 
purely  bound,  as  was  widely  expected.  Such  leakage  can  produce 
cross  talk  between  neighboring  circuits  and  deteriorate  performance. 
Furthermore,  these  leaky  modes  form  a  new  class  of  leaky  modes  since 
the  polarization  of  the  leakage  energy  is  opposite  to  that  of  the  excita¬ 
tion  energy,  in  contrast  to  all  other  known  cases. 

An  invited  paper23  has  recently  been  presented  which  discusses 
the  characteristics  of  these  new  leaky  modes,  explains  why  leakage 
occurs,  discusses  how  one  predicts  which  modes  are  leaky,  and  pre¬ 
sents  some  numerical  results.  Two  other  papers  are  in  the  process  of 
being  written,  one  stressing  additional  basic  material  and  another  on  a 
simple  graphical  method  for  predicting  when  modes  on  rib  waveguides 
will  leak. 

Previously,  our  numerical  values  for  leakage  were  computed  by 
including  only  the  surface  waves  and  their  interactions,  and  neglecting 
the  contributions  from  the  continuous  spectrum,  which  would  change  the 
numerical  values  somewhat  but  would  not  introduce  any  new  physical 
effects  or  alter  existing  ones.  Nevertheless,  now  that  accurate  results 
for  the  dielectric  step  junction  are  available  (see  the  previous  section), 
we  have  made  use  of  those  results  to  obtain  accurate  numerical  values 
which  include  all  modal  contributions.  We  therefore  obtained  during 
this  past  year  reliable  and  accurate  results  for  those  waveguides  for 
which  we  previously  had  only  approximate  values. 

A  paper24  by  Japanese  authors  appeared  in  March  1979  which 
observed  that  our  published  results  neglected  the  continuous  spectrum; 
they  then  presented  results  based  on  an  accurate,  but  still  approxi¬ 
mate,  method  for  taking  into  account  all  modal  contributions  and  solving 
for  the  leakage.  They  compared  their  results  with  our  admittedly  - 
approximate  ones,  and  concluded  that  ours  were  too  small  by  an  order 
of  magnitude.  In  their  paper,  they  criticized  the  usefulness  of  our 
approximate  procedure  because  of  the  large  discrepancy  However, 
they  later  informed  us  that  they  found  a  numerical  error  in  their  re¬ 
sults,  making  them  smaller  by  a  factor  of  (1.6n)2,  which  is  about  25. 
They  published  a  correction23  to  their  numerical  values,  and  presented 
a  new  curve  which  shows  much  closer  agreement  with  our  approximate 
data,  but  they  unfortunately  did  not  remove  their  earlier  critical  state¬ 
ment.  In  any  case,  when  we  compare  our  new  accurate  results  to  their 
published  ones,  we  find  very  good  agreement  except  in  the  vicinity  of 
the  resonance,  or  cancellation,  points.  Our  results  show  very  sharp 
nulls,  whereas  their  nulls  are  filled  in  substantially,  because  of  the 
approximate  method  they  used  to  compute  the  attenuation  constant  which 
yields  the  leakage. 

A  recent  paper26  by  two  Japanese  authors  describes  an  experiment 
they  performed  to  detect  leakage  from  a  rib-type  dielectric  structure. 
The  shape  of  the  measured  loss  curve  bears  striking  resemblance  to  the 
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calculated  leakage  curve;  the  numerical  values  differ  because  of  material 
loss  and  coupling  losses  at  the  feed  and  detection  ends.  The  authors 
conclude  that  their  measurements  verify  the  leakage  theory. 

D.  Anomalous  Absorption  Effects  in  Multilayered  and  Periodic 
Structures 


This  investigation  of  anomalous  absorption  effects  had  been  initi¬ 
ated2  by  the  observation21  that  total  absorption  can  occur  in  a  dielec¬ 
tric  grating  having  small  losses  even  though  all  of  the  incident  energy 
would  have  been  scattered  back  into  free  space  if  these  losses  were 
absent.  We  have  so  far  obtained  the  following  results. 

(1)  By  using  rigorous  methods,  we  have  verified6'12  preliminary 
indications2'22  that  the  anomalous  absorption  effect  is  due  to  a  wave 
interaction  phenomenon,  whereby  the  incident  field  couples  energy  to  a 
leaky  wave  that  can  be  supported  by  the  structure.  In  this  context, 
we  have  also  shown  that  the  phenomenon  can  then  also  occur  on  multi¬ 
layered  configurations  that  need  not  necessarily  include  any  periodicity. 

(2)  We  have  initiated  an  experimental  effort  to  verify  the  absorp¬ 
tion  phenomenon  in  multilayered  dielectric  media,  because  that  effect 
had  been  previously  observed  experimentally  only  in  periodic  structures 
involving  a  metal  substrate. 

(3)  We  have  provided  the  basis  for  a  qualitative  and  quantitative 
description14  of  the  absorption  phenomenon  in  terms  of  a  large  lateral 
beam  shift  of  the  Goos-Haenchen  type.  This  lateral  shift  is  due  to  the 
wave-interaction  process  that  occurs  when  the  finite  (Gaussian)  incident 
beam  excites  a  leaky-wave  field  which  forms  part  of  the  scattered 
energy . 

(4)  We  have  explored  the  possibility  of  using  the  anomalous 
absorption  effect  to  improve  the  response  of  photodiodes  by  maximizing 
its  quantum  efficiency  while  maintaining  minimal  response  time.  Our 
calculations  for  silicon  photodiodes  show  that  this  possibility  exists,  but 
the  conditions  for  its  realization  are  very  critical.  A  relaxation  of  such 
critical  conditions  could  be  obtained  by  imposing  a  periodicity  on  the 
diode  configuration,  but  this  necessitates  further  exploratory  work. 

(5)  We  have  concluded  our  experimental  verification15' 16  of  the 
anomalous  absorption  phenomenon  in  (non-periodic)  multilayered  con¬ 
figurations.  By  using  a  microwave  set-up  as  an  appropriate  model,  we 
have  obtained  nearly  total  (>  99%)  absorption  at  10.5  GHz  in  a  slightly 
lossy  (tan  6  =  0.07)  wood  layer,  which  showed  negligible  absorption  at 
frequencies  that  differ  by  more  than  0.5  GHz  away  from  that  central 
frequency.  This  has  been  the  first  observation  of  such  effects  in 
multilayered  media. 

(6)  We  have  developed  a  procedure17  for  exploring  the  total 
absorption  effect  in  terms  of  a  pole-zero  analysis  of  the  planar  scat¬ 
tering  coefficients  in  complex  wavenumber  space.  While  similar  analytic 
tools  have  been  used  in  the  past  in  network-theory  problems,  this 
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procedure  is  novel  in  the  context  of  wave  scattering  phenomena.  A 
preliminary  application  of  this  method17  has  already  provided  the  in¬ 
teresting  result  that  the  total  absorption  effect  can  be  obtained  in 
simple  canonic  structures  involving  not  more  than  a  single  layer  placed 
between  two  semi-infinite  spaces,  of  which  one  is  a  substrate.  More 
interestingly,  the  superstrate  may  then  either  be  a  denser  or  a  rarer 
medium.  The  former  possibility  was  the  only  one  that  had  been  pre¬ 
viously  known  for  practical  conditions,  but  the  latter  possibility  is  by 
far  more  useful  because  it  allows  incidence  from  air  rather  than  from  a 
denser  (additional)  medium,  such  as  a  prism. 

E.  Propagation  in  Inhomogeneous  Waveguides  by  Evanescent 
Waves  and  Complex  Ray  Techniques 

Difficulties  occur  when  one  attempts  to  generalize  the  evanescent 
wave  theory,  which  has  been  shown  to  be  highly  successful  from  non- 
leaky  modal  propagation  in  longitudinally  uniform,  lossless  graded  index 
waveguides,  to  curved  or  longitudinally  non-uniform  configurations 
where  leakage  occurs.  To  better  understand  the  leakage  problem  due 
to  curvature,  we  have  looked  at  the  canonical  problem  of  a  two- 
dimensional  cylindrically  curved  reactive  surface  layer  or  dielectric 
slab.  In  this  event,  exact  solutions  are  available  for  comparison. 

As  is  well  known,  a  mode  that  is  completely  trapped  in  a  straight 
waveguide  begins  to  radiate  on  the  convex  side  when  the  waveguide  is 
curved.  Loosely  speaking,  radiation  takes  place  when  the  local  phase 
velocity,  which  is  less  than  the  speed  of  light  in  the  exterior  medium 
near  the  surface,  approaches  the  speed  of  light  as  the  observer  moves 
away  from  the  surface.  Phrased  in  terms  of  evanescent  waves,  the 
angularly  propagating  local  evanescent  plane  wave  fields  become  less 
strongly  evanescent  with  distance  from  the  waveguide  surface.  In  the 
transition  region  to  the  radiative  regime,  their  character  changes  from 
evanescent  on  one  side  to  essentially  non-evanescent  (i.e.,  almost 
ordinary  local  plane  wave  or  ray  fields)  on  the  other  side. 

It  had  been  hoped  that  a  perturbation  approach  to  the  straight 
non-radiating  waveguide  mode  would  lead  to  a  correct  description  of  the 
(small)  leakage  when  the  waveguide  is  weakly  curved.  A  thorough 
study  revealed  that  this  is  not  the  case.  It  is  indeed  possible  to  find 
corrections  to  the  phase  paths  and  phase  fronts  on  either  side  of  the 
transition  region,  but  it  is  not  possible  to  connect  these  uniquely 
across  the  transition  region. 

One  way  to  avoid  this  difficulty  is  to  employ  complex  rays,  i.e., 
local  plane  wave  fields  with  complex  phase  which,  in  a  homogeneous 
medium,  proceed  along  straight  lines  in  a  complex  coordinate  space.  In 
curved  geometry,  complex  ray  fields  are  generated  by  complex  caustics 
(surfaces  of  revolution  in  complex  coordinate  space)  and,  depending  on 
the  point  of  origin  of  a  ray  on  the  caustic,  may  represent  fields  in  real 
space  that  are  strongly  or  weakly  evanescent.  Thus,  the  disjoint 
family  of  evanescent  waves  produced  by  direct  tracking  in  real  space 
becomes  connected  across  the  transition  region  when  described  in  terms 
of  complex  rays.  Complex  ray  theory  therefore  unifies  the  class  of 
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phenomena  associated  with  "evanescent  tunneling"  of  a  wave  field  to  the 
radiating  regime  because  of  curvature  of  the  initial  (caustic)  surface. 
Once  the  complex  ray  fields  have  been  determined,  they  may  be  inter¬ 
preted  in  real  space  in  terms  of  the  physically  meaningful  local  evanes¬ 
cent  waves,  and  one  may  now  chart  the  complete  and  unique  phase 
fronts  and  phase  paths  of  a  leaky  field.  Numerical  calculations  have 
been  made  to  exhibit  these  quantities  for  wave  fields  with  varying 
degrees  of  leakage. 

Having  understood  the  ray  fields  generated  by  a  complex  caustic, 
we  have  then  employed  these  fields  to  synthesize  the  guided  modes  that 
may  exist  on  circularly  curved  waveguides  with  prescribed  surface 
reactance  or  slab  parameters .  One  now  requires  both  outgoing  and 
incoming  ray  fields  with  respect  to  the  caustic  to  describe  modal  tun¬ 
neling.  Invoking  ray  field  closure  in  complex  space,  analogous  to  the 
"transverse  resonance"  relation  in  microwave  network  theory,  one 
determines  the  self-consistent  caustic  parameters  for  "source-free" 
propagation .  These  parameters  yield  the  eigenvalues  of  the  guided 
modes . 

We  may  therefore  conclude  that  complex  ray  theory  provides  the 
unifying  tool  for  asymptotic  analysis  of  modal  propagation  on  curved 
waveguides.  The  analysis  requires  extension  of  the  rays  as  well  as  the 
guiding  surfaces  into  a  complex  coordinate  space.  The  analytic  con¬ 
tinuation  is  unique,  and  complex  ray  tracing  numerically  tractable  (for 
complex  ray  tracing  applied  to  beam  fed  parabolic  reflector  antennas, 
see  reference  27).  When  leakage  is  strong  so  that  the  transition  region 
is  near  the  surface,  one  must  employ  uniform  transition  functions. 
These  can  be  determined  from  rigorous  canonical  prototypes. 

Having  understood  the  complex  ray  description  for  circularly  bent 
waveguides,  we  have  also  sought  to  generalize  these  results  to  accom¬ 
modate  weak  deviations  of  the  guide  axis  from  the  circular  case.  We 
have  found  in  preliminary  studies  that,  in  general,  the  permissible 
deviations  are  tied  to  the  geometry  and  do  not  appear  to  scale  with 
frequency;  i.e.,  increasing  the  frequency  does  not  permit  larger  devia¬ 
tions.  This  is  in  contrast  to  what  is  known  for  ordinary  geometric 
(real)  ray  fields,  where  surface  features  must  change  slowly  with 
respect  to  the  local  wavelength;  thus,  "strong"  geometric  changes  are 
permitted  for  sufficiently  small  wavelengths.  For  complex  rays,  "local¬ 
ity"  may  exist  in  complex  space  but  may  translate  into  "non-locality"  in 
real  space.  This  observation  has  also  been  made  in  connection  with 
evanescent  wave  and  complex  ray  tracing  from  large  apertures.28 

In  summary,  the  investigation?  described  above  have  unified  the 
asymptotic  treatment  of  modal  propagation  on  curved  waveguides,  and 
have  also  clarified  the  difficulties  encountered  when  one  attempts  to 
generalize  leakage  phenomena  from  canonical  (circular)  to  non-canonical 
configurations.  The  lesson  learned  is  that  such  generalizations,  usually 
done  by  perturbation,  must  be  approached  with  caution.  More  work 
needs  to  be  done  in  this  important  problem  area. 
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F.  Surface  Waves  Propagating  Obliquely  on  Periodic  Planar  Dielec- 
trie  Waveguides 

We  have  this  past  year  made  calculations  on  a  structure  not  pre¬ 
viously  included  in  our  three-year  program  because  of  the  sudden 
availability  of  very  accurate  experimental  data  and  our  development  of  a 
somewhat  similar  theoretical  procedure29  in  another  context  for  milli¬ 
meter  waves.  The  topic  is  that  of  optical  surface  waves  propagating  at 
an  angle  on  a  periodically-corrugated  dielectric  layer  on  a  dielectric 
substrate.  Such  structures,  with  obliquely-incident  waves,  have  been 
finding  increasing  application  lately  in  such  components  as  multiplexers, 
filters,  and  mode  deflectors. 

The  guidance  of  optical  surface  waves  propagating  in  a  direction 
normal  to  the  grooves  on  a  periodically-corrugated  planar  dielectric 
waveguide  is  well  known,  and  various  theoretical  treatments  are  avail¬ 
able.  Such  structures  have  been  used  in  grating  couplers  and  in 
distributed-feedback  lasers,  for  example.  The  electromagnetic  boundary 
value  problem  for  this  case  is  scalar,  and  TE  and  TM  surface  waves 
remain  independent  of  each  other.  When  the  surface  wave  propagates 
at  an  oblique  angle  with  respect  to  the  grooves,  the  TE  and  TM  surface 
waves  are  now  coupled  together,  and  the  boundary  value  problem 
becomes  a  three-dimensional  vector  one.  The  problem  is  one  for  which 
an  exact  solution  has  not  previously  been  available. 

A  corresponding  guidance  problem  was  developed  as  a  basic  inter¬ 
mediary  step  in  the  analysis  of  a  class  of  millimeter-wave  antennas  in 
Work  Unit  EMI-1.  It  was  originally  intended  there  to  employ  the  solu¬ 
tion  only  in  the  range  of  parameters  for  which  a  single  leaky  wave  is 
produced,  consistent  with  the  antenna  application.  In  the  optical  case, 
the  geometry  is  somewhat  different,  in  that  the  metal  ground  plane  must 
be  replaced  by  a  dielectric  substrate,  and  the  applications  of  interest 
are  usually  in  the  non-radiating  region,  especially  near  or  in  the  stop 
bands . 

There  was  originally  no  thought  to  apply  the  millimeter-wave 
antenna  solution  to  the  optical  domain,  but  the  redirection  was  stimu¬ 
lated  by  the  appearance  of  accurate  experimental  data  in  the  stop  band 
regions.  We  first  learned  about  these  experimental  results  at  a  Work¬ 
shop30  and  then  further  at  an  Optical  Society  Topical  Meeting.31  The 
speaker.  Dr.  Reinhard  Ulrich  of  Germany,  indicated  that  no  accurate 
theoretical  data  were  available  with  which  he  could  compare  his  mea¬ 
surements,  and  that  a  suitable  theoretical  expression  would  be  very 
valuable . 

We  therefore  revised  the  expressions  developed  for  the  millimeter 
wave  case  by  replacing  the  ground  plane  by  the  substrate,  and  then 
appropriately  modified  the  computer  program.  Detailed  calculations  were 
made  for  the  dispersion  curves  corresponding  to  the  specific  structure 
measured  by  Dr.  Ulrich,  and  the  agreement  between  the  measured  and 
theoretical  results  were  remarkably  good.  It  is  important  to  recognize 
that  these  calculations  were  made  in  and  near  the  stop  band  regions, 
which  are  the  regions  most  challenging  to  the  theory  since  the  space 
harmonics  must  appropriately  pair  off  properly  there. 
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These  dispersion  curves  in  the  stop- band  region  show  that,  due  to 
the  TE-TM  coupling,  two  extra  stop  bands  appear.  The  implications  of 
these  extra  stop  bands  were  recently  investigated  by  us  in  certain 
millimeter  wave  dielectric  devices  which  operate  in  the  Bragg  reflection 
region  to  produce  strong  reflections  with  negligible  radiation.  The 
TE-TM  coupling  turns  out  to  produce  an  extra  (unwanted)  transmission 
dip  or  reflection  peak,  and  thus  to  impair  the  device  performance.  The 
relevance  here  is  that  a  paper32  on  this  study  was  recently  presented 
which  also  includes  the  comparison  discussed  above  between  Dr. 
Ulrich's  measurements  and  our  theoretical  calculations. 

It  is  also  of  interest  that  further  examination  of  our  rigorous 
solution  for  the  oblique  guidance  of  surface  waves  reveals  that  a  rich 
variety  of  physical  effects  arise  as  a  result  of  oblique  guidance,  which 
are  not  present  at  all  for  normal  guidance.  Among  these  effects  are 
anisotropy  phenomena,  beam  steering,  radiation  of  beams  at  peculiar 
skew  angles,  and  strong  cross -polarization  effects  in  portions  of  the 
radiation  region.  We  hope  to  explore  these  effects  further;  however, 
some  of  these  effects  were  presented  by  us  in  a  recent  invited  talk.33 
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5.  DoD  AND  OTHER  INTERACTIONS 

(a)  Professor  Peng  spent  two  summers  (1977,  1978),  at  the  US 
Army  CORADCOM  at  Fort  Monmouth,  New  Jersey,  working  with  Dr. 
Felix  Schwering  on  dielectric  grating  and  taper  antennas,  under  an 
LRCP  arrangement  with  the  US  Army  Research  Office.  These  interac¬ 
tions  have  continued  in  the  form  of  a  Post-LRCP  contract  and  successor 
contracts. 

(b)  Professor  Peng  is  currently  (1980)  performing  a  task  on 
waveguides  for  integrated  optics  for  Dr.  John  Zavada  of  the  US  Army 
ARRADCOM,  Dover,  New  Jersey,  under  arrangements  made  through  the 
Battelle  Memorial  Institute. 

(c)  We  are  discussing  arrangements  for  a  collaborative  program 
with  Professor  R.  Ulrich  of  the  University  of  Hamburg-Harburg  on  the 
topic  of  new  physical  effects  due  to  oblique  guidance  on  optical  grating 
structures,  where  we  would  do  the  theory  and  his  group  would  perform 
careful  measurements. 

(d)  Collaboration  with  Dr.  G.  Jacobsen,  Technical  Univ.  of  Den¬ 
mark,  Lyngby,  Denmark,  is  continuing  on  hybrid  methods  applied  to 
beam  propagation  in  optical  fibers. 

(e)  Collaboration  with  Professor  A.M.  Scheggi,  Institute  for 
Electromagnetic  Wave  Propagation,  Florence,  Italy,  is  continuing  on 
hybrid  methods  applied  to  optical  fibers. 

(f)  The  studies  on  optical  beam  couplers  and  infrared  gratings 
have  led  to  collaboration  and  interaction  with  Dr.  P.K.  Cheo's  group  at 
United  Technologies,  Hartford,  Connecticut,  to  whom  we  have  supplied 
analytical  solutions  and  numerical  data. 

(g)  The  use  of  blazing  for  optimizing  the  performance  of  periodic 
beam  couplers,  which  had  originally  been  proposed  by  us,  has  been 
explored  by  Dr.  Hammer's  group  at  RCA  Princeton  Laboratories  and  has 
stimulated  other  efforts  for  designing  and  fabricating  such  blazed 
gratings . 

(h)  Comparisons  of  analytical  procedures  and  numerical  data  on 
beam  and  periodic-structure  problems  have  occurred  on  a  continuing 
basis  with  Dr.  Streifer  of  Xerox  Corp.,  Palo  Alto,  California. 
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E.  EVANESCENT  WAVE  TRACKING  -  A  NEW  APPROACH  TO  THE 
ANALYSIS  OF  LARGE  REFLECTOR  AND  LARGE  APERTURE  AN¬ 
TENNAS 


Professor  L.B.  Felsen 

Unit  EMI -5 

1.  OBJECTIVE(S) 

To  develop  a  new  method  for  analyzing  the  radiation  properties  of 
large  aperture  and  large  reflector  antennas,  based  on  the  local  tracking 
of  evanescent  wave  fields  from  the  feed  to  subreflector(s)  and  main 
reflector  and  then  to  the  far  field,  without  intervening  integrations 
over  an  equivalent  aperture  plane. 

The  radiation  characteristics  of  large  reflector  and  aperture  an¬ 
tennas  are  analyzed  conventionally  by  the  method  of  physical  optics  or 
the  Geometrical  Theory  of  Diffraction  (GTD).  By  these  procedures, 
one  may  establish  equivalent  fields  on  the  reflector  surfaces  or  in  an 
aperture  plane,  and  then  perform  Kirchhoff  type  integrations  to  obtain 
the  far  field.  These  integrations,  when  done  numerically,  are  costly 
and  time-consuming.  While  GTD  (with  uniform  corrections)  may  be 
employed  to  yield  the  near  and  far  field  directly  for  certain  special 
configurations,  this  is  not  possible  for  strongly  tapered  illuminations. 

To  avoid  the  need  for  integration  over  one  or  several  equivalent 
aperture  planes,  a  new  method,  based  on  local  evanescent  wave  track¬ 
ing  will  be  employed. 

2.  APPROACH 

The  determination  of  the  phase  paths  poses  a  principal  difficulty  in 
applying  the  evanescent  wave  method,  especially  for  large  distances 
from  the  initial  surface.  However,  one  may  show  that  the  phase  paths 
at  all  ranges  are  approximately  hyperbolic  in  the  weakly  evanescent 
region  near  the  field  maximum.  For  symmetrical  field  distributions  with 
Gaussian  amplitude  taper,  all  phase  paths  in  the  paraxial  region  belong 
to  a  confocal  elliptic  coordinate  system,  and  this  statement  is  equivalent 
to  saying  that  the  paraxial  evanescent  field  in  real  space  is  generated 
by  a  line  source  (2-dimensional  field)  or  point  source  (3-dimensional 
field)  displaced  into  a  complex  coordinate  space.  It  is  proposed  to 
approximate  amplitude  tapers  that  depart  from  the  symmetrical  Gaussian 
shape  locally  by  a  Gaussian  profile,  thereby  providing  a  global  evanes¬ 
cent  field  generated  by  a  complex  caustic.  The  corresponding  phase 
paths  in  real  space  would  then  still  be  hyperbolic  but  not  confocal. 
This  feature  should  systematize  and  simplify  the  tracking  problem. 

3.  PROGRESS 

Referring  to  the  previous  proposal,  it  had  been  hoped  that  hyper¬ 
bolic  paths,  which  correspond  to  a  particular  "canonical"  tapered  aper¬ 
ture  distribution,  could  be  matched  locally  to  aperture  fields  that  de- 
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part  from  the  globally  canonical  hyperbolic  form,  at  least  in  the 
(paraxial)  weakly  evanescent  region  near  the  maximum  of  beam- type 
fields.  However,  by  comparison  with  independently  calculated  fields  for 
various  non-canonical  aperture  profiles,  we  have  found  that  this  is  not 
generally  applicable,  although  the  procedure  is  adequate  and  advan¬ 
tageous  for  certain  types  of  initial  distributions.  The  deviation  of  the 
true  phase  paths  from  a  hyperbolic  shape  can  generally  not  be  ignored 
for  accurate  calculation  of  the  field. 

The  difficulties  may  be  overcome  by  calculating  the  phase  paths 
and  fields  via  complex  ray  theory,  which  also  describes  the  propagation 
properties  of  local  plane  wave  fields  with  complex  phase.  Complex  rays 
follow  straight  trajectories  in  a  complex  coordinate  space.  They  ori¬ 
ginate  on  a  complex  initial  surface  with  complex  initial  conditions 
generated  by  analytic  continuation  of  the  specified  real  initial  data.  An 
observable  field  is  associated  with  the  real-space  intersection  of  a 
complex  ray;  the  true  phase  paths  of  evanescent  wave  theory  (EWT) 
correspond  to  those  real-space  intersections  of  complex  ray  fields  satis¬ 
fying  the  constraint  of  constant  exponential  amplitude.  While  ray 
tracing  in  a  complex  coordinate  space  is  generally  quite  involved,  the 
procedure  does  become  tractable  for  "slightly  complex"  rays  character¬ 
izing  the  weakly  evanescent  fields  near  the  beam  maximum.  However, 
the  complex  ray  fields  per  se,  unlike  the  real-space  fields  of  EWT,  do 
not  have  a  simple  physical  interpretation.  It  has  therefore  been  con¬ 
cluded  that  complex  ray  theory  should  be  regarded  essentially  as  an 
algorithm  for  solving  the  EWT  trajectory  and  field  equations.  This 
combined  complex-ray  and  EWT  approach  is  tractable,  both  analytically 
and  numerically,  and  furthermore  describes  the  complete  evolution  of 
the  field  from  the  near  zone  to  the  far  zone  in  terms  of  physically 
significant  wave  processes. 

The  basic  EWT  and  complex  ray  studies  outlined  above  have  been 
presented  in  a  manuscript  that  has  been  accepted  for  publication . 1  To 
further  test  the  validity  of  the  conclusions  from  this  investigation,  the 
EWT -complex  ray  approach  has  been  applied  to  the  direct  tracking  of  an 
initial  beam-type  field  from  the  source  plane  to  the  far  zone  via  an 
intervening  reflection  at  a  paraboloidal  surface.  Here,  the  complex  ray 
tracing  required  extension  of  the  initial  conditions  as  well  as  the  re¬ 
flector  geometry  into  a  complex  coordinate  space.  The  success  of  the 
procedure  for  the  single  reflector  problem  is  documented  in  another 
paper  accepted  for  publication.2 

4.  PUBLICATIONS 

1.  P.D.  Einziger  and  L.B.  Felsen,  "Evanescent  Waves  and  Com¬ 
plex  Rays,"  accepted  for  publication  in  IEEE  Trans.  Ant.  and 
Propag . 

2.  F.J.V.  Hasselman  and  L.B.  Felsen,  "Asymptotic  Analysis  of 
Large  Parabolic  Reflector  Antennas,"  accepted  for  publication 
in  IEEE  Trans.  Ant.  and  Propag. 
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5.  DoD  AND  OTHER  INTERACTIONS 

Collaboration  with  Professor  R.  Zich,  Technical  Univ.,  Turin, 
Italy,  and  Professor  F.  Hasselman,  Catholic  Univ.,  Rio  de  Janeiro, 
Brazil  on  evanescent  wave  and  complex  ray  methods  for  large  reflector 
antennas . 


SECTION  II:  SOLID  STATE  ELECTRONICS 


A.  X-RAY  GUIDED  WAVE  ELECTRONICS 


Professors  B.  Post  and  H.J.  Juretschke 

Unit  SSI-1 

1.  OBJECTIVE(S) 

To  use  the  anomalous  propagation  of  x-rays  in  good  crystals  for 
developing  guided-wave  ipodes  with  optimal  properties  in  thick  crystals; 
to  apply  these  modes  to!  both  the  study  of  the  solid  state  of  the  host 
crystal  and  the  coupling  of  x-rays  to  other  waves  within  the  solid;  and 
to  examine  coupled-mode  techniques  to  determine  their  suitability  for 
refining  structure  determinations  in  crystallography. 

Dynamical  interactions  of  multiple  x-rays  within  nearly  perfect 
crystals  are  intended  to  lead  eventually  to  special  waveguides  and 
components  that  can  form  the  basis  for  possible  future  coherent  x-ray 
circuitry.  A  second  portion  of  the  program  involves  the  coupling  of 
x-rays  to  other  wave  types,  such  as  optical  or  acoustic  waves,  within 
the  crystal.  If  the  interactions  prove  to  be  sufficiently  strong,  these 
other  waves  can  be  used  to  modulate  the  x-ray  beam,  or  to  switch  it 
from  one  state  to  another.  Thus,  the  x-ray  beam  could  be  controlled 
electronically ;  otherwise ,  it  would  be  necessary  to  turn  the  crystal 
mechanically  to  achieve  the  desired  changes.  The  proposed  program  on 
x-rays  thus  involves  a  number  of  fundamental  investigations  which 
could  make  possible  a  systematic  approach  to  x-ray  circuitry. 

2.  APPROACH 

Based  on  the  dynamical  theory  of  n-beam  interactions  of  x-rays  in 
crystals,  and  on  the  availability  of  good  single  crystals  of  ever  increas¬ 
ing  variety,  we  will  carry  out  a  combined  program  of  theoretical  ex¬ 
ploration  and  experimental  realization  of  the  behavior  of  x-ray  modes 
traveling  along  directions  close  to  those  specifying  two  or  more  Bragg 
conditions  within  the  crystal.  Since  strong  dynamical  effects  occur  only 
within  minutes  of  arc  or  less  of  these  directions,  the  experimental 
methods  require  high  spatial  and  energy  resolution.  Coupling  to  other 
types  of  waves  will  employ  the  general  techniques  of  nonlinear  optics 
already  well-developed  in  other  frequency  domains. 

3.  PROGRESS 

A.  Improvements  in  Instrumentation 

(1)  The  new  5kVA  x-ray  generator  has  been  set  up  for  very  high 
intensity  output,  using  special  filaments,  and  it  has  been  fitted  with  a 
grooved  crystal  Ge  monochromator  that  has  reduced  the  incident  beam 
divergence  to  about  10  seconds  of  arc. 

(2)  For  Renninger  type  experiments,  where  the  divergence  must 
be  small  in  any  direction,  we  have  set  up  pin  holes  about  1  m  from  the 
0.1  by  0.1  mm  target  to  obtain  less  than  30"  divergence. 
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(3)  The  use  of  divergent  beams  transmitting  x-rays  through  the 
crystal  in  order  to  locate  readily  regions  of  many  beam  interactions, 
using  the  same  type  of  target  as  in  (2),  has  reduced  the  observation 
time  on  the  rotating  anode  generator  by  a  factor  of  10  relative  to  the 
previously  used  stationary  microbeam  target. 

(4)  We  have  acquired  a  Cr  rotating  target  to  extend  the  high 
intensity  range  of  the  instrument  to  other  wavelengths. 

B.  N-Beam  Interactions 

(1)  Efforts  to  extend  phase  determination  from  simple  relatively 
perfect  crystals1'2  to  mosaic  organic  crystals  have  not  been  successful, 
when  Cu  K  radiation  was  used.  Phase  effects  are  there,  but  smear 
out.  The  use  of  Cr  seemed  to  help.  It  increased  the  absorption  by 
the  specimen  and,  as  discussed  in  the  last  report,  such  increases 
facilitate  the  observation  of  dynamic  effects  on  the  phase.  Neverthe¬ 
less,  it  appears  at  present  unlikely  that  the  use  of  Cr  will  make  pos¬ 
sible  routine  exploration  of  the  phase  determination.  Still  longer  wave¬ 
lengths  would  be  needed,  and  for  that  we  will  use  synchrotron  radia¬ 
tion,  probably  at  Brookhaven. 

(2)  Studies  of  polarization  effects  in  3-beam  diffraction  have  been 
completed.3  This  will  be  included  in  a  comprehensive  review  of  3 -beam 
diffraction  to  be  published  in  Z.  fuer  Naturforschung,  in  early  1982. 

(3)  We  have  collaborated  with  a  group  at  the  University  of  Am¬ 
sterdam  in  a  study  of  the  relative  merits  of  calculating  phase  proba¬ 
bilities  based  explicitly  on  dynamical  diffraction  theory  compared  with 
statistical  methods  now  widely  used.  This  study  indicates  that,  even  in 
its  present  primitive  form,  the  dynamic  approach  is  as  successful  in 
predicting  probable  phases  as  are  these  other  methods.  Results  of  this 
work  are  being  prepared  for  publication  in  Acta  Cryst. 

(4)  A  modification  of  our  experimental  approach  to  phase  deter¬ 
mination  has  been  tried  recently .  *  It  shows  great  promise  for  the 
determination  of  phases.  In  essence,  it  involves  the  application  of  the 
methods  developed  at  the  Polytechnic  to  the  analysis  of  the  intensity 
distribution  of  simultaneous  reflections,  but  now  recorded  in  the  Bragg 
mode,  under  conditions  of  high  resolution.  The  technique  is  simple  and 
rapid.  It  has  been  successfully  applied  to  the  analysis  of  patterns  of 
several  simple  substances.  Work  is  proceeding  to  optimize  the  proce¬ 
dure  and  to  apply  the  method  to  the  determination  of  phases  in  mosaic 
crystals.  The  restriction  of  phase  determination  to  cases  where  all 
structure  factors  involved  are  essentially  of  equal  magnitude  does  not 
apply  to  the  newer  procedure.  Any  range  of  F's  is  usable. 

C.  Interaction  of  X-Rays  with  Other  Waves 

The  progress  in  understanding  the  dynamic  interaction  of  x-rays 
with  phonons,  or  with  the  phonon  portion  of  polaritons  in  polar  crys¬ 
tals,  has  been  very  satisfactory.  The  major  new  advance  has  been  in 
treating  this  problem  as  an  x-ray  fully  interacting  with  many  phonons. 
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This  approach  is  obviously  needed  when  dealing  with  a  thermal  phonon 
distribution,  but  it  was  also  found  to  be  essential  for  the  interpretation 
of  an  experiment  of  x-rays  coupled  to  intense  "monochromatic"  phon¬ 
ons.5  In  both  cases,  the  n-phonon  problem  still  reduces  in  practice  to 
a  set  of  n  one-phonon  problems,  as  far  as  the  satellites  are  concerned. 
However,  the  line  shape  and  the  amplitude  of  the  main  reflection  are 
influenced  simultaneously  by  all  phonons,  and  can  be  drastically  al¬ 
tered.6  It  is  this  discrepancy  between  observed  and  calculated  inten¬ 
sities  of  the  main  reflection,  already  alluded  to  in  the  last  progress 
report,  that  has  delayed  the  publication  of  the  two  papers  that  have 
been  in  preparation. 

The  specific  areas  of  progress  are: 

(1)  We  have  found  a  fully  analytic  solution  of  the  six-beam  case 
of  an  x-ray  interacting  dynamically  with  monochromatic  phonons  of 
intensity  which  may  be  small  or  large  moving  parallel  to  the  reflecting 
planes  of  the  crystal.  In  this  formulation,  the  problem  factors  into 
cubic  and  quadratic  equations ,  and  gives  explicit  solutions  for  the 
three-sheet  dispersion  surface  and  the  sixtuplets  of  x-ray  field  ampli¬ 
tudes  associated  with  each  point  on  this  surface.  We  can  therefore 
solve  the  boundary  problem  exactly,  and  determine  explicitly  all  ex¬ 
ternal  reflected  fields,  for  arbitrary  intensity  of  the  phonon  beam. 
This  exact  solution  has  been  extremely  useful  in  exploring  the  conse¬ 
quence  of  changing  the  various  parameters  in  the  problem,  and  es¬ 
pecially  under  what  condition  the  usual  weak  coupling  approximation 
holds.  We  have  studied  solutions  of  both  model  problems  and  of  the 
specific  interaction  of  Mo  K-radiation  reflected  from  the  (004)  planes  of 
InSb.7 

(2)  We  have  formulated  the  dynamic  n-beam  phonon  problem,  and 
set  up  the  explicit  equations  for  the  2n+l  sheeted  dispersion  surface, 
and  all  the  associated  fields.  The  procedure  follows  closely  the  method 
of  item  1.  Introduction  of  the  boundary  conditions  leads  to  a  2n  by  2n 
determinant  for  the  relative  incident  fields,  and  also  to  the  correspond¬ 
ing  expressions  for  the  reflected  fields.  When  the  coupling  is  weak, 
these  results  can  be  cast  into  explicit  modifications  of  the  uncoupled 
solutions  for  all  but  the  main  reflection.  This  reflection  has  to  be 
solved  using  a  modified  dispersion  relation.  Our  expressions  show  so 
far  that  this  modification  contains  parts  equivalent  to  the  Debye-Waller 
factor,  but  also  other  contributions  that  appear  only  when  the  full 
dynamical  problem  is  considered  before  making  the  weak  coupling  ap¬ 
proximation. 

(3)  We  have  evaluated  the  modification  of  the  main  reflected 
intensity  for  many  thermal  phonons,  when  the  reflection  occurs  away 
from  the  Bragg  angle,  and  find  that  the  multi-phonon  coupling  leads  to 
a  significant  reduction  from  that  expected  for  the  uncoupled  two-beam 
case,  or  the  case  with  the  traditional  correction  for  thermal  diffuse 
scattering.8  Numerical  calculation  is  now  in  progress  to  determine  this 
reduction  quantitatively. 
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(4)  Reanalysis  of  the  experimental  results  of  the  x-ray  intense 
phonon  beam  interaction  in  InSb  in  light  of  the  exact  solution  of  item  1, 
above,  has  shown  that  the  phonon  beam  used  was,  in  fact,  far  from 
monochromatic,  but  had  a  spread  in  energy  larger  than  that  of  the 
original  Bragg  peak,  although  its  direction  was  still  well  defined.  We 
have  obtained  specific  expressions  for  the  influence  of  this  phonon 
distribution  on  the  main  line ,  where  experimentally  the  intensity  was 
seen  to  be  reduced  by  5  to  20%. 

(5)  The  interaction  of  x-rays  with  laser-induced  polaritons  in 
polar  crystals  has  been  analyzed  from  the  point  of  view  of  energy 
transfer,  and  it  is  found  that  under  optimal  conditions,  using  normal 
laser  intensities,  the  weak  coupling  approximation  for  the  phonon  part 
of  the  polariton  is  justified.  This  has  made  it  possible  to  adapt  directly 
the  formalisms  worked  out  under  1,  above,  to  single  phonons  traveling 
in  arbitrary  directions,  but  with  the  additional  condition  that  energy 
transfer  between  polariton  and  x-ray  is  not  negligible.  This  implies 
that  the  coupled  modes  contain  fields  of  differing  frequencies,  that 
refer  to  different  dispersion  surfaces.  However,  our  formalism  lends 
itself  readily  to  incorporate  this  complication .  One  immediate  conse¬ 
quence  is  that  the  expected  six-beam  coupling  reduces  to  four-beam 
cases,  i.e.,  polaritons  are  predominantly  either  absorbed  or  emitted, 
but  not  both.  We  are  in  the  process  of  determining  the  optimum  ex¬ 
perimental  configuration  for  observing  the  predicted  laser-induced 
switching  of  x-rays. 

The  work  reported  here  under  items  (1)  to  (5)  represents  a  Ph.D. 
thesis  in  progress  that  is  expected  to  be  completed  during  the  fall 
semester  of  1981. 6  Various  parts  of  this  thesis  will  be  prepared  for 
publication . 

4.  PUBLICATIONS 

1.  B.  Post,  Acta  Cryst.  A  35,  17  (1979). 

2.  P.  Wang,  "Direct  Experimental  Detection  of  X-Ray  Phases 
from  Intensity  Measurements,"  Ph.D.  Thesis,  Polytechnic 
Institute  of  New  York  (1979). 

3.  T.  Horn,  "Three  Beam  Dynamical  Interactions,"  Ph.D.  Thesis, 
Polytechnic  Institute  of  New  York  (1979). 

4.  P.  Gong,  "Phase  Determination  from  X-Ray  Intensities," 
Ph.D.  Thesis,  Polytechnic  Institute  of  New  York,  in  pro¬ 
gress  . 

5.  S.D.  LeRoux,  R.  Colella  and  R.  Bray,  Phys.  Rev.  Lett.  35, 

230  (1975).  ~ 

6.  F.  Robbins,  "X-Ray-Phonon  Interaction  as  a  Dynamical  N- 
Beam  Problem,"  Ph.D.  Thesis,  Polytechnic  Institute  of  New 
York,  in  progress. 
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7.  F.  Robbins  and  H.J.  Juretschke,  "Dynamic  Effects  in  X-Ray 
Monochromatic  Intense  Phonon  Interaction  in  Symmetric  Bragg 
Reflection,"  to  be  published. 

8.  F.  Robbins  and  H.J.  Juretschke,  "Dynamic  Effects  in  X-Ray 
Thermal  Phonon  Interaction  in  Symmetric  Bragg  Reflection," 
to  be  published. 

5.  DoD  AND  OTHER  INTERACTIONS 

(a)  The  American  Crystallographic  Association  will  hold  a  special 
session  honoring  Professor  Post’s  contributions  to  crystallography  at  the 
Spring  1982  meeting  at  Gaithersburg,  MD.  Both  Professor  Juretschke 
and  Professor  Post  have  been  invited  to  present  papers. 

(b)  Some  of  the  above  proposed  work  may  involve  specific  collabo¬ 
ration  with  laboratories  in  Brazil  (Professor  Chang),  Greece  (Professor 
Alexandropoulos),  and  Australia  (Professor  Wagenfeld).  The  collabora¬ 
tion  is  being  proposed  by  them  in  order  to  participate  in  our  research. 

(c)  Professor  Post  has  been  presenting  invited  talks  on  aspects  of 
the  problem  of  phase  determination  in  Amsterdam,  Oslo  and  Trondheim, 
during  the  summer  of  1980.  He  also  presented  a  paper  on  the  statis¬ 
tical  aspects  of  phase  determination,  discussed  above,  at  the  meeting  of 
the  American  Crystallographic  Association,  Ottawa,  in  August  1981. 


« 
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B.  INTERACTION  OF  MILLIMETER  WAVES  AND  SEMICONDUCTORS 


Professor  B.  Senitzky 

Unit  SSI-2 

1.  OBJECTIVE(S) 

The  object  of  this  study  is  to  investigate  the  interaction  of  milli¬ 
meter  wavelengths  with  mobile  charge  carriers  in  semiconductors. 

The  effect  of  conduction  carriers  on  guided  energy  is  character¬ 
ized  by  a  complex  dielectric  constant  which,  in  appropriate  structures, 
attenuates  and/or  changes  the  phase  of  the  guided  radiation.  This 
principle  is  the  basis  for  the  familiar  class  of  devices  which  includes 
modulators,  phase  shifters,  and  switches.  It  is  therefore  important  to 
have  a  quantitative  understanding  of  the  interaction  of  millimeter  wave 
radiation  and  the  non -equilibrium  charge  carriers  which  are  injected 
from  the  terminals.  In  addition,  it  is  also  important  to  understand  the 
interaction  of  the  equilibrium  charge  carriers  with  the  millimeter  wave 
radiation  since  this  will  determine  the  electromagnetic  properties  of  the 
active  structure  in  the  absence  of  injected  carriers. 

At  present,  the  commonly  accepted  approach  is  to  estimate  the 
dielectric  constant  by  assuming  a  simple  Drude  model  with  a  single 
collisional  relaxation  time  which  is  determined  from  DC  resistivity  vs. 
concentration  characteristics.  The  validity  of  this  assumption  has  been 
questioned  both  theoretically  and  experimentally,  at  microwaves  and 
infra-red,  but  not  at  millimeter  wavelengths.  It  is  our  primary  goal  to 
investigate,  both  theoretically  and  experimentally,  the  interaction  of 
radiation  with  free  carriers  in  semiconductors  at  millimeter  wavelengths. 

2.  APPROACH 

Our  original  approach  was  to  study  the  interaction  of  millimeter 
wave  energy  and  semiconductors  in  a  slotted  metallic  waveguide.  Be¬ 
cause  of  the  very  small  sizes  involved,  and  the  influence  of  discon¬ 
tinuity  structures  that  extended  the  fields  into  regions  outside  the 
waveguide,  this  approach  failed  to  yield  reliable  results.  We  are  now 
studying  the  interaction  by  measuring  the  reflection  and  transmission 
properties  of  silicon  wafers  in  free  space  with  transmitting  and  receiv¬ 
ing  horns. 

3.  PROGRESS 

Let  us  first  recall  that  semiconductor  material  is  used  to  guide  and 
control  millimeter  wave  radiation.  The  millimeter  wave  energy  is  guided 
in  devices  such  as  waveguides  and  antennas ,  and  controlled  by 
switches,  electronic  phase  shifters  and  modulators. 

The  control  of  the  millimeter  wave  energy  by  semiconductors  is 
usually  achieved  by  changing  the  number  of  conduction  electrons  or 
holes  in  the  electromagnetic  field.  In  most  applications  the  guiding 
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structure  for  the  electromagnetic  energy  is  designed  so  that  the  radia¬ 
tion  can  interact  effectively  with  the  semiconductor.  At  microwave 
frequencies,  matching  sections  have  been  developed  which  can  achieve 
the  required  interaction.  At  millimeter  wave  frequencies  the  matching 
structures  become  lossy,  and  factors  such  as  junction  capacity,  lead 
impedance,  and  encapsulation  further  degrade  device  performance.  New 
types  of  matching  structures  and  semiconductor  devices  must  therefore 
be  developed  for  this  frequency  range.  To  accomplish  this,  a  knowl¬ 
edge  of  the  complex  dielectric  constant  of  the  semiconductor  as  a  func¬ 
tion  of  the  mobile  carrier  density  is  required. 

Before  proceeding,  it  is  important  to  note  that  a  unique  dielectric 
constant  can  only  be  defined  under  certain  conditions.  Mobile  charge 
carriers  can  have  an  equilibrium  or  non-equilibrium  energy  distribution. 
For  example,  carriers  which  arise  from  a  background  impurity  level  and 
are  subject  to  a  weak  electromagnetic  field  are  in  an  equilibrium  distri¬ 
bution  and  can  be  characterized  by  a  unique  dielectric  constant.  On 
the  other  hand  carriers  which  are  injected  by  an  electric  field  from 
external  contacts  are,  strictly  speaking,  not  in  thermal  equilibrium  and 
may  or  may  not  be  characterized  by  the  same  dielectric  constant.  We 
will  be  studying  the  properties  of  carriers  in  thermal  equilibrium. 

The  dielectric  constant  of  semiconductors  at  microwave  frequencies 
is  usually  studied  in  guided  structures  whereas  at  optical  frequencies 
these  studies  are  conducted  in  free  space.  In  the  millimeter  wave 
region  one  could,  in  principle,  take  either  approach.  We  have  been 
able  to  obtain  more  reliable  results  by  free  space  measurements  which 
we  will  describe  below.  Our  measurements  were  performed  on  n-type 
single  crystal  silicon  at  a  frequency  of  107  GHz. 

We  measured  the  reflectance,  R,  and  transmittance,  T,  of  a  silicon 
wafer  of  uniform  thickness  for  plane  wave  incidence.  The  solution  of  R 
and  T  as  a  function  of  the  relative  permittivity  of  the  medium, 
k  =  k'  -jK~,  the  thickness,  t,  of  the  wafer  and  the  polarization  of  the 
incident  radiation  is  well  known .  We  determined  the  relative  permit¬ 
tivity  from  a  measurement  of  R,  T  and  t  using  the  system  shown  in 
Figure  1.  The  millimeter  wave  source,  manufactured  by  Varian,  emits 
radiation  at  107.3  GHz.  To  determine  the  reflectance,  R,  we  measured 
the  signal  at  receiving  horn  2  with  the  silicon  wafer  in  place  and  with 
the  silicon  wafer  replaced  by  an  aluminum  slab.  Assuming  the  aluminum 
slab  behaved  as  a  perfect  mirror,  the  reflectance  of  the  silicon  wafer 
could  be  determined  as  the  ratio  of  the  two  measurements.  The  trans¬ 
mittance  was  determined  as  the  ratio  of  the  signal  received  at  horn  1 
with  the  silicon  wafer  in  place  and  the  silicon  wafer  removed. 

When  the  silicon  wafer  is  positioned  so  that  its  surface  is  perpen¬ 
dicular  to  the  direction  of  propagation  of  the  incident  radiation,  reso¬ 
nances  between  the  horns  and  the  wafer  can  occur.  These  interactions 
have  the  form  of  a  residual  standing  wave  pattern  such  that  the  re¬ 
ceived  millimeter  wave  signal  is  a  sinusoidal  function  of  position  with  a 
period  of  X/2  (where  \  is  the  free  space  wavelength).  The  standing 
waves  were  practically  eliminated  by  rotating  the  silicon  wafer  so  that 
the  surface  normal  is  at  an  angle,  0,  with  the  direction  of  incidence. 
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In  our  experiment  we  used  an  angle  0  of  about  17  degrees.1  Although 
residual  standing  wave  interactions  can  still  occur,  they  are  minor,  and 
the  measurement  is  performed  by  averaging  the  received  signals  with 
the  receiving  horn  at  two  longitudinal  positions  spaced  \/4  apart.2'3 

Another  experimental  problem  that  we  encountered  was  a  deviation 
from  plane  wave  conditions  due  to  radiation  from  the  edges  of  the  wafer 
and  wafer  holder.  This  deviation  gave  rise  to  a  transverse  interference 
pattern.  This  effect  was  greatly  reduced  by  using  a  tapered  millimeter 
wave  absorber  to  cover  the  edges  of  the  wafer  and  wafer  holder. 

Our  measurements  were  performed  over  two  resistivity  ranges. 
Both  transmittance  and  reflectance  measurements  were  obtained  in  the 
5-46  ohm-cm  range  whereas  in  the  0.007-5  ohm-cm  range  only  reflec¬ 
tance  measurements  could  be  obtained  because  the  wafers  became  too 
absorbant  to  obtain  reliable  transmission  data.  In  Figure  2  the  points 
represent  the  values  of  k  '  and  k~  determined  from  our  experimental 
measurements  of  reflectance  and  transmittance.  In  Figure  3  the  points 
represent  our  measured  reflectance. 


RESISTIVITY.  C.  in  OHM-CM 

Fig.  2  The  real  and  imaginary  parts,  k'  and  k",  of  the  relative 
permittivity  as  a  function  of  resistivity  p.  The  points  are 
measured  and  the  solid  line  is  the  theoretical  prediction  based 
on  the  Drude  model. 
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Fig.  3  The  reflectance  R  as  a  function  of  resistivity  p.  The  points 
are  measured,  and  the  so’ d  line  is  the  theoretical  prediction 
based  on  the  Drude  model. 

Our  experimental  data  in  these  two  graphs  are  compared  with  the 
theoretical  predictions  of  a  simple  Drude  model  with  an  energy- 
independent  collision  frequency,  o.  For  an  applied  electric  field, 
varying  harmonically  as  exp(juit),  we  treat  the  electron  as  a  classical 
point  charge  subject  to  collisional  damping  and  write  the  relative  per¬ 
mittivity,  k,  of  silicon  as  a  sum  of  a  wavelength  independent  contri¬ 
bution  of  the  intrinsic  crystal,  and  a  varying  contribution  of  free 
carriers  arising  from  impurities, 

K  =  Kp  -  (U.2/U.2)  (1  -  j  ^y1  (1) 

where  uj2  =  Ne2/m*e0,  N  is  the  donor  impurity  concentration,  e  is  the 

electronic  charge,  m*  is  the  electron  effective  mass  and  e0  is  the 
vacuum  permittivity.  The  above  model  also  yields  the  expression  for 
resistivity 


P-1  =  <*»2  £o  u"1  (2) 

To  determine  k  we  first  measured  the  resistivity,  p,  of  the  silicon  wafer 
using  a  colinear  four  point  probe.  The  impurity  concentration,  N,  is 
determined  from  National  Bureau  of  Standards  data4  (N  is  a  function  of 
p).  The  collision  frequency,  u,  can  then  be  determined  from  (2)  above 
and  the  results  substituted  in  (1)  to  determine  the  relative  permittivity, 
k,  as  a  function  of  resistivity  as  shown  in  Figure  2.  From  these  values 
we  can  also  compute  the  reflectance,  R,  as  a  function  of  resistivity  as 
shown  in  Figure  3. 

The  high  resistivity  data  shown  in  Figure  2  indicates  that  the 
experimentally  determined  values  of  k'  and  are  in  good  agreement 
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with  the  predictions  of  the  Drude  model.  Our  lower  resistivity  data  in 
Figure  3  indicates  that  our  measured  reflectivity  is  also  in  good  agree¬ 
ment  with  the  Drude  model.  This  result  was  also  found  in  the  infrared 
region5  over  a  frequency  range  of  200-1000  cm"1  (1  cm  1  =  30  GHz). 

On  the  basis  of  the  above  data  we  conclude  that  we  can  use  a 
simple  Drude  model  to  predict  the  dielectric  constant  of  n-type  silicon 
over  a  wide  resistivity  range  at  a  wavelength  of  3mm. 
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5.  DoD  AND  OTHER  INTERACTIONS 

Measurements  on  silicon  similar  to  the  ones  described  here,  but  in 
the  far  infrared  region  (2-50  pm),  were  made  in  collaboration  with  S.P. 
Weeks  of  the  Bell  Telephone  Laboratories.  Their  interest  was  in  devel¬ 
oping  a  reliable  technique  for  thin  epitaxial  silicon  layers. 
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C.  ACOUSTOELECTRIC  AND  ACOUSTOOPTIC  DEVICES 


Professors  W.C.  Wang  and  H.  Schachter 

Unit  SSI-3 

1.  OBJECTIVE(S) 

The  objectives  are: 

(a)  To  study  the  physical  phenomena  underlying  the  performance 
of  surface  acoustic  wave  devices  employing  nonlinear  effects  in  semicon¬ 
ductors,  and  in  the  process  to  improve  device  performance  or  invent 
new  devices. 

(b)  To  examine  a  new  method  for  diagnosing  the  surface  and 
interfacial  conditions  of  a  thin-film  semiconductor  by  acoustoelectric 
means . 

(c)  To  study  thin  film  deposition  processes  and  to  develop  ZnO 
and  InSb  thin  films  of  high  quality  for  acoustoelectric  and  acoustooptic 
applications . 

2.  APPROACH 

(a)  As  part  of  the  study  of  physical  phenomena  underlying  the 
performance  of  SAW  devices,  a  new  type  of  FM  demodulator  has  been 
developed  based  on  the  strong  space  charge  nonlinearity  induced  by 
SAW's  in  an  adjacent  semiconductor  and  on  the  filtering  action  genera¬ 
ted  by  spatial  integration  over  the  length  of  the  nonlinear  interaction. 
This  type  of  FM  demodulator  is  structurally  very  simple  and  can  be 
operated  at  high  carrier  frequencies;  its  performance  is  presently  under 
study . 

(b)  A  theory  based  on  the  nonlinear  interaction  between  induced 
charges  and  fields  has  been  developed  in  order  to  describe  the  opera¬ 
tion  of  a  semiconductor  thin-film  convolver  at  flat  band.  In  the  past, 
various  convolver  theories  have  been  proposed  and  presented,  but 
those  theories  were  valid  only  for  semiconductors  of  thickness  large  in 
comparison  with  the  acoustic  wavelength.  One  of  the  reasons  for  devel¬ 
oping  such  a  thin-film  theory  is  that  it  will  help  in  the  diagnoses  of 
thin-film  surfaces  and  interfacial  properties.  Using  SAW  as  a  diagnostic 
tool,  the  advantages  are  (i)  the  depth  of  electric  field  penetration  can 
be  controlled  by  varying  the  sonic  frequency;  thus,  the  interfacial 
properties  of  the  thin  film  can  be  investigated,  (ii)  the  diagnostic 
process  is  contactless  and  nondestructive,  and  (iii)  the  electric  field 
associated  with  the  SAW's  has  both  tangential  and  normal  components, 
so  that  some  of  the  thin-film  orientational  properties  can  be  revealed. 
Computer  results  for  thin-film  convolvers  at  flat  band  have  been  ob¬ 
tained.  Theories  and  models  for  the  non-flat  band  case  will  be  devel¬ 
oped.  The  results  will  be  used  to  evaluate  the  experimental  outcomes. 
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(c)  We  are  developing  a  new  method  for  depositing  ZnO  films  of 
high  quality  using  a  specially-modified  sputter  S-gun  system.  Using 
this  new  method,  a  systematic  study  will  be  made  of  such  films  de¬ 
posited  on  various  substrates  at  various  temperature  levels  to  optimize 
the  technique  and  the  conditions  of  operation.  We  are  also  studying 
techniques  for  improving  the  deposition  of  InSb  films  and  avoiding  a 
layer  of  Sb-InSb-In.  In  this  connection,  we  are  examining  a  flash 
evaporation  scheme  with  a  special  preheating  procedure. 

3.  PROGRESS 

During  the  past  three  years  we  have  made  the  following  major 
accomplishments . 

(a)  We  have  successfully  designed  and  constructed  a  novel  and 
superior  thin  film  deposition  system  utilizing  a  modified  sputter  gun. 
This  new  system  permits  us  to  fabricate  thin  piezoelectric  or  semicon¬ 
ducting  films  of  exceptionally  high  quality  at  both  low  temperature  and 
low  pressure.  The  new  system  has  many  addi  tionaFTd  vantages :  it  is 
cost  effective  (low  cost  both  to  construct  and  to  operate),  it  has  a  high 
deposition  rate  (comparable  to  that  of  the  more  costly  and  less  advan¬ 
tageous  planar  magnetron  system),  it  reduces  greatly  the  problem  of 
thermal  stresses,  and  it  permits  the  fabrication  of  clear,  well-oriented, 
high  quality  films  on  many  different  substrates,  in  contrast  to  other 
systems . 

(b)  Using  our  newly-developed  sputtering  system,  we  have  fabri¬ 
cated  excellent  thin  films  of  ZnO  and  AIN  on  various  substrates.  Good 
InSb  films  have  also  been  made  using  evaporation  techniques.  The  ZnO 
films,  which  came  first,  are  as  good  as  any  made  elsewhere,  to  our 
knowledge,  and  have  received  high  praise  from  industry.  We  have  also 
fabricated  them  on  glass,  silicon  and  quartz  substrates,  with  equal 
facility  and  quality.  The  low  pressure  associated  with  our  sputtering 
system  permits  us  to  produce  AIN  films  of  exceptional  quality,  colorless 
and  with  an  extremely  smooth  surface.  The  InSb  films  we  have  pro¬ 
duced  are  not  yet  of  top  quality,  but  to  our  knowledge  they  are  the 
best  in  the  world  outside  of  the  Soviet  Union,  and  we  have  achieved 
bidirectional  SAW  amplification  with  them,  which  indicates  that  they  are 
already  quite  good. 

(c)  Because  of  discrepancies  which  consistently  appeared  in 
experimental  data  involving  the  acoustoelectric  (AE)  effect,  we  theoreti¬ 
cally  examined  the  Weinreich  relation,  which  relates  the  acoustoelectric 
current  to  the  acoustic  power  loss  and  the  carrier  mobility.  A  two- 
dimensional  analysis  was  carried  out  on  a  configuration  consisting  of  a 
semiconductor  film  on  a  SAW  substrate,  taking  into  account  both  the 
transverse  and  tangential  dc  acoustoelectric  effects.  It  was  found  that 
the  Weinreich  relation,  which  is  widely  used  and  assumed  to  be  correct, 
is  not  valid  most  of  the  time  it  is  used  in  SAW  configurations.  Its 
range  of  validity  is  difficult  to  express  simply,  but  we  are  currently 
clarifying  it. 
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(d)  Since  monolithic  thin  film  semiconductor  convolvers  are  of 
substantial  interest  and  should  become  feasible  now  that  very  good  thin 
films  can  be  fabricated,  we  examined  the  properties  of  several  possible 
types  of  such  convolvers.  The  stress  in  the  study  was  on  the  semi¬ 
conductor  film  on  a  piezoelectric  substrate,  which  should  have  a  time- 
bandwidth  product  comparable  to  that  of  the  successful  but  inconvenient 
separated-medium  type.  We  made  calculations  on  the  dispersion,  atten¬ 
uation,  and  efficiency,  as  a  function  of  film  thickness,  conductivity  and 
frequency,  for  two  types  of  output  terminal  arrangement.  These  anal¬ 
yses  involved  new  considerations,  not  taken  into  account  by  available 
theories . 

These  accomplishments  have  led  to  a  series  of  publications  and 
talks,  the  most  recent  of  which  are  presented  in  Section  4. 

Some  of  the  details  regarding  our  progress  and  accomplishments 
are  briefly  summarized  below. 

A.  ZnO  Films  Grown  at  Low  Temperature  by  Our  Newly  Developed 
"Modified  Sputter  Gun"  System 

The  deposition  of  C-axis  oriented  ZnO  films  by  a  sputtering  tech¬ 
nique  was  initiated  in  the  mid-1960s.  Since  then,  ZnO  films  have  been 
considered  to  be  one  of  the  most  promising  materials  for  optical  wave- 
guiding,  acoustoelectric  applications,  acoustooptic  interaction  media,  and 
SAW  transducers,  based  on  its  high  electromechanical  coupling  coeffi¬ 
cient  and  its  superior  optical  qualities.  However,  ZnO  film  did  not  gain 
momentum  until  the  group  at  Kyoto  University,  Japan,  reported  several 
years  ago  that  high  quality  ZnO  films  can  be  obtained  on  7059  glass 
substrates  at  a  relatively  high  growth  rate  (=  2pm/hr)  by  utilizing 
planar  magnetron  sputtering.  (The  thermal  expansion  coefficients  of 
ZnO  and  7059  glass  are  nearly  identical.)  It  is  of  interest  that  Dr.  S. 
Onishi,  one  of  the  collaborators  on  this  work  unit,  was  a  key  member  of 
that  group  under  Professor  Shiosaki,  and  that  the  work  on  planar 
magnetron  sputtering  was  the  topic  of  his  doctoral  dissertation. 

However,  the  films  on  other  substrates,  such  as  silicon,  quartz, 
etc . ,  are  normally  not  optically  transparent  due  to  their  mismatch  in 
thermal  expansion  coefficients.  In  order  to  improve  the  quality  of  ZnO 
films  on  all  substrates,  we  studied  these  past  techniques  and  recognized 
that  in  order  to  achieve  clear  films  one  must  essentially  avoid  the 
problem  of  thermal  stress  caused  by  electron  bombardment.  We  then 
found  how  to  avoid  such  stress  by  suitably  modifying  the  "research 
S-gun"  sputtering  system.  Presently,  this  modified  S-gun  provides  us 
with  films  of  the  highest  quality.  The  reason  why  will  become  clear  as 
we  review  the  past  sputtering  methods. 

(1)  Using  a  conventional  sputtering+unit:  the  substrate  holder  is 
of  positive  polarity.  While  the  ionized  A  ions  strike  the  targets  (ZnO 
or  Zn)  to  remove  material  for  deposition,  the  electrons  are  bombarding 
the  substrate  since  it  is  at  positive  potential.  As  a  result,  the  sub¬ 
strate  temperature  is  arbitrarily  high  and  can  not  be  controlled.  Films 
thus  obtained  are  of  poor  quality. 
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(2)  Utilizing  a  planar  magnetron  sputtering  system:  With  this 
syo-an,  the  majority  of  the  electrons  are  trapped  in  the  magnetic  field. 
Nevertheless,  the  substrate  will  still  be  bombarded  by  some  of  the 
electrons,  since  the  substrate  holder  is  at  positive  potential.  As  stated 
earlier,  the  planar  magnetron  system  has  significantly  improved  the  film 
deposition  technique.  However,  ZnO  films  which  are  well  oriented, 
clear  and  of  high  quality  are  usually  difficult  to  obtain  on  other  than 
the  7059  glass  substrate. 

With  the  so-called  "research  S-gun"  sputtering  system,  it  was 
necessary  for  us  to  augment  and  specially  arrange  the  magnetic  field  so 
that  the  thin-film  surface  is  completely  outside  of  the  plasma  column  and 
free  from  electron  bombardment.  When  these  modifications  have  been 
made,  the  S-gun  system  is  the  system  that  has  provided  us  with  the 
best  quality  of  film  on  various  substrates,  for  the  following  reasons: 

(a)  The  substrate  holder  is  at  neutral  potential,  and  the 
substrate  surface  is  free  from  electron  bombardment. 

(b)  The  water  cooling  system  is  very  effective,  since  it  is 
right  next  to  the  target.  The  target  will  therefore  be  at 
a  lower  temperature  and  the  substrate  surface  will  not 
be  subjected  to  excessive  heat  radiation  from  the  target. 

(c)  The  magnetic  field  is  several  times  stronger  than  that 
associated  with  a  planar  magnetron.  The  plasma  column 
is  therefore  more  tightly  confined. 

(d)  The  deposition  rate  of  an  S-gun  system  is  comparable  to 
that  of  a  planar  magnetron,  whose  deposition  rate  is 
good. 

Therefore,  the  substrate  holder  is  at  a  low  temperature,  around  150°C, 
rather  than  about  500°C  as  with  the  planar  magnetron.  This  feature  is 
of  key  significance,  since  it  greatly  reduces  the  problem  of  thermal 
stresses.  As  a  result,  with  this  system,  we  are  able  to  obtain  clear, 
well-oriented,  high  quality  ZnO  films  on  aU  substrates  available  to  us, 
such  as  glass,  silicon  and  quartz  substrates. 

In  addition  to  the  key  features  of  low  pressure  and  low  tempera¬ 
ture  sputtering,  other  advantages  of  this  newly- developed  sputter  gun 
system  include:  (a)  a  much  larger  area  of  high  quality  thin  film  can  be 
obtained,  and  (b)  the  system  is  cheaper  to  construct  and  cheaper  to 
run. 


B.  AIN  Films  Grown  at  Low  Pressure  and  Low  Temperature,  Util¬ 
izing  Qur  "Modified  Sputter  Gun"  System 

.  The  utilization  of  our  modified  sputter  gun  has  also  been  extended 
for  the  deposition  of  aluminum  nitride  (AIN)  films  on  various  sub¬ 
strates.  The  strong  magnetic  field  associated  with  this  gun  system 
makes  it  possible  to  maintain  a  stable  plasma  at  a  total  sputtering  gas 
pressure  (TSGP)  as  low  as  1  pTorr.  This  low-pressure  sputtering  pro- 
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vides  films  with  denser  and  finer  grains  and  it  improves  c-orientatior , 
surface  smoothness,  film  color,  transparency  and  deposition  rate.  Our 
experience  shows  that  the  qualities  of  AIN  films  are  extremely  sensitive 
to  the  level  of  total  sputtering  gas  pressure,  the  lower  the  better. 
The  total  pressure  in  our  system  varies  from  1  pTorr  to  8  pTorr.  The 
deposition  rate  of  the  sputter  gun  is  three  times  as  high  as  that  of  a 
conventional  planar  magnetron,  since  the  target  cooling  system  makes  it 
possible  to  increase  rf  or  dc  power.  It  is  2.4  pm/hr.  at  an  input 
power  of  300  watts.  The  film  is  transparent  and  colorless  only  at  low 
pressure  sputtering.  High  quality  films  are  deposited  at  temperatures 
around  350°C,  which  is  low  for  AIN  deposition.  The  substrate  surface 
in  this  system,  as  stated  in  A,  is  free  from  charged  particle  bombard¬ 
ment  and  excessive  heat  radiation. 

C.  InSb  Films 

The  deposition  of  InSb  thin  films  on  LiNb03  substrates  has  been 
accomplished  at  different  centers  since  the  early  1970s,  and  the  process 
has  been  usi  a  with  little  variation  in  Japan,  France,  Russia  and  the 
U.S.  However,  due  to  difficulties  in  obtaining  high  quality  films, 
acoustoelectric  devices  using  InSb  films  have  not  been  proven  to  be 
practical.  For  AE  applications  it  is  required  that  the  InSb  film  pos¬ 
sesses  high  drift  mobility  and  high  resistance.  The  best  data  reported 
to  <iate  are  by  Kotelyanskii ,  et  al,  in  Russia;  they  can  obtain  a  film  of 
500A  thickness  and  a  drift  mobility  of  1600  cm2/v-sec.  This  mobility  is 
over  three  times  better  than  what  has  been  reported  in  the  western 
world.  We  have  been  engaging  in  the  development  of  InSb  films  for 
three  years  and  have  achieved  bidirectional  SAW  amplification  with  those 
films  (necessary  for  reducing  the  insertion  loss  of  a  convolver).  The 
drift  ^nobility  we  obtained  here  is  =  700  cm2/v-sec  at  a  film  thickness 
-  500A,  which  is  better  than  what  is  reported  in  this  country  but  is 
still  low  in  comparison  with  that  of  Kotelyanskii,  et  al.  Private  discus¬ 
sions  with  Kotelyanskii  during  his  recent  visit  to  us  here  revealed  that 
their  set-up  is  much  more  elaborate  and  costly,  and  they  have  invested 
10  years  in  this  effort. 

We  have  measured  the  Hall  mobility  as  a  function  of  the  conduc¬ 
tivity  thickness  product.  The  film  thickness  is  about  450A.  Our 
process,  especially  on  forming  the  electric  contacts,  we  feel  is  superior. 
The  highest  Hall  mobility  obtained  by  us  is  about  2500  cm/volt-sec, 
which  means  that  if  the  traps  can  be  eliminated  in  this  film  a  drift 
mobility  of  that  same  magnitude  can  be  achieved.  The  significance  of 
these  measurements  is  that  they  tell  us  how  good  a  film  can  be  achieved 
ultimately,  when  the  Hall  and  drift  mobilities  are  equal.  Thus,  these 
measurements  offer  us  substantial  hope. 

D.  Examining  the  dc  Acoustoelectric  Current  Produced  by  SAW 

•  It  is  generally  believed  that  surface  acoustic  waves  will  be  useful 
in  determining  the  surface  and  interfacial  properties  of  semiconductor 
wafers  and  semiconductor  films;  however,  in  the  process  of  our  corre¬ 
lating  data  from  SAW  measurements  with  those  obtained  from  conven¬ 
tional  methods,  some  discrepancies  consistently  appear.  This  led  us  to 
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examine  the  validity  of  the  commonly- used  Weinreich  relationship  for  a 
configuration  consisting  of  a  semiconductor  on  a  SAW  substrate.  The 
Weinreich  formula  relates  the  acoustoelectric  current  to  the  acoustic 
power  loss  and  the  carrier  mobility.  Since  both  the  AE  current  and  the 
power  loss  are  measurable  quantities,  the  carrier  mobility  can  thus  be 
determined.  Comparing  the  mobility  measured  here  with  the  Hall  mo¬ 
bility  measurement,  trapping  dynamics  can  also  be  revealed. 

The  Weinreich  relation  was  first  derived  for  phonon-drag  phenom¬ 
ena  in  Ge  due  to  the  deformation  potential.  Its  validity  in  bulk  piezo¬ 
electric  semiconductors  and  in  the  presence  of  acoustic  gain  was  experi¬ 
mentally  verified  by  one  of  us  (W.C.  Wang)  in  1962,  and  subsequently 
analyzed  by  H.N.  Spector.  In  examining  the  Weinreich  formula  for 
thin  films,  however,  we  have  carried  out  a  two-dimensional  analysis  for 
a  configuration  consisting  of  a  semiconductor  film  on  a  SAW  substrate, 
taking  into  account  both  the  transverse  and  tangential  dc  AE  effects. 
The  analysis  shows  that  the  range  of  validity  of  the  Weinreich  relation 
in  such  a  system  is  very  limited  and  that  in  general  it  cannot  be  ap¬ 
plied  here.  (Its  range  of  validity  cannot  be  simply  stated,  but  we  are 
currently  assessing  the  validity  conditions . )  In  addition ,  the  analysis 
reveals  that  (a)  the  field  associated  with  the  transverse  AE  current  in 
general  tends  to  accumulate  the  charge  at  the  semiconductor  surface, 
(b)  the  polarity  of  the  transverse  AE  current  will,  in  general,  not 
change  sign  when  the  film  thickness  is  larger  than  the  SAW  wavelength, 
and  (c)  the  diffusion  effect  cannot  be  ignored  in  the  calculation  of 
tangential  AE  current. 

E.  Monolithic  Thin  Film  Semiconductor  Convolver 

Many  forms  of  convolver  have  been  designed  and  tested.  The  best 
performance  to  date  is  the  one  given  by  the  separated-medium  Si  on 
LiNb03  normal-component  convolver  fabricated  at  Lincoln  Laboratory  and 
Texas  Instruments.  However,  looking  ahead  and  taking  into  considera¬ 
tion  eventual  cost  effectiveness  and  smaller  size,  the  monolithic  con¬ 
volver  still  remains  a  strong  candidate.  There  are  three  possible 
monolithic  AE  convolver  structures:  (i)  piezoelectric  film  on  semicon¬ 
ductor  substrate;  (ii)  semiconductor  film  on  piezoelectric  substrate;  (iii) 
piezoelectric  film  on  semiconductor  film.  The  time  bandwidth  (TB) 
product  for  case  (i)  is  usually  rather  low  due  to  the  dispersion  intro¬ 
duced  by  the  piezoelectric  film  (operating  in  the  Sezawa  mode).  The 
TB  product  for  case  (ii)  should  be  comparable  to  that  of  the  separated 
medium  type.  Therefore,  we  started  to  analyze  its  performance.  Case 
(iii)  will  be  investigated  in  the  future. 

Two  types  of  output  terminal  arrangement  for  case  (ii)  have  been 
investigated.  In  one  type,  the  output  terminal  plate  is  in  direct  ohmic 
contact  with  the  top  of  the  semiconductor  surface,  and  in  the  other  an 
insulating  dielectric  layer  is  inserted  between  the  semiconductor  surface 
and  . the  convolver  on  the  output  plate.  Therefore,  the  former  has  a 
short-circuit  condition  and  the  latter  an  open-circuit  condition  on  the 
semiconductor  surface.  Since  the  boundary  conditions  are  different, 
the  velocity,  attenuation  constant  and  convolver  efficiency  are  found  to 
be  quite  different  in  the  two  configurations.  Calculations  on  the  dis- 
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persion,  attenuation  and  efficiency  have  been  carried  out  as  a  function 
of  film  thickness,  conductivity  and  frequency. 

4.  PUBLICATIONS 

1.  "Transparent  and  Highly  Oriented  ZnO  Films  Grown  at  Low 
Temperature  by  Sputtering  with  a  Modified  Sputter  Gun,"  by 

S.  Onishi,  M.  Eschwei  and  W.C.  Wang,  Applied  Phys.  Lett., 
vol .  38(6),  p.  319  (1981). 

2.  "Colorless-Transparent,  c-Oriented  Aluminum  Nitride  Films 
Grown  at  Low  Temperature  by  a  Modified  Sputter  Gun,"  by 
S.  Onishi,  M.  Eschwei,  S.  Bielaczy,  and  W.C.  Wang,  sched¬ 
uled  to  appear  in  October  15,  1981  issue  of  Applied  Phys. 
Lett. 

3.  "Properties  of  InSb  Film  for  SAW  Amplifier,"  by  M.  Eschwei, 
S.  Onishi  and  W.C.  Wang,  accepted  for  presentation  in 
American  Vacuum  Society  28th  National  Vacuum  Symposium 
(1981). 

4.  "Monolithic  Thin  Film  Acoustoelectric  Convolver,"  by  Z.S.  Wu, 

S.  Onishi,  K.S.  Chen  and  W.C.  Wang,  to  be  presented  in 
1981  Ultrasonics  Symposium. 

5.  "Acoustic,  Optical  Properties  of  AIN  Films  Grown  at  Low 
Temperature  by  Modified  Sputter  Gun,"  by  S.  Onishi,  M. 
Eschwei,  S.  Bielaczy  and  W.C.  Wang,  to  be  presented  in  1981 
Ultrasonics  Symposium. 

6.  "D-C  Acoustoelectric  Current  Induced  in  a  Semiconductor  Film 
Adjacent  to  a  SAW  Substrate,"  by  K.S.  Chen,  Z.S.  Wu,  H. 
Schachter  and  W.C.  Wang,  to  be  presented  in  1981  Ultra¬ 
sonics  Symposium. 

5.  DoD  AND  OTHER  INTERACTIONS 

Industry  has  expressed  substantial  interest  in  our  thin  film  fabri¬ 
cation  capability  and  its  potential  for  producing  excellent  films  in  a 
cost-effective  manner.  Three  examples  are  given: 

(a)  Fairchild  Camera  and  Instrument  is  discussing  with  us  a 
program  which  would  utilize  an  InSb  thin  film  amplifier  in  the  construc¬ 
tion  of  a  compact  and  low  loss  SAW  delay  line  for  repeater-jammer 
applications.  Negotiations  with  Mr.  E.  Magill  are  in  progress. 

(b)  Hazeltine  Corporation  is  discussing  with  us  a  program  to 
develop  ZnO  thin  films  on  quartz  for  wideband  transducers.  Negotia¬ 
tions  with-  Dr.  R.  LaRosa  are  in  progress. 

(c)  The  earliest  public  recognition  of  our  capability  was  made  by 
Dr.  F.  Hickernell  of  Motorola,  Arizona,  with  whom  we  were  in  contact. 
During  his  invited  talk  on  ZnO  thin  films  at  the  1979  IEEE  Ultrasonics 
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Symposium,  he  held  up  one  of  our  high  quality  (transparent  and  color¬ 
less)  ZnO  films  on  glass  and  stated  that  this  is  an  example  of  what  a 
really  good  film  should  be  like. 

In  addition,  we  have  had  visits  and  arrangements  to  visit  us  from 
a  large  variety  of  sources.  Some  examples  are: 

(d)  Drs.  J.M.  Kotelyanskii  and  V.T.  Potapov  of  the  Institute  of 
Radioengineering  and  Electronics,  Academy  of  Sciences  of  the  USSR, 
Moscow,  visited  us  in  Brooklyn  (and  therefore  did  not  see  our  facility, 
which  is  in  Farmingdale) ,  and  we  exchanged  information  on  InSb  films. 
Dr.  Kotelyanskii's  laboratory  produces  the  best  InSb  films  in  the  world. 
It  is  of  interest  that  about  one  month  later,  at  an  URSI  General  Assem¬ 
bly  in  Washington,  DC,  Professor  Oliner  was  approached  by  a  "young 
scientist"  delegate,  Dr.  V.I.  Anisimkin,  from  the  same  laboratory  in  the 
USSR,  who  referred  to  that  meeting  and  asked  if  he  could  spend  up  to 
six  months  in  Professor  Wang's  laboratory. 

(e)  Professors  P.  Kornreich  and  S.  Kowel  of  Syracuse  University 
are  scheduled  to  visit  us  to  discuss  our  AIN  and  ZnO  deposition  tech¬ 
niques. 

(f)  Dr.  H.  Gautier  and  Dr.  P.  Defranould  of  Thomson-CSF, 
Cagnes  sur  Mer,  France,  are  scheduled  to  visit  us  just  after  this  year's 
IEEE  Ultrasonics  Symposium. 

(g)  Among  others  who  have  expressed  interest  in  visiting  our 
facility  are  Mr.  E.  Stern  of  the  MIT  Lincoln  Laboratory  and  Professor 
P.  Das  of  the  Rensselaer  Polytechhnic  Institute. 
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D.  ELECTRIC  AND  MAGNETIC  INTERACTIONS  IN  THIN  FILMS  AND 

SURFACE  REGIONS  OF  SOLIDS 

Professor  H.J.  Juretschke 

Unit  SSI-4 

1.  OBJECTIVE(S) 

To  understand  the  equilibrium  properties,  the  approach  to  equilib¬ 
rium,  and  the  response  to  applied  fields  and  waves  of  the  surface 
regions  of  selected  solids,  both  conductors  and  insulators,  and  to 
develop  new  methods  for  their  investigation. 

2.  APPROACH 

The  experimental  approach  concentrates  on  the  change  of  prop¬ 
erties  of  solids  when  the  surface/volume  ratio  of  their  geometry  becomes 
significant.  Transport  properties  are  examined  for  their  size  effects, 
and  on  the  influence  of  specific  parameters,  such  as  electrostatic  charg¬ 
ing,  changes  of  composition  or  of  order,  or  mechanically  induced  sur¬ 
face  strains,  on  these  size  effects.  Emphasis  is  on  those  special  inter¬ 
actions  that  lead  to  nonlinear  responses  or  to  resonances,  and  on  the 
coupling  between  the  electromagnetic  and  elastic  fields  in  the  surface 
region.  The  theory  will  emphasize  use  of  the  methods  developed  in 
recent  years  for  the  study  of  defects  in  solids  for  the  investigation  of 
the  surface,  considered  as  an  extended  defect. 

3.  PROGRESS 

The  main  emphasis  during  the  past  year  has  been  in  the  two 
areas : 

A.  Charge- induced  surface  stresses  at  metal  surfaces 

B.  Surface  layer  self -diffusion. 

Except  for  item  (1),  below,  our  progress  will  be  summarized  sepa¬ 
rately  for  areas  A  and  B. 

(1)  Work  in  both  areas  has  been  backed  by  the  development  of  a 
new  theoretical  model  of  conduction-electron  surface  scattering,  that  is 
based  on  scattering  in  depth  within  a  surface  layer  of  electrical  prop¬ 
erties  differing  from  those  of  the  underlying  material.  This  theory, 
incorporating  effects  of  thermal  surface  phonons  as  well  as  of  geo¬ 
metrical  rearrangements,  is  now  ready  for  publication.  Its  predictions 
both  overlap  and  differ  from  those  of  other  treatments  of  surface  scat¬ 
tering.  Thus,  as  expected,  surface  scattering  becomes  more  important 
as  either  the  sample  thickness  decreases,  or  the  carrier  mean  free  path 
increases,  but  the  universal  behavior  in  terms  of  the  ratio  of  these  two 
lengths  characteristic  of  most  other  models  is  no  longer  true.  In  addi¬ 
tion,  surface  phonons  may  predict  spurious  asymptotic  values,  and 
temperature- independent  surface  structures  may  simulate  temperature- 
dependent  surface  scattering,  if  interpreted  in  terms  of  those  models. 

A.  Charge-induced  surface  scattering 

(2)  The  effect  of  surface  charging  on  electron  scattering  at 
well-annealed  metal-vacuum  interfaces  has  now  been  measured  exten¬ 
sively  on  samples  of  Ag,  Au,  and  Cu,  in  the  range  1.8  <  T  <  300°K. 
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The  data,  though  individually  ranging  over  factors  larger  than  100,  are 
interpreted  successfully  in  terms  of  the  nearly  universal  scaled  behavior 
predicted  by  our  theory,  for  values  of  the  scaling  parameter  K  (=  sam¬ 
ple  thickness/bulk  mean  free  path)  differing  by  as  much  as  103.  The 
temperature-independent  surface  charge  density  dependence  of  the 
surface  strain  is  about  0.1/(C/m2),  for  Ag. 

(3)  Preliminary  data  showing  anomalies  at  very  low  temperatures 
in  surface  scattering  at  very  low  temperatures  have  not  been  confirmed, 
largely  because  of  the  uncertainty  in  the  bulk  properties  needed  for 
reference ,  which ,  in  this  temperature  range ,  are  very  sensitive  to 
purity  and  sample  preparation. 

(4)  The  mean  penetration  depth  of  surface  strains  at  these  well- 
annealed  surfaces  has  been  determined  to  be  about  10  to  15  A.  Charged- 
induced  scattering,  which  alters  this  depth  by  a  very  small  amount 
(~  10~5)  can  be  suppressed  if  the  metal  surface  is  covered  by  a  layer 
of  graphite  of  this  order  of  thickness.  Graphite  has  been  found  to 
present  less  problems  as  a  very  thin  layer  than  the  previously  at¬ 
tempted  gold. 

(5)  For  both  Ag  and  Au,  positive  charging  increases  the  surface 
strain,  in  accord  with  the  increase  in  perpendicular  lattice  constant 
observed  on  most  metal  surfaces.  The  opposite  sign,  seen  in  Cu,  has 
been  found  to  arise  from  the  rapid  oxidation  of  any  exposed  copper 
surface,  and  reflects  the  influence  of  adsorbed  gases  on  surface 
strains. 

(6)  The  strains  induced  in  all  samples  by  Maxwell  stresses  (i.e., 
proportional  to  q2,  where  q  is  the  surface  charge  per  unit  area)  also 
produce  changes  in  (volume)  scattering,  in  the  same  range  of  tempera¬ 
tures  for  which  we  have  studied  surface  strain  scattering.  This  volume 
scattering  shows  first  an  anomalously  large  increase  as  the  temperature 
is  lowered,  and  is  then  followed  by  a  precipitous  drop,  often  to  com¬ 
plete  vanishing,  at  temperatures  below  about  7°K.  Both  aspects  of  this 
behavior  are  unexpected,  and,  so  far  unexplained.  One  obstacle  to 
understanding  this  is  the  complete  absence  of  experimental  or  theoreti¬ 
cal  information  on  the  strain  sensitivity  of  the  bulk  conductivity  of 
metals  at  low  temperatures.  Another  is  the  uncertainty  in  the  stress 
distribution  and  therefore  the  strain  distribution  caused  by  finite  sam¬ 
ple  size,  which  is  very  sensitive  to  sample  geometry.  An  additional 
possibility  that  remains  to  be  explored  is  that  surface  strains  also 
contain  a  contribution  quadratic  in  q  at  very  low  temperatures. 

(7)  Charge-modulated  surface  scattering  at  the  Ag-Mica  interface 
follows  roughly  the  same  behavior  as  at  the  free  surface,  for  T  >  60°K. 
However,  preliminary  results  below  this  temperature  show  that  surface 
strain  scattering  reverses  sign,  and  is  accompanied  by  a  contribution 
quadratic  in  q  peaked  at  the  cross-over,  but  of  opposite  sign  to  that 
induced  by  Maxwell  stresses.  A  similar  effect  had  been  seen  earlier  in 
Au,  but,  occurring  at  a  much  higher  temperature,  it  is  far  less  drama¬ 
tic,  and  far  less  detailed.  If  confirmed,  the  low  temperature  data  for 
the  Ag-Mica  interface  indicate  that  this  interface  is  much  more  complex 
than  the  free  surface,  and  that  it  probably  includes  interfacial  strains 
arising  from  a  temperature-dependent  epitaxial  mismatch  of  the  two 
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adjacent  materials.  If  properly  interpreted,  the  data  should  become  a 
rich  source  of  information  on  such  interfacial  strains. 

B.  Disorder-induced  surface  scattering 

(8)  Experiments  on  the  initial  response  of  surface  scattering  to  a 
disordered  Ag  overlay  deposited  on  a  flat  silver  surface  at  low  tempera¬ 
ture  show  that  surface  scattering  saturates  when  the  layer  thickness 
exceeds  about  20  A,  indicating  that  above  that  thickness  the  overlay 
acts  like  a  perfect  absorber  of  incident  electrons.  Consequently,  the 
use  of  surface  scattering  to  investigate  the  layer  structure  is  most 
sensitive  for  smaller  layer  thicknesses. 

(9)  For  each  surface  layer  structure,  we  have  found  a  tempera¬ 
ture  Tc  below  which  surface  scattering  is  fully  reversible  as  the  temper¬ 
ature  T  moves  up  or  down.  But  once  T  exceeds  T  ,  surface  scattering 

decreases  irreversibly,  first  very  fast,  and  then  more  slowly.  By 
proper  choice  of  the  excess  temperature  and  the  time  interval  during 
which  it  is  maintained  before  the  temperature  is  reduced  again,  a  new 
reversible  temperature  range  for  surface  scattering,  and  a  new  Tc  for 

this  next  stage  of  annealing,  can  be  established.  These  characteristics 
follow  very  much  those  of  a  typical  glass-crystalline  transition. 

(10)  By  dividing  this  process  of  partial  annealing  into  about  ten 
stages,  and  then  analyzing  the  reversible  region  below  each  Tc  in  terms 

of  our  theory  (item  1,  above),  we  have  been  able  to  interpret  the 
scattering  data  in  terms  of  an  effective  disturbed  layer  thickness  and 
an  effective  disorder  parameter  within  that  layer,  as  a  function  of  T  . 

The  first  results  of  such  analysis  show  that  the  layer  thickness  first 
reduces  slowly  and  then  more  rapidly,  and  stabilizes  at  between  10  and 
15  A,  and  that  the  disorder  parameter  first  decreases  rapidly,  and  then 
more  slowly  as  annealing  proceeds,  and  also  stabilizes  at  a  finite  value. 
Both  trends  are  in  good  agreement  with  the  numbers  obtained,  in  a 
completely  different  experimental  manner,  under  item  4  above.  How¬ 
ever,  other  implications  of  the  theory  have  to  be  confirmed  before  this 
interpretation  can  be  fully  accepted. 

(11)  We  have  been  able  to  establish  satisfactory  experimental 
conditions  for  following  the  change  of  surface  scattering  as  a  function 
of  time  at  a  given  T,  and  are  thus  able  to  determine  the  effective 
activation  energies  connecting  any  two  stages  of  annealing.  As  ex¬ 
pected,  the  activation  energy  increases  as  the  annealing  approaches 
completion . 

(12)  In  the  final  stages  of  annealing,  Tc  increases  rapidly,  though 

still -stayirig  well  below  room  temperature,  indicating  that  there  is  an 
extended  tail  of  high  activation  energies  that  have  to  be  overcome 
before  the  surface  is  stable  and  reversible  under  changes  of  tempera¬ 
ture.  We  hope  to  be  able  to  establish  by  a  more  detailed  study  whether 
these  high  energy  sites  belong  to  the  outer  surface,  or  to  the  overlay- 
base  interface. 
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Items  2  to  7  represent  work  contained  in  a  Ph.D.  thesis  to  be 
completed  by  September  1981.*  Its  results  will  appear  in  a  number  of 
papers.  Items  8  to  12  are  part  of  another  Ph.D.  thesis  currently  in 
progress.** 

4.  PUBLICATIONS 

1.  D.J.  Lischner  and  H.J.  Juretschke,  "Maxwell  Stresses  at 
Charged  Metal  Surfaces  from  Thin-Film  Elastoresistance , "  J. 
Appl.  Phys.  51.  474  (1980). 

2.  H.J.  Juretschke  and  R.  Pimpinella,  "Size  Effects  in  Thin 
Metallic  Films  Resulting  from  Modified  Surface  Transport 
Properties,"  submitted  to  Surface  Science. 

3.  R.  Pimpinella  and  H.J.  Juretschke,  "Low  Temperature  Anoma¬ 
lies  in  the  Surface  Charge  Modulation  of  Thin  Metallic  Film 
Conductance,"  to  be  published. 

4.  R.  Pimpinella  and  H.J.  Juretschke,  "Thin  Film  Surface  Struc¬ 
ture  Determined  by  Surface  Scattering  of  Conduction  Elec¬ 
trons,"  to  be  published. 

5.  Y.  Amani  and  H.J.  Juretschke,  "Transport  in  Thin  Films  with 
Rough  Overlays,"  to  be  published. 

5.  DoD  AND  OTHER  INTERACTIONS. 

In  view  of  our  extensive  experience  with  the  physics  of  strains  on 
metal  surfaces,  gained  largely  under  JSEP  support  on  this  work  unit 
and  its  predecessor,  we  were  asked  by  AFOSR  about  three  or  four 
years  ago  to  advise  it  on  the  RAD  AM  problem,  involving  the  modulation 
of  radar  scattering  by  moving  or  vibrating  contacts.  This  consultation 
led  to  a  separate  contract  with  AFOSR  to  determine  the  effects  due  to 
the  intrinsic  nonlinear  electrical  properties  of  weathered  and  rusty  metal 
surfaces.  In  particular,  we  investigated  the  role  of  microscopic  metal- 
oxide  barriers .  We  found  that  they  allow  tunneling  of  electrons 
through  surprisingly  thick  barriers,  and  at  low  energies,  and  that  the 
sensitivity  of  such  barriers  to  macroscopic  strain  is  large,  but  not  quite 
large  enough  to  explain  the  modulation  effects  on  that  basis  alone. 


*  R.  Pimpinella,  "Metallic  Field  Effect  and  Surface  Scattering  at  Free 
-Metal  Surfaces." 

**  Y.  Amani,  "Equilibrium  and  Non-equilibrium  Structures  of  Silver 
Surfaces . " 
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E.  NEW  SOLID  STATE  MATERIALS 


Professor  E.  Banks 

Unit  SSI-5 

1.  OBJECTIVE(S) 

To  synthesize,  characterize  and  measure  the  physical  properties  of 
new  solid  state  materials.  Such  materials  are  potentially  ferroelectric 
(and  optically  nonlinear),  ferromagnetic,  fluorescent  or  semiconducting, 
all  of  which  have  considerable  technical  interest  and  value  if  the  prop¬ 
erties  investigated  can  be  applied  to  practical  devices.  Even  when  they 
are  not  of  immediate  practical  value,  the  study  of  their  properties  can 
lead  to  understanding  that  may  point  the  way  to  the  design  of  improved 
materials . 

2.  APPROACH 

Attempts  to  synthesize  new  materials  with  interesting  electronic 
properties  are  based  either  on  modifying  existing  materials  (e.g.,  by 
forming  solid  solutions)  or  by  attempting  the  synthesis  of  materials 
based  on  theories  which  attempt  to  explain  the  dependence  of  the  rele¬ 
vant  property,  e.g.,  ferromagnetism,  on  features  of  crystal  structure 
and  chemical  bonding.  New  materials  are  characterized  by  crystallo¬ 
graphic  methods  and  measurement  of  relevant  electronic  properties  such 
as  optical  spectroscopy,  fluorescence,  magnetic  properties,  Mossbauer 
.spectra,  conductivity  and  dielectric  properties. 

3.  PROGRESS 

(a)  Complex  Fluorides  of  Bivalent  Rare  Earths:  Potential  Laser 
and  Scintillation  Counters 


Experiments  on  crystal  growth  of  the  new  compounds,  EuMgF4  and 
SrMgF4,  failed  previously  because  the  compounds  melt  incongruently , 
decomposing  to  the  simple  fluorides.  During  the  past  year,  however, 
we  have  found  that  mixtures  of  SrMgF4  with  BaMgF4  can  be  melted 
without  decomposition  up  to  70  mole  %  SrMgF4.  Boules  of  these  mix- 
tures  have  been  found  to  have  a  uniform  composition.  This  means  that, 
under  suitable  conditions,  single  crystals  of  such  solid  solutions  should 
be  capable  of  growth.  Polycrystalline  boules  of  90%  BaMgF4,  10% 
EuMgF4  have  been  observed  to  fluoresce  bright  blue  under  x-ray  exci¬ 
tation.  Since  the  optimum  brightness  in  the  SrMgF4-EuMgF4  system 
occurs  at  75%  EuMgF4,  we  can  expect  more  efficient  x-ray  detection  at 
compositions  from  50-74%  EuMgF4  in  the  BaMgF4  system.  With  such  a 
high  concentration  of  heavy  elements  (Ba,  Eu)  the  x-ray  sensitivity 
should  be  very  high.  Crystals  of  such  properties  should  be  very 
useful  detectors  for  x-rays  and  other  ionizing  radiation.  They  would 
have  the  advantage  over  current  crystals,  based  on  Csl,  of  being 
resistant  to  moisture,  and  even  very  small  crystals  would  be  useful  for 
array  devices  and  medical  applications,  such  as  CAT  scanners. 
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An  x-ray  structure  analysis  of  the  three  new  compounds  has  been 
completed,  and  they  are  all  centrosymmetric ,  unlike  BaMgF4  which  is 
non-centro  and  hence  piezoelectric  at  room  temperature.  The  solid 
solutions  should  have  composition-dependent  ferroelectric  transition 
temperatures,  creating  the  possibility  of  preparing  compositions  whose 
fluorescence  can  be  modulated  by  applied  electric  fields. 

In  the  continuation  of  this  program,  major  emphasis  will  be  placed 

on  the  single  crystal  growth  problem,  as  this  is  crucial  to  all  experi¬ 

mental  work,  both  of  the  basic  and  applied  sort.  The  concentration  will 
be  on  direct  melting  and  Bridgman  crystal  growing  experiments  on  the 
solid  solutions  of  the  BaMgF4-SrMgF4  system,  and  the  BaMgF4-EuMgF4 
system,  focusing  on  the  higher  (50  mole  percent  and  above)  concen¬ 
trations  of  Sr  and  Eu  compounds. 

Initial  studies  of  these  crystals  will  deal  with  the  transition  from 

non-centric  to  centric  crystal  symmetry  as  a  function  of  composition  and 

temperature,  followed  by  studies  of  ultraviolet- excited  fluorescence  of 
the  Eu- doped  materials  over  the  range  from  cryogenic  temperatures  to 
the  temperature  of  fluorescence  quenching  (expected  to  be  above  room 
temperature).  If  the  expected  ferroelectric  phase  transitions  are  de¬ 
tected,  studies  of  the  effect  of  such  transitions  on  the  luminescence  will 
be  made.  Other  experimental  measurements  such  as  19F  NMR  studies, 

1SlEu  MSssbauer  studies  and  Eu  EPR  studies,  will  be  taken  to  obtain 
a  more  detailed  understanding  of  the  nature  of  the  transitions  and  their 

2+ 

effect  on  the  local  environment  of  the  Eu  ion  in  the  crystals. 

2+ 

If  additional  time  is  available,  we  plan  to  add  Sm  to  the  crystals 
in  various  concentrations.  The  purpose  is  to  study  the  line  spectrum 
2+ 

of  Sm  as  a  function  of  relative  concentration  of  Sm  and  Eu  and  to 

2+  2+ 

study  the  interesting  energy  transfer  from  Eu  to  Sm  in  these  mater- 

2+ 

ials.  The  dependence  of  the  fluorescent  lifetimes  of  the  donor  Eu  ion 

2+ 

and  the  acceptor  Sm  ion  on  concentration  will  provide  evidence  of  the 
rates  of  energy  transfer  in  these  materials.  In  addition  to  the  basic 
scientific  interest,  these  measurements  will  provide  information  as  to  the 
feasibility  of  lasers  prepared  from  these  materials.  Some  of  the  energy 
transfer  experiments  can  be  done  on  powdered  materials  and  some 
samples  have  already  been  prepared.  If  time  permits,  we  also  expect  to 
conduct  tests  of  x-ray  sensitivity  of  these  materials,  and,  if  single 
crystals  are  available,  some  tests  of  scintillation  counter  efficiency  will 
be  undertaken. 


(b)  Attempted  Preparation  of  Ferromagnetic  Oxides 


Attempts  to  prepare  ordered  perovs  kites  containing  equal  numbers 
of  d3  and-  d5  transition  metal  ions,  whose  nearest-neighbor  exchange 
interactions  are  ferromagnetic,  have  continued.  Single  phase  perovskite 
materials  have  been  prepared  in  a  narrow  range  of  compositions  around 

3+ 


4+  3+ 

BaQ  5G<1q  jMn  q  5Fe  o  5°3  anc*  same 
ing  GdJ  ,  have  been  prepared.  Magnetic  susceptibility  studies 


with  Eu  replac- 
show 
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the  Eu  compounds  to  undergo  apparent  antiferromagnetic  ordering  at 
18°K,  but  the  Gd  compounds  behave  like  paramagnetic  materials  down  to 
4.2°K.  Mossbauer  spectroscopy  has  been  unsuccessful,  owing  to  the 
high  absorption  by  Ba  and  the  rare  earths  for  the  14  KeV  gamma  rays 
from  s7Fe.  In  both  series  of  compounds,  long-range  magnetic  ordering 
is  either  absent  or  very  weak.  A  major  difference  between  the  two 
series  of  materials  is  seen  in  the  behavior  of  the  high- temperature 
dependence  of  the  reciprocal  susceptibility-temperature  data.  In  the 
case  of  the  Gd  series,  this  curve  extrapolates  to  zero  at  about  4°K, 
close  to  absolute  zero,  while  the  Eu  compounds  have  an  intercept  at 
120°K,  indicating  strong  antiferromagnetic  interactions  between  nearest 

3+ 

neighbors.  Since  the  magnetic  moment  of  Eu  is  very  small,  this 
suggests  that  the  magnetic  coupling  is  dominated  by  antiferromagnetic 
Mn-Mn  and  Fe-Fe  interactions  over  the  ferromagnetic  Mn-Fe  interac- 

4+  3+ 

tions.  In  short,  these  materials  behave  as  if  the  Mn  and  Fe  ions 
are  disordered  in  the  crystals. 

On  the  other  hand,  the  Gd  series  contain  a  rare  earth  ion  with  a 
high  magnetic  moment  (S  =  7/2).  This  is  close  to  the  sum  of  the  spin 

quantum  numbers  of  Fe3+  and  Mn4+  (|  +  \  ~  §  =  4).  If  the  net  moment 

of  the  octahedral  site  ions  coupled  antiferromagnetically  to  the  Gd  ion 
moment,  and  no  long  range  magnetic  ordering  occurred,  but  the  mag¬ 
netic  moments  became  frozen  in  random  directions,  we  would  observe 
the  behavior  described  above.  This  has  come  to  be  known  as  "spin- 
glass"  behavior.  Such  behavior  is  characterized  by  the  detection  of 
hyperfine  splitting  in  MSssbauer  spectra,  coupled  with  the  absence  of 
long-range  magnetic  ordering  as  observed  by  neutron  diffraction. 

Unfortunately,  both  Gd  and  Eu,  the  rare  earths  that  formed  these 
materials,  have  exceedingly  high  absorption  cross-sections  for  thermal 
neutrons .  The  next  rare  earths ,  Tb  and  Dy ,  have  high  magnetic 
moments,  and  their  nuclei  are  not  strong  neutron  absorbers.  Prelim¬ 
inary  experiments  to  form  perovskite  phases  have  been  unsuccessful. 

We  will  continue  attempting  the  preparation  of 
BaQ  5TbQ  5MnQ  5FeQ  and  its  dysprosium  analog.  Some  effort  will 

be  placed  on  the  use  of  low-melting  fluxes  such  as  Ba-0-B203,  both  for 
crystal  growth  and  for  enhancing  the  tendency  of  the  atoms  to  assume 
an  ordered  arrangement  in  the  structure.  Preparation  of  the  Gd  com¬ 
pound  using  57Fe-enriched  starting  material  will  be  attempted  in  order 
to  obtain  a  sample  suitable  for  Mfissbauer  spectroscopy.  Attempts  will 
be  made  to  prepare  garnet  phases  where  the  octahedral  sites  contain  d3 

3+  3+ 

ions  such  as  Cr  while  the  tetrahedral  sites  contain  Fe  .  Early 

3+ 

attempts  failed  to  get  complete  substitution  of  Cr  in  the  octahedral 

3+ 

sites  of  Y3Fe5012,  apparently  because  the  Cr  ions  are  too  small;  we 
will  attempt  this  in  a  garnet  with  a  smaller  unit  cell,  e.g.,  YbFe5012  or 
LuFe5012. 
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4.  PUBLICATIONS 

1.  E.  Banks,  Y.  Okamoto  and  Y.  Ueba,  "Synthesis  and  Charac¬ 
terization  of  Rare  Earth  Metal-Containing  Polymers.  I. 

Fluorescence  Properties  of  Ionomers,"  J.  Appl.  Polymer  Sci., 
25,  359  (1980). 

2.  Y.  Ueba,  E.  Banks  and  Y.  Okamoto,  "Synthesis  and  Charac¬ 
terization  of  Rare  Earth  Metal -Containing  Polymers.  II. 

3+ 

Fluorescence  Properties  of  Eu  -Polymer  Complexes  Containing 
p-Diketone  Ligand,"  J.  Appl.  Polymer  Sci.,  25,  2007  (1980). 

3.  Y.  Avigal  and  E.  Banks,  "Conductivity  of  Reduced  Solid 

Solutions  in  the  System  Na20-Al203-Ti02,"  Material  Research 
Bulletin  15,  1381  (1980). 

4.  E.  Banks,  S.  Nakajima  and  M.  Shone,  "New  Complex  Fluor¬ 

ides  EuMgF4,  and  Their  Solid  Solutions:  Photoluminescence 
and  Energy  Transfer,"  J.  Electrochem.  Soc.,  127,  2234 
(1980). 

5.  R.  Sacks,  Y.  Avigal  and  E.  Banks,  "New  Solid  Electrolytes 
Based  on  Cubic  ZrP207,"  submitted  to  J.  Electrochem.  Soc. 

6.  E.  Banks,  R.  Jenkins  and  B.  Post,  "Crystal  Structures  of 

EuMgF4,  SmMgF4  and  SrMgF4,"  presented  at  15th  Rare  Earth 
Research  Conference,  Rolla,  MO,  15  June  1981;  manuscript 
submitted  for  publication  in  Conference  Proceedings. 

7.  E.  Banks,  Y.  Ueba  and  Y.  Okamoto,  "Lanthanide  Ion  Fluores¬ 
cence  as  a  Probe  of  Ionomer  Structure,"  Annals  of  N.Y. 
Acad.  Sciences  366,  356  (1981). 

5.  DoD  AND  OTHER  INTERACTIONS 

(a)  Dr.  J.S.  Prener,  General  Electric  R&D  Laboratory,  Schenec¬ 
tady,  NY:  interested  in  crystals  of  EuMgF4  to  use  as  scintillation 
counters  in  CAT  scanners. 

(b)  Professor  R.  Reisfeld  and  Dr.  A.  Loewenschuss  of  the 
Hebrew  University,  Jerusalem:  collaboration  on  study  of  spectra  of 

+2 

Sm  in  the  compounds  explored  in  the  JSEP  work  unit  under  laser 
excitation  from  15°K  to  room  temperature. 
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F.  ELECTRONIC  PROPERTIES  OF  HYDROGEN  IN  SOLIDS 

Professor  D.C.  Mattis 

Unit  SSI-6 

1.  OBJECTIVE(S) 

To  develop  first  principles  and  theoretical  techniques  for  relating 
the  observed  behavior  of  the  hydrides  of  transition  and  rare  earth 
metals  and  their  alloys  to  their  underlying  electronic  properties.  The 
mathematics  will  be  applied  to  the  bulk,  surface  and  transport  prop¬ 
erties  of  these  materials .  A  study  of  the  broader  question  of  the 
electronic  properties  of  hydrogen  in  metals  will  be  conducted. 

2.  APPROACH 

The  behavior  of  hydrogen  in  solids  is  fast  becoming  one  of  the 
most  important  and  unsolved  problems  of  modern  science.  Despite 
astronomical  efforts  resulting  in  a  vast  body  of  mostly  empirical  knowl¬ 
edge,  little  has  been  achieved  in  the  nature  of  comprehensive  theory  or 
an  understanding  that  can  be  used  for  predictions.  Currently,  interest 
in  the  study  of  the  electronic  properties  of  hydrogenated  solids  is 
increasing,  especially  because  of  the  existence  of  major  unsolved  prob¬ 
lems  whose  resolution  bears  on  the  use  of  solids  in  electronic  energy 
conversion  and  transmission.  On  a  fundamental  level,  the  electronic 
structure  of  defects  lies  at  the  forefront  of  current  theoretical  interest. 
Substitutional  or  interstitial  hydrogen  is  one  of  the  simplest  possible 
solid-state  defects,  a  near  ideal  testing  ground  for  theories  of  impuri¬ 
ties  in  solids.  Its  simple  electronic  structure  also  accounts  for  the 
many  studies  of  its  interaction  with  surfaces,  catalysis  and  chemisorp¬ 
tion.  Recent  theoretical  work  has  laid  the  groundwork  for  attacking 
many  of  the  important  roles  of  hydrogen  in  solids  from  a  fundamental 
point  of  view. 

The  combination  of  the  approach  of  J.  Friedel  and  co-workers, 
relating  the  observed  behavior  of  transition  metals  to  many-body  ef¬ 
fects,  and  methods  developed  by  us  for  studying  the  response  of  inter¬ 
acting  particles  to  various  defects,  will  form  the  foundation  for  a  theory 
of  hydrogen  in  transition  and  rare-earth  metals  that  has  direct  applica¬ 
bility  to  the  many  and  startling  effects  observed  in  these  systems. 
Ground  state  properties  of  the  unperturbed  metal  system  will  be  des¬ 
cribed  by  a  variational  Gutwiller-type  wave  function,  and  schemes  will 
be  developed  for  applying  modern  techniques  to  the  study  of  its  exci¬ 
tation  spectrum.  The  effect  of  hydrogen  impurities  is  then  treated  by 
many-body  perturbation  theory  of  both  the  ground  state  and  its  exci¬ 
tations,  in  order  to  permit  the  description  of  time- independent  and 
time- varying  properties  of  the  solids  resulting  from  the  presence  of 
hydrogen  impurities . 
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3.  PROGRESS 

The  main  task  within  this  unit  during  the  past  year  has  been  an 
extensive  investigation  of  the  conditions  under  which  a  neutral  hydro¬ 
gen  atom  can  bind  to  a  positively  charged  particle.  This  problem  is 
also  that  of  an  exciton  trapped  at  a  positive  defect  inside  a  solid,  which 
represents  one  possible  state  of  a  proton  within  the  solid.  The  first 
analysis  has  shown  that  for  such  binding  to  occur,  the  mass  of  the 
positive  particle,  or  hole,  must  exceed  1.46  of  the  mass  of  the  free 
electron.  This  is  an  important  finding  because  it  defines  necessary 
conditions  for  the  electronic  properties  of  hydrogen  in  solids.  The 
first  analysis  used  approximate  wave  functions  for  a  three-body  system 
related  to  those  of  the  H2  ion.  However,  if  the  mass  of  the  positive 
particle  is  comparable  to  that  of  the  electron,  its  dynamic  effects  are 
also  of  the  same  order,  and  a  many -particle  wave  function  of  both  light 
particles,  and  sensitive  to  their  interaction,  is  called  for.  The  ground 
state  of  the  system  must  be  found  by  variational  techniques  that  mini¬ 
mize  the  energy  of  the  system  with  respect  to  several  unconstrained 
parameters.  The  generalization  of  the  first  analysis  to  include  these 
aspects  is  now  under  study,  using  different  forms  of  the  trial  wave 
functions .  A  preliminary  account  of  this  work  has  already  been  given . 1 
Additional  publications  that  resulted  from  earlier  work  under  this  unit 
are  also  listed  below.2*4 

Professor  Mattis  was  grante  i  a  leave  of  absence  for  part  of  last 
year;  however,  his  research  on  this  unit  continued  with  the  collabora¬ 
tion  of  a  colleague  and  students,  but  with  reduced  effort.  This  unit 
will  now  be  terminated,  however,  because  his  leave  of  absence  (at  the 
University  of  Utah)  is  continuing. 

4.  PUBLICATIONS 

1.  M.  Farid  and  D.C.  Mattis,  "Does  the  Hydrogen  Atom  Bind  a 
Proton?"  Bull.  Am.  Phys.  Soc.  26,  617  (1981). 

2.  R.E.  Pena  and  D.C.  Mattis,  "Optical  Lineshapes  and  Elec¬ 
tronic  Eigenstates  of  Substitutional  Defects  and  Vacancies,"  J. 
Phys.  C:  Solid  State  Physics  14  ,  649  (1981). 

3.  A.  Theumann,  "Ferromagnetic  and  Spin-Glass  Behavior  in  the 
Random  Hierarchical  Model,"  Phys.  Rev.  B  22,  5441  (1980). 

4.  P.S.  Riseborough  and  D.C.  Mills,  "Frequency  Dependent 
Conductivity  of  the  Anderson  Lattice,"  in  Valence  Fluctuations 
in  Solids,  ed.  L.  Falicov,  North  Holland  Publishina  Co. 
(1581). 
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A.  REDUCED-ORDER  MODELS  OF  LARGE  SCALE  SYSTEMS 

Professors  D.C.  Youla,  F.  Kozin  and  J.J.  Bongiorno,  Jr. 

Unit  IE1-1 

1.  OBJECTIVE(S) 

Our  objective  is  to  produce  reduced-order  models  for  large-scale 
systems  by  a  combination  of  time-domain,  statistical  and  aggregation-like 
methods.  The  objectives  can  be  divided  into  two  major  efforts. 

(1)  Study  and  develop  useful  computationally  stable  and  efficient 
order-reduction  techniques.  These  techniques  must  yield  good  approxi¬ 
mations  of  the  dynamical  quantities  to  be  monitored. 

(2)  Study  and  assess  the  effectiveness  of  reduced-order  model 
control  designs  when  applied  to  the  original  large  scale  system. 

2.  APPROACH 

Society  is  comprised  of  a  multitude  of  large-scale  systems,  Com¬ 
munications,  transportation,  power,  defense,  waste  disposal,  and  agri¬ 
culture  are  but  a  few  of  the  system  substructures  upon  which  our 
society  is  based.  The  heavy  demands  put  upon  these  large  systems  by 
a  society  that  continues  to  develop,  as  well  as  the  desire  to  keep  the 
systems  themselves  developing  and  operating  in  some  reasonably  effi¬ 
cient,  safe,  and  cost-effective  fashion,  of  necessity,  requires  intensive 
monitoring  of  each  of  the  systems.  However,  the  magnitude  of  these 
systems  makes  in-detail  monitoring  prohibitively  costly  even  with  cur¬ 
rent  high-speed  high-capacity  computers.  Reduced-order  models, 
therefore,  must  be  developed. 

Problems  to  be  studied  for  reduced-order  models  include:  On  what 
basis  shall  we  construct  a  reduced-order  model?  What  criteria  shall  be 
employed  to  arrive  at  a  reduced-state  model  of  the  original  system? 
What  significant  dynamical  characteristics  of  large-scale  systems  can  be 
predicted  from  lower-order  models?  Finally,  how  effective  are  design 
studies  generated  by  lower-order  models  when  applied  to  the  original 
large-scale  system?  Parameters  and  state  evolution  may  not  be  known 
exactly.  Hence,  there  is  uncertainty,  either  stochastic  or  deterministic, 
in  our  knowledge  of  the  system.  This  fact  must  be  accounted  for  in 
our  development  of  order-reducing  procedures,  and  their  related  com¬ 
putational  algorithms.  Thus,  for  example,  the  methods  must  be  stable 
in  the  sense  that  small  observation  noise  will  not  lead  to  large  errors. 

In  our  proposed  research  program,  we  shall  be  concerned  with 
three  major  problem  areas: 

(a)  Analytical  and  computational  methods  for  order  reduction. 

(b)  Applicability  of  control  design  concepts,  generated  from 
reduced-order  models  to  the  original  large-scale  system. 

(c)  Implementation  to  real  large-scale  systems. 
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3.  PROGRESS 

The  usual  controllers  for  linear  large-scale  systems  have  been 
restricted  to  one-degree-of-freedom  control  methodology.  As  a  result, 
it  has  been  difficult  to  develop  meaningful  reduced-order  models  for 
such  large-scale  systems.  During  this  past  year,  we  have,  by  general¬ 
izing  the  permissible  topology,  developed  a  rigorous  theory  of  two- 
degree-of-freedom  controller  design.  We  now  present  the  background 
of  the  problem,  derive  some  of  the  theoretical  results,  and  present  by 
example  some  of  the  exciting  possibilities  of  this  new  theory. 

A.  Introduction  and  Background 

By  and  large,  the  design  of  optimal  Wiener-Hopf  controllers  for 
linear  large-scale  time-invariant  plants  has  been  restricted  to  one- 
degree-of-freedom  control  laws.  In  turn,  this  has  practically  precluded 
the  use  of  reduced-order  modeling  and  has  also  made  it  impossible  to 
incorporate  other  "hard"  design  specifications  imposed  by  sensitivity 
constraints  and  phase-gain  stability  margin  requirements,  into  a  bona- 
fide  optimality  theory. 

In  general,  every  control  problem  possesses  two  natural  degrees  of 
freedom:  one  corresponds  to  the  availability  of  the  exogeneous  set- 

point  input  u  and  the  other  to  the  availability  of  the  sensor  output  w. 
However ,  in~  the  conventional  single-loop  configuration ,  the  controller 
generates  the  plant  input  r  by  processing  only  the  difference  u-w,  and 
one  degree  of  freedom  is  completely  lost.  To  exploit  the  two  "inherent 
degrees  of  freedom  fully,  the  controller  must  process  u  and  w  indepen¬ 
dently  . 

Now,  apart  from  all  other  considerations,  the  design  of  a  linear 
controller  has  two  clearly  defined  objectives:  one,  the  transient  acqui¬ 
sition  of  a  prescribed  set-point  input  u(s)  by  the  plant  output  y(s) 
must  be  accomplished  satisfactorily  in  the^  presence  of  measurement  noise 
m(s)  and  a  saturation  constraint  on  the  plant  input  r(s);  two,  regula¬ 
tion  of  the  face  of  load  disturbance  d(s)  which  enters  through  the  plant 

must  be  acceptable In  the  most  general  situation  both  requirements 
are  present  simultaneously,  and  the  problem  is  one  of  servo-regulation. 
Available  to  the  designer  are  the  exogenous  input  u(s)  and  the  feed¬ 
back  sensor  output  w(s)  so  that  the  overall  configuration  may  be  de¬ 
picted  schematically  as  shown  in  Figure  1. 

To  take  advantage  of  both  degrees  of  freedom,  the  controller  can 
generate  an  output 

£  =  cuJi  -  Cw£  CD 

in  which  the  two  (open-loop)  transfer  matrices  Cu(s)  and  Cw(s)  may  be 

shaped  arbitrarily,  subject  solely  to  the  demand  of  closed-system  sta_ 
bility.  On  the  other  hand,  because  of  its  restrictive  topology,  the 


(1)  All  matrix  functions  of  the  complex  variable  s  =  a  +  jw  are  real  and 
rational . 
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d(s) 


m(a) 


Fig.  1:  General  multivariable  servo-regulator  configuration 

classical  single-loop  structure  of  Fig.  2  is  able  to  process  only  the 
difference  u(s)-w(s)  and  one  degree  of  freedom  is  completely  lost.  It 
is  therefore  not  ^surprising  to  find  that  under  this  limitation  the  two 
design  objectives  often  conflict  and  appear  mutually  contradictory. 
Even  worse,  usually  little  or  no  latitude  remains  for  the  incorporation  of 

sensitivity  constraints  and  phase-gain  stability  margins. ^ 


d(s) 


m(s) 


Fig.  2:  The  classical  single-loop  servo- regulator  configuration 

A  primary  goal  of  this  year's  research  has  been  the  development  of 
a  rigorous  theory  of  two-degree-of-freedom  optimal  controller  design 
which  subsumes  these  additional  design  requirements  into  an  overall 
analytic  framework.  Preliminary  investigations  indicate  that  such  a 
theory  may  actually  be  within  reach  and  a  simple  example  presented  in 
Section  C  reveals  some  of  the  exciting  possibilities.  Section  B  contains 
the  necessary  background  material. 


(2)  The  text  by  I.  Horowitz1  is  notable  for  its  appreciation  of  many  of 
these  factors. 
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B.  The  Pure  Two-Degree-of-Freedom  Servo  Problem 

In  Fig.  1  the  transfer  matrices  P(s)  of  the  plant  and  Ft(s)  of  the 
feedback  sensor  are  prescribed  and  1 


Tc(s>  =  [Cu(s)  -  Cw(s)]  (2) 

is  the  desired  controller  transfer  matrix  from  inputs  u(s),  w(s)  to 
r(s).  By  straightforward  analysis  we  obtain 


r  =  Ruu  -  Rw(Ftd  +  m)  ,  (3) 

where 

Ru  *  O  *  c„Ftp>'lcu  <4> 

and 

Rw  =  a  *  cwFtP)-lc„  .  (5) 

Evidently,  Ru(s)  and  Rw(s)  are  closed-system  transmission  matrices 
from  exogenous  inputs  u(s)  to  r(s)  and  -m(s)  to  r(s),  respectively. 


Let  Gfs)  denote  the  spectral  density  of  the  error 

C 

e(s)  =  u(s)  -  y(s)  (6) 


and  Gr(s)  that  of  r(s).  Then,  in  the  standard  quadratic-cost  setting,2 
transient  error  and  saturation  are  accounted  for  by  the  functional 
E  =  Et  +  kEs,  where(3) 


E  =  — — 
Lt  2nj 


i* 

/  Tr  G  ds 

-jo. 


(7) 


is  the  cost  incurred  by  transient  acquisition  of  u(s)  and 


E 


s 


/  Tr  (Q  G  )  ds 

-jo. 


(8) 


is  the  penalty  for  excessive  excursions  in  plant  input  r(s).  The  ad¬ 
justable  nonnegative  constant  k  permits  tradeoff  between  Et  and  Eg, 

and  Q(si  is  a  suitably  chosen  nonnegative-definite  parahermitean  weight- 

( 4  ) 

matrix.  Let  the  respective  spectral  densities  of  u(s),  d(s)  and  m(s) 

be  denoted  by  G  (sj,  G.isj  and  G  (si. 

_  u  u  m 

(3)  Tr  A  =  trace  A. 

(4)  Q(s)  =  Q#(s)  --  Q  (-si  and  W<jw)  -  0  for  all  real  i u.  (The  prime 
denotes  matrix  transpose.) 
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For  the  pure  servo  problem  with  negligible  sensor  noise,  G^fs) 
and  Gm(s)  are  both  zero  and  the  total  cost  is  given  by 

E  =  jjr  TrKi  -  PRU)  Gu(l  -  PRU)*  +  kQRu  Gy  R^]  ds  .  (9) 

J  —  JOO 

Thus,  we  are  led  to  the  following  definition. 

Definition  1:  Ry(s)  1S  said  t0  be  acceptable  for  (P,Ft)  if  there 

exists  a  controller  which  realizes  Ru(s)  as  the  designated  closed-system 

transmission  matrix  in  an  asymptotically  stable  configuration  of  generic 
type  shown  in  Figure  1. 

Clearly  then,  in  the  minimization  of  E,  Eq.  (9),  only  acceptable 
Ru's  qualify  as  admissible  candidates  and  their  parametrization  is  an 

important  first  step.  Note  also  that  our  definition  of  acceptability  does 
not  restrict  the  controller  to  be  dynamical,  so  that  Tc(s)  need  not  be 

proper.  Some  reasons  for  lifting  this  constraint  have  already  been 
offered  in  reference  2,  but,  as  a  rule,  optimal  designs  usually  exhibit 
this  property  automatically.  In  any  case,  the  matter  can  always  be 
investigated  after  the  parametrization  of  the  larger  class  of  acceptable 


Assumption  1:  The  plant  and  feedback  sensor  are  free  of  unstable 
hidden  modes,  Ft(s)  is  analytic  in  Re  s  >  0,  and  the  finite  pole  of  P(s) 

and  Ft(s)P(s)  in  Re  s  >  0  coincide  and  have  identical  McMillan  degrees. 

Except  for  the  analyticity  of  Ft(s)  in  Re  s  >  0,  all  parts  of  as¬ 
sumption  1  are  necessary  for  the  existence  of  a  stabilizing  controller.2 
Let 

A-1B  =  FtP  =  B1Al1  (10) 

where  (A,B)  is  any  left-  and  (Blf  At)  any  right-coprime  pair  of  poly¬ 
nomial  matrices.  There  exist1'2  polynomial  matrices  X(s),  Y(s),  Xx(s), 
Yi(s)  such  that 


AX  +  BY  =  1  , 

XlA1  +  Y  lB1  =  1 

(ID 

det  X  •  det  Xt  t  0 

• 

(12) 

Theorem  1  (New): 

Under  assumption  1 , 

Ru(s)  is  acceptable  for 

(P.F^)  iff  Ru(s)  and  P(s)Ru(s)  are  both  analytic  in  Re  s  >  0.  Equi¬ 
valently,  Ru(s)  is  acceptable  iff  it  is  of  the  form 

Ru(s)  =  Ai(s)  Hu(s)  (13) 

where  Hu(s)  is  any  matrix  analytic  in  Re  s  >  0. 
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Hence,  to  minimize  E  we  replace  Ru(s)  in  (9)  by  (13)  and  vary 
over  the  free  parameter  Hu(s)  analytic  in  Re  s  >  0;  optimization  there¬ 
fore  fixes  Ru(s)  but  not  Cu(s)  and  Cw(s)  individually,  as  is  seen  from 

Eq.  4.  This  latitude  enables  the  designer  to  satisfy  other  requirements 
without  sacrificing  optimality.  To  make  further  progress  it  is  necessary 
to  introduce  three  additional  assumptions  which  appear  absolutely  indis¬ 
pensable.2 

Assumption  2:  AGuA*,  A1:|t(P*P  +  kQ)Ax  and  Q  are  analytic  and 
nonsingular  on  the  finite  s  =  jui-axis. 


Assumption  3:  (Ft  -  1)P  is  analytic 

Assumption  4:  The  spectral  density 
fast  as  1/s2  as  s  ■*  »,  i.e.  ,v  ' 


on  the  finite  s  =  jw-axis. 
Gu(s)  vanishes  at  least  as 


Gu(s)  <  Ots'2) 


(14) 


Theorem  2  (New):  Suppose  that  assumptions  1-4  are  satisfied  and 

we  let  A(s)  and  fi  (s)  denote  the  Wiener-Hopf  spectral  solutions^  of 
the  two  equations  w 


A^#(P*P  +  kQ)Ax  =  A^A  , 

Gu  =  Qufiux  • 

Let 

f  =  A*  A^P*^  . 

(7) 

Then,  the  unique  acceptable  design  that  minimizes  E  is  given  byv  ' 


(15) 

(16) 

(17) 


R 


u 


A1A'1{r}+n’1 


(18) 


Comment:  The  compactness  of  (18)  is  somewhat  deceptive  because 
it  depends  both  explicitly  and  implicitly  on  the  coprime  polynomials 
Ax(s)  and  Bx(s).  It  is  obviously  desirable  to  eliminate  these  factors 
from  all  final  formulas  whenever  possible. 


(5)  A(s)  <  O(s^)  means  that  no  entry  in  A(s)  grows  faster  than  s^  as 
s  ■»  ». 

(6)  A(s)  and  flu(s)  are  square  real  rational  matrices  analytic  together 

with  their  inverses  in  Re  s  >  0;  they  are  uniquely  determined  up 
-to  constant  orthogonal  matrix  multipliers.3 

(7)  In  the  partial  fraction  expansion  of  r(s)  the  contribution  made  by 
all  the  finite  poles  in  Re  s  <  0  is  denoted  by  {T}+.  Clearly  {T}+  is 
analytic  in  Re  s  >  0  and  <  ff(l/s). 
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Corollary  1:  Let  P(s)  be  free  of  finite  poles  in  Re  s  >  0  and  let 
Ar(s)  denote  the  Wiener-Hopf  spectral  solution  of  the  "reduced"  equa¬ 
tion 


P*P  +  kQ  =  Ar*Ar  .  (19) 

Then,  under  the  conditions  of  Theorem  2,  the  optimal  design  is  given 
by 

Ru  =  'V1f'Yip«VAl  <20) 

and  all  coprime  factorizations  are  by-passed. 


This  corollary  to  Theorem  2  has  a  very  wide  range  of  practical 
applications  because  the  number  and  location  of  the  jui-axis  poles  of 
P(s)  are  totally  unrestricted.  These  "infinite  gain"  frequencies  play  an 
important  role.  For  example,  good  scalar  servo  design  utilizes  a  pole  of 
P(s)  at  s  =  0  to  produce  a  system  that  can  track  a  shape-deterministic 
step-input  with  zero  steady-state  error.  More  generally,  in  the  multi- 
variable  case,  a  design  preserves  the  steady-state  acquisition  potential 
of  the  plant  if  every  finite  ju>-axis  pole  of  P(s)  is  a  zero  of  the  deter¬ 
minant  of  the  error  transmission  matrix  1  -  P(s)R  (s)  of  sufficiently 
high  multiplicity.2 

Corollary  2:  Any  acceptable  design  Ru(s)  whose  associated  cost  E 

is  finite  preserves  the  steady- state  acquisition  potential  of  the  plant. 
A  fortiori,  so  does  the  optimal  design  (18). 

C.  An  Example 

A  single-input-output  plant  has  the  transfer  function 


p<s>  *  i?&) 


(21) 


Under  the  conditions  k  =  Q(s)  =  Ft(s)  =  1,  use  Theorem  2  to  design  a 

servo  that  tracks  a  step-input  u(s)  with  spectral  density  G  (s) 

=  -1/s2,  optimally.  (Sensor  noise  is  negligible.) 


Solution :  Assumptions  1-4  are  satisfied 
A  =  Ax  =  s(s-2)  ,  B  =  Bx  =  s  -  1 

A=s2+V7s+1 


fiu  =  i 


Thus  (details  omitted),  (18)  yields  the  optimal  design 


H. ,  = 


-1 


R..  =- 


-s(s-2) 


u  s2+V7s+l  '  u  s2+/7s+l 


3.65  . 


(22) 

(23) 

(24) 
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In  any  stable  closed-system  realization  of  R  (s),  the  tracking 
error  e  =  u-y  =  (1  -  PRu)u  where 


=  s(s  ♦  (1  ♦  V7)) 
u  s2  +  JT  s  +  1 


(25) 


Hence,  because  of  the  presence  of  the  factor  s  in  the  numerator  of 
(25),  any  such  servo  can  track  a  step-input  with  zero  steady-state 
error,  as  predicted  by  corollary  2. 

The  single- loop  realization  of  R  (s)  shown  in  Fig.  3  cannot  work, 

since 


Ru  =  (1  +  CP)"1C  , 

C  =  R  (1  -  PR,,)'1  = - 2-~-S-  - 

u  u  s  ♦  1  +  p 


(26) 

(27) 


and  stability  is  precluded  by  the  Re  s  >  0  pole-zero  cancellation  be¬ 
tween  P(s)  and  C(s)  at  s  =  2. 


I  I 


j  Contrnlu*r  I 

I - j 


Fig.  3.  Standard  single-loop  realization  of  the  optimal  design  Ru(s). 

However,  if  we  insist  on  a  single-loop  topology,  the  smallest  pos¬ 
sible  realizable  cost  is  E  =  131  which  is  achieved2  by  the  design 

R  -S(s-2)((11  t  4V7)s-2)  c  =  (11  +  4/7  s)-2  (2i 

u  (s  +  2)(s2  +  /Ts  +  2)  '  s-(7  +  3/7) 

To  accept  this  design  is  to  accept  an  unnecessary  36- fold  deqrada- 

ZO\  ' 

tion.in  performance! K  J  Actually,  the  cost  E  =  3.65  is  easily  attained 
by  means  of  minor-loop  rate-feedback  around  the  plant  (Fig.  4). 

(8)  The  degradation  is  even  more  dramatic  for  k  >  1. 
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I _ I 


C  =  C  ;  C  +  F  =  C 
m  u'  m  m  w 

Fig.  4.  Realization  of  the  cost  E  =  3.65  with  minor-loop  feedback. 


To  ensure  closed-system  stability,  the  gain  constant  a  must  be  re¬ 
stricted  to  lie  in  the  interval 


-2  <  a  <  -  1  ,  (29) 

and  the  obvious  midpoint  candidate  a  =  -1.5  corresponds  to  the  first- 
order  lag  compensator 


C 


m 


1  +  s 

2(s  +  1  +  -JT) 


(30) 


But  as  is  shown  in  Fig.  5,  one  can  also  achieve  the  cost  E  =  3.65  by 
introducing  additional  compensation  in  the  feedback  path  of  the  single¬ 
loop  configuration  of  Figure  3. 


Fig.  5.  Modified  single-loop  realization  of  the  cost  E  =  3.65. 
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Evidently,  to  decide  between  the  two  designs,  other  criteria  must  be 
invoked . 

The  results  of  a  detailed  study  of  the  sensitivity  of  the  input- 
output  transfer  function  T(s)  =  P(s)Ru(s)  may  be  summarized  briefly  as 

follows .  ^ 

1)  From  0  to  83.57  Hz,  the  sensitivity  of  T  with  respect  to 
F  in  Fig.  5  exceeds  that  of  T  with  respect  to  Fm  in 

Figure  4.  In  particular,  over  the  range  0  to  11.44  Hz, 
it  exceeds  it  by  at  least  17.3  db! 

2)  From  0  to  42.73  Hz,  the  sensitivity  of  T  with  respect  to 
P  (to  C)  in  Fig.  5  exceeds  that  of  T  with  respect  to  P 
(to  Cm)  in  Figure  4. 

Since  the  3-db  bandwidth  of  T(s)  is  less  than  6.28  Hz,  it  is 
concluded  that  the  minor-loop  design  of  Fig.  4  offers  better  sensitivity 

performance  with  respect  to  all  parameters  of  interest.  Neverthe¬ 

less,  we  are  not  suggesting  that  the  minor-loop  structure  is  invariably 
superior:  On~tRe  contrary,  the  subject  is  wide  open,  especially  in  the 
multivariable  case! 


Over  the  bandwidth  of  the  plant,  the  ratio  of  the  magnitude  of  (31)  to 

that  of  (32)  reaches  a  maximum  of  =  6,  but  for  large  tu  it  grows  without 

bound.  Apparently  then,  the  different  sensitivity  behavior  exhibited 

by  the  two  designs  is  intimately  related  to  the  different  choice  of  closed- 

system  transmission  function  R  (s). 

_  w 

(9)  a  =  -1.5. 

(10)  In  our  opinion,  the  improperness  of  the  minor-loop  compensator  is 
of  no  great  practical  consequence. 

(11)  In  Fig.  4,  C  =  Cm  and  C  +  Fm  =  C  but  in  Fig.  5,  C  =  C„  and 

ill  III  III  W  U 

‘  CF  =  Cw-  Thus,  (31)  and  (32)  are  obtained  by  using  the  formula 
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5.  DoD  AND  OTHER  INTERACTIONS 

(a)  Under  an  arrangement  that  the  Rome  Air  Development  Center, 
New  York,  has  with  Syracuse  University,  New  York,  known  as  SCEE, 
Professor  Youla  has  performed  consulting  work  for  Dr.  Heywood  Webb 
of  RADC  each  summer  for  the  past  three  summers.  The  topics  which 
he  first  explored  under  JSEP  sponsorship,  were:  Engineering  applica¬ 
tions  of  the  Prony  technique;  Recovery  of  spectral  information  from 
analogue  covariance  data;  and  Diagnostics  for  nonuniform  media. 

(b)  On  an  extension  of  work  originally  supported  by  JSEP, 
Professor  Youla  has  recently  had  a  contract  with  RADC  on  the  topic  of 
image  restoration. 

(c)  Professor  Youla  was  invited  to  present  a  talk  on  his  new 
theory  of  spectral  estimation,  developed  under  this  work  unit,  at  the 
meeting  on  "Techniques  of  Spectral  Estimation,"  held  at  the  Naval 
Postgraduate  School,  Monterey,  California,  during  1980,  and  sponsored 
by  the  Office  of  Naval  Research.  Many  Navy  researchers  were  present, 
and  very  good  interaction  resulted  with  Professor  Youla  concerning 
their  special  problems. 

(d)  Professor  Youla  was  invited  to  present  talks  on  the  subject  of 
spectral  estimation  at  Texas  Tech  and  at  the  University  of  Maryland. 
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B.  ENHANCEMENT,  EXTRACTION  AND  RECONSTRUCTION  TECH- 

NIQUES  IN  IMAGE  PROCESSING 

Professors  A.  Papoulis  and  L.  Kurz 

Unit  IE1-2 

1.  OBJECTIVE(S) 

The  main  purpose  of  this  research  effort  is  to  investigate  some 
new  and  promising  approaches  to  the  general  problem  of  image  enhance¬ 
ment,  extraction,  and  reconstruction.  The  investigation  will  encompass 
two  philosophically  different  classes  of  techniques:  spectral-estimation 
and  statistical. 

In  connection  with  the  spectral-estimation  approach,  it  is  proposed 
to  develop  a  new  method  of  filtering,  combining  frequency  domain 
techniques  with  the  advantages  of  adaptive  and  recursive  estimators. 
The  method  is  based  on  the  representation  of  an  arbitrary  signal  f(t)  in 
terms  of  the  samples  of  its  "running  FFT,"  and  it  can  be  used  without 
any  knowledge  of  prior  statistics.  A  special  form  of  the  filter  leads  to 
an  extension  of  the  gradient-seeking  algorithm  used  in  a  variety  of  LMS 
estimators.  The  study  will  include  the  development  of  the  properties  of 
running  transforms  and  their  use  in  spectral  estimation. 

The  statistical  class  of  techniques  will  include  the  study  of  several 
types  of  statistical  maks,  edge  detection  procedures,  and  image  recon¬ 
struction  algorithms.  In  addition,  feature  extraction  and  object  identi¬ 
fication  based  on  various  versions  of  factor  analysis  techniques  will  be 
studied . 

The  prime  objective  of  the  research  effort  is  to  obtain  general 
results  with  broad  implications  in  such  areas  as  interpretation  and 
classification  of  pictorial  data,  processing  of  radar  data  of  all  forms, 
geological  and  underwater  sounding,  meterological  and  oceanographic 
data,  space  and  communications  data,  etc.  Whenever  feasible,  compara¬ 
tive  studies  of  the  two  classes  of  approaches  for  a  given  application  will 
be  introduced. 

Among  others,  the  following  aspects  will  be  considered:  various 
features  of  image  processing  with  emphasis  on  noise  and  clutter  reduc¬ 
tion;  detection  of  quasimonochroma  tic  sonar  signals  in  the  presence  of 
echoes  and  engine  interference;  rapid  measurement  of  the  instantaneous 
frequency  of  FM  signals;  adaptive  equalizers;  recursive  estimation  of 
covariance  matrices. 

2.  APPROACH 

The  two  different  basic  approaches  are  discussed  separately  in  the 
presentation  below. 
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A.  Spectral-Estimation  Approach 

To  place  the  proposed  spectral  method  into  familiar  perspective,  we 
shall  relate  one  of  its  applications  to  adaptive  LMS  estimation.  We  shall 
use  an  illustration  the  estimation  of  a  discrete  process  y[n]  in  terms  of 
the  output  (see  Fig.  1). 


N-l  T 

?[n]  =  I  a.  [n]  x[n-k]  =  Xi[n]  A[n]  (1) 

k=0  K 

of  a  non-recursive  time-varying  filter  of  length  N  with  the  data  process 
X[n]  as  the  input.  Our  objective  is  to  determine  the  weights  a.[n]  so 
as  to  minimize  the  MS  value  of  the  estimation  error. 

e[n]  =  y[n]  -  y[n]  (2) 

If  the  statistics  of  the  processes  x[n]  and  y[n]  are  known,  then 
(Wiener  filter)3 

A[n]  =  R_1r  (3) 

where 

R  =  E{X[n]XT[n]}  T  =  E{X[n]  y[n]}  (4) 

In  the  absence  of  this  information,  we  can  either  estimate  the 
covariance  matrix  R  and  the  cross-covariance  vector  r,  or  we  can 
determine  A[n]  adaptively.  A  simple  method  of  adaptation  based  on  an 
instantaneous  version  of  a  gradient- seeking  technique  for  minimizing  the 
MS  error  is  the  algorithm 

ak[n+l]  =  akfn]  +  p  e[n]  x[n-k]  (5) 

introduced  by  Widrow1'4  and  used  in  several  areas  of  signal  process¬ 
ing.5 


The  Widrow  filter  is  simple;  it  requires  no  prior  knowledge  of  any 
statistics,  and  it  can  be  applied  to  stationary  and  non-stationary  pro¬ 
cesses.  It  has,  however,  a  number  of  disadvantages:  The  number  N 
of  coefficients  that  are  adaptively  controlled  equals  the  length  N  of  the 
filter.  This  often  introduces  unnecessary  complications  if  N  is  large. 
Filtering  of  noise  is  best  performed  in  the  frequency  domain.  This 
leads  to  a  small  number  of  unknown  parameters  if  they  control  directly 
the  frequency  response  of  the  filter.  The  filter  of  Fig.  1  does  not 
have  that  property.  The  solution  of  (5)  is  a  random  set  of  weights 
A[n]  that  does  not  approach  the  Wiener  optimum6'7  [see  (3)].  The 
estimation  can  be  simplified  if  the  adaptively  controlled  coefficients  are 
associated  with  nearly  uncorrelated  data  (Karhunen-Lofeve).3  This  is 
not  the  case,  in  general,  with  the  data  x[n-k]  in  (1).  The  length  N  of 
the  filter  is  fixes;  this  is  often  undesirable. 
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*[n]  x[n-Q  xQpi-2]  x[n-N+l] 


Fig.  2.  Frequency  domain  adaptive  estimator. 
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To  overcome  these  difficulties,  we  impose  the  restriction  that  the  N 
weights  A^[n]  are  linearly  dependent  on  MQ  parameters  bm [n] : 

A[n]  =  WB[n] 

In  the  above,  B[n]  is  a  vector  consisting  of  M  elements  and  W  is  an 
NXMq  matrix  to  be  determined 

y[n]  =  XT[n]  W  B[n]  =  ZT[n]  B[n]  (6) 


where 


Z[n]  =  WTX[n]  (7) 

Equations  (6)  and  (7)  specify  the  system  of  Fig.  2,  consisting  of 
the  matrix  W  transforming  the  N-vector  X[n]  into  the  MQ-vector  Z[n], 

followed  by  the  Mq  weights  B(n]  and  an  adder.  The  vector  B[n]  is 

adaptively  controlled: 

B [n+1]  =  B[n]  +  p  e[n]  Z[n]  (8) 

as  in  (5).  However,  unlike  (5),  we  have  only  Mq  unknown  weights 
and,  as  we  shall  see,  Mq  can  be  much  smaller  than  N. 

The  selection  of  the  matrix  W  depends  on  the  properties  of  the 
processes  x[n]  and  y[n].  Suppose,  for  example,  that 

x[n]  =  y[n]  +  u[n]  (9) 

where  u[n]  is  noise.  The  vectors  X[n]  and  Y[n]  and  the  elements 
x[n-k]  and  y[n-k],  respectively  are  points  in  an  N -dimensional  space. 

If  the  signal  vector  Y[n]  belongs  to  a  subspace  of  dimensionality 

Mq,  then  we  choose  for  W  the  projection  operator  into  S^. 

In  the  selection  of  W,  we  are  guided  by  the  following:  small  Mq; 

data  vector  Z[n]  with  uncorrelated  components;  frequency  domain 
properties  of  Z[n].  To  meet  these  requirements,  we  select  for  W  a 
matrix  whose  elements  are  the  roots  of  unity: 

W1™,  W  =  e]2n/N,  0  <  k  <  N-l ,  -M  <  m  <  M  (10) 

where  MQ  =  2M+1  and  the  index  is  suitably  changed.  The  elements  of 
the  resulting  vector  Z[n]  are  given  by  [see  (7)] 

N-l  -if 

F[n,m]  =  I  x[n-k]  WmK  (11) 

k=0 


that  is,  they  are  the  discrete  Fourier  series  (DSF)  of  the  N  samples 
x[n] , . . .  ,x[n-N+l]  of  the  input  x[n]. 
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With  this  choice  of  W,  the  output  of  the  filter  of  Fig.  2  is  given 
by  [see  (6)] 

M 

y[n]  =  I  b  [n]  F[n,m]  (12) 

m=-M  m 

and  the  coefficients  bm[n]  are  adaptively  controlled: 

bm[n+l]  =  bm[n]  +  e[n]  F[n,m]  (13) 

B.  Various  Statistical  Approaches 

The  masking  techniques  and  techniques  based  in  ANOVA  (analysis 
of  variance)  used  by  the  senior  investigator  (Professor  Kurz)  in  the 
past  may  be  unified  within  the  framework  of  the  statistical  linear 
models.  In  a  recent  paper,8  he  introduced  a  collection  of  simple 
procedures  based  on  masks  to  be  used  in  computer  image  processing. 
Small,  carefully  selected  regions  of  an  image  are  viewed  at  one  time. 
The  process  is  repeated  throughout  the  image.  Many  different  types  of 
masks  were  introduced  which  accomplish  efficiently  the  tasks  of  edge 
detection,  noise  reduction,  shrinking,  expansion,  object  recognition, 
measuration,  etc.  In  a  separate  paper9  the  effectiveness  of  ANOVA 
techniques  were  demonstrated  in  image  processing  and  enhancement. 
Simultaneously,  classes  of  Robbins-Munro  type  stochastic  approximation 
minimum  variance  least  squares  (SAMVLS)  estimators  were  devel¬ 
oped10'11'12  for  which  the  resultant  regression  function  is  interpreted 
as  a  recursive  correlator.  Robustness  (relative  insenitivity)  to  varia¬ 
tions  distribution  and  insensitivity  to  contamination  of  data  is  achieved 
by  batch  preprocessing  of  the  data  either  by  using  linear  rank  tests 
(batch-nonlinear-linear  or  BNL  approach)  or  by  using  adaptive  gain 
coefficients  with  properly  preselected  nonlinearity.  By  combining  the 
theory  of  ANOVA,  experimental  design  and  robustized  SAMVLS,  new 
three-dimensional  recursive  masks  may  be  introduced  for  real-time 
parallel  processing  and  detection  of  stationary  or  moving  images  in 
unspecified  noise  environment  from  observations  taken  of  a  noisy  scene. 
The  method  has  the  potential  for  development  of  techniques  for  esti¬ 
mation  of  velocity  and  acceleration  of  moving  objects.  In  addition,  the 
theory  of  Latin  Square  experimental  designs  is  particularly  attractive  as 
an  approach  to  the  design  of  three-dimensional  masks. 

In  a  separate  paper,13  it  was  demonstrated  that  two-dimensional 
partition  tests  are  effective  in  achieving  excellent  quality  edge  detection 
in  the  presence  of  severe  and  undeformed  noise.  The  philosophy  of 
this  approach  needs  further  extension  to  include  nonlinear  processing 
which  is  particularly  useful  in  the  processing  of  images  appearing  in 
radar  and  underwater  sounding  data.  The  extended  model  will  include 
the  problem  of  moving  objects,  velocity,  acceleration,  etc. 

Using  masking  techniques  in  conjunction  with  occupancy  vectors  of 
partition  tests  represents  an  interesting  marriage  of  the  two  approaches 
outlined  above.  The  full  impact  of  this  hybrid  version  of  image  pro¬ 
cessing  problems  is  not  clear  at  this  time  but  preliminary  results  indi- 
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cate  its  potential  for  efficiently  solving  numerous  problems  in  image 
processing. 

The  SAMVLS  algorithm  mentioned  above12  permits  one  also  to 
attack  the  problem  of  image  reconstruction.  Preliminary  results  in  this 
area  are  available  in  reference  14.  The  preliminary  results  of  this 
approach  are  promising  and  justify  a  major  effort  in  this  area. 

Two  classes  of  recursive  factor  analysis  techniques  as  a  means  of 
achieving  feature  extraction  and  image  identification  were  introduced  in 
reference  15.  This  approach  is  particularly  useful  in  processing  mass¬ 
ive  arrays  of  data.  Further  development  and  refinement  of  this  ap¬ 
proach  will  result  in  the  attainment  of  an  efficient  software  library  for 
reducing  and  identifying  severely  corrupted  and/or  unknown  images 
imbedded  in  an  unknown  background  noise. 

3.  PROGRESS 

Progress  will  be  reported  separately  for  the  spectral  estimation  and 
statistical  approaches. 

A.  Spectral  Estimation  Approach 

During  the  past  year,  we  have  treated  the  following  problems. 

(1)  A  new  approach  to  the  bandlimited  extrapolation  algorithm 

As  we  noted  in  last  year's  progress  report,  the  extrapolation 
algorithm  introduced  earlier16,17  for  the  reconstruction  of  finite  objects 
in  terms  of  their  diffraction  limited  images,  is  effective  primarily  for 
point  sources18'19  For  extended  objects,  the  convergence  is  slow  and 
the  results  are  strongly  affected  by  noise.  To  eliminate  these  prob¬ 
lems,  we  have  introduced  a  modification  of  the  algorithm  based  on  a  new 
iteration  whose  limit  is  no  longer  bandlimited  but  whose  transform  is 
close  to  the  aperture  field  of  the  unknown  object.  By  a  suitable  selec¬ 
tion  of  the  iteration  parameters,  we  achieve  a  balance  between  accuracy 
of  image  restoration  and  noise  effects.  We  are  in  the  process  of  testing 
the  method  numerically. 

(2)  Application  of  running  FFT’s 

The  method  of  running  FFT’s  developed  earlier20  has  been  applied 
to  a  variety  of  problems  involving  filtering  and  image  restoration. 
Recent  results  involve  the  application  of  the  method  to  the  problem  of 
detecting  a  target  in  the  presence  of  strong  clutter.  We  have  devel¬ 
oped  a  time-varying  filter  wrose  parameters  are  adaptively  controlled 
according  to  the  local  statistics  of  the  environment.  The  processing  is 
carried  out  in  the  frequency  domain  and  the  cut-off  frequencies  of  the 
filter  are  estimated  with  various  threshold  techniques.  Unlike  other 
methods,  this  approach  does  not  assume  prior  knowledge  of  background 
statistics.  The  method  is  particularly  effective  for  the  clutter  reduction 
problem  because  the  filtering  is  carried  out  directly  in  the  frequency 
domain . 
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To  improve  the  effectiveness  of  the  method,  we  are  currently 
considering  prediction  techniques  based  on  focal  plane  processing.21 

(3)  Maximum  entropy 

We  have  considered  various  aspects  of  spectral  estimation  based  on 
the  method  of  maximum  entropy.  Our  results  center  on  the  following 
two  areas:  (a)  A  critical  examination  of  the  method  and  its  relationship 
to  the  Wiener  theory  of  prediction,  autoregressive  processes,  the  Levin¬ 
son  algorithm,  lattice  filters,  all-pole  and  all-pass  systems,  and  stabil¬ 
ity.22'24  (b)  Reexamination  of  the  method  based  on  prior  information 

about  the  bandwidth  of  the  unknown  spectrum.  Introduction  of  win¬ 
dows23  for  reducing  the  effects  of  oversampling  and  the  presence  of 
spurious  lines  in  the  all-pole  estimate  of  the  unknown  spectrum. 

B.  Various  Statistical  Approaches 

In  connection  with  the  work  on  statistical  approaches  to  the  prob¬ 
lems  of  image  enhancement,  extraction  and  reconstruction,  several 
important  new  results  were  obtained  and  some  of  the  past  results  were 
extended  and  sharpened. 

(1)  Radar,  sonar  and  earthsounding  images  result. from  particu¬ 
larly  weak  signals.  In  addition,  the  underlying  distributions  in  the 
latter  class  of  problems  are  poorly  defined.  Faced  with  these  difficul¬ 
ties,  numerous  investigations  in  this  area  concentrated  their  research 
efforts  on  array  processors,  preferably  of  the  robust  type.26,27'28  In 
the  last  year  some  new  results  in  this  area  were  presented  in  a 
paper.25  The  typical  assumptions  of  gaussianity  and  stationarity  of 
additive  noise  were  removed.  In  addition,  it  was  assumed  that  the 
noise  at  the  sensor  inputs  are  inter-correlated.  To  robustize  the  array 
processing  procedure,  the  optimal  nonlinearities  at  the  correlator's 
inputs  were  approximated  by  simplified  nonlinear  functions  based  on 
parameters  which  can  be  easily  estimated.  The  effect  of  the  inter¬ 
sensor  noise  dependence  on  the  asymptotic  relative  efficiency  (ARE)  of 
a  multiple-input  correlation  detection  procedure  was  examined.  Easy  to 
implement  detector  realizations  were  shown  to  be  efficient  when  com¬ 
pared  with  the  optimum  detector  configurations. 

(2)  In  last  year’s  report  three  classes  of  problems  —  object 
reconstruction  algorithms,29  extraction  of  moving  objects  by  recursive 
parallelopiped  masks30  and  edge  detection  using  Latin  Square 
masks31  --  were  described  which  depended  strongly  on  a  good  robust 
recursive  vector  parameter  estimator.  This  year  a  refined  version  of 
such  an  estimator  was  generated.32  An  improved  rate  of  convergence, 
guaranteed  low  variance  of  estimation  and  high  degree  of  robustness 
(insensitivity  to  variations  in  underlying  noise  distributions)  result  from 
the  use  of  rank-type  preprocessors.  The  latter  class  of  preprocessors 
act  as  efficient  guassianizers  of  poorly  conditioned  data. 

(3)  The  statistical  mask  procedures  presented  in  references  30 
and  31  may  be  considered  as  belonging  to  a  class  of  masks  based  on 
experimental  designs.  Another  important  class  of  problems  was  consid- 
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ered  in  the  last  year  --  trajectory  detection  --  which  may  also  be 
treated  using  statistical  masks  based  on  the  theory  of  experimental 
designs.  In  particular,  masks  based  on  the  symmetrical  balanced  incom¬ 
plete  block  designs  (SBIB)  were  studied  in  detail. 

A  trajectory  is  identified  by  a  continuous  and  narrow  region  of 
high  intensities  which  is  surrounded  by  a  region  of  low  intensities  on 
each  side.  When  an  image  containing  such  a  trajectory  is  processed  by 
a  localized  operator  (window  or  mask),  the  grey  level  distribution 
within  the  mask  behaves  like  either  a  pulse  function  or  a  step  function, 
depending  on  whether  the  mask  is  placed  perpendicularly  or  in  parallel 
to  the  trajectory-  The  traditional  square-shaped  edge  detector  is 
ineffective  in  this  application  regardless  of  the  relative  position  of  the 
mask  to  the  trajectory.  A  small  edge  detector  may  not  contain  enough 
variations  due  to  trajectory  elements  to  cause  significant  block  effects. 
On  the  other  hand,  a  large  edge  detector  may  cover  too  wide  an  area 
so  that  adjacent  trajectory  segments  are  included  in  the  same  mask. 
What  is  needed  is  an  experimental  design  in  the  shape  of  a  stretched 
rectangle.  When  this  design  is  placed  perpendicularly  to  the  trajec¬ 
tory,  it  allows  for  inclusion  of  three  component  parts  of  a  pulse  func¬ 
tion.  When  the  mask  is  placed  in  parallel  to  the  trajectory,  it  allows 
for  the  inclusion  of  substantial  elements  from  only  one  trajectory  seg¬ 
ment  to  form  a  step  function. 

A  trajectory  is  further  characterized  by  two  inherent  properties: 
the  block  nature  and  strong  contrast  of  trajectory  elements.  Firstly, 
the  elements  of  the  actual  trajectory  tend  to  be  adjacent  to  one  another. 
Their  grey  level  or  texture  tend  to  be  similar  when  compared  with  the 
surrounding  elements  on  both  sides  of  the  trajectory.  Secondly,  the 
grey  level  distribution  changes  drastically  on  both  sides  of  the  trajec¬ 
tory.  Thus  only  those  block  effects  which  are  associated  with  strong 
contrast  can  be  considered  as  indicative  of  a  "true"  trajectory.  Based 
on  these  properties,  an  image  data  can  be  considered  as  a  three- 
dimensional  array  with  the  row  and  column  number  describing  the 
planar  position  and  the  grey  level  distribution  being  described  by  the 
third  dimension.  In  terms  of  experimental  designs,  the  analysis  of  an 
image  containing  a  trajectory  calls  for  a  model  which  can  determine  the 
presence  of  row,  column  or  grey  level  effects  independently.  An  SBIB 
design  satisfies  all  the  required  conditions:  it  is  shaped  as  a  stretched 
rectangle,  its  analysis  determines  if  there  exists  significant  block  or 
grey  level  effects  independently,  and  the  adjusted  block  means  can  be 
used  as  unbiased  estimators  of  the  true  block  means.  In  its  application 
to  the  trajectory  detection  problem,  the  hypothesis  that  a  SBIB  mask 
contains  a  trajectory  portion  is  accepted  if  the  following  three  condi¬ 
tions  are  satisfied  simultaneously:  (a)  there  exist  significant  block 
effects;  (b)  there  exist  significant  grey  level  effects;  (c)  the  adjusted 
block  means  behave  like  a  pulse  function.  To  satisfy  the  randomization 
criterion,  a  SBIB  design  is  selected  from  a  set  of  existing  designs 
randomly  and  then  its  grey  level  designations  are  further  randomized 
by  permutations  of  the  row  and  column  assignments.  The  first  two 
conditions  can  be  easily  implemented  by  comparing  the  F  statistics  of  a 
SBIB  design  to  a  threshold  value  corresponding  to  a  preselected  con¬ 
fidence  level.  The  shape  of  the  grey  level  distribution  can  be  de- 
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termined  by  thresholding  the  adjusted  block  means  by  their  overall 
means . 

The  ideas  outlined  above  have  been  implemented  and  found  to 
produce  excellent  results  even  if  the  trajectories  are  severely  corrupted 
by  noise.  A  research  paper  based  on  these  ideas  is  in  preparation  and 
should  be  submitted  for  publication  in  the  fall  of  1981. 

(4)  Detecting  edges  separating  grey  level  is  an  important  aspect 
of  image  processing.  A  new  procedure  for  detecting  edges  based  on 
radial  F -statistics  has  been  developed  in  the  last  year.  Though  the 
procedure  is  uniformly  most  powerful  (UMP)  for  a  gaussian  noise  back¬ 
ground,33  it  remains  robust  (insensitive  to  changes  in  noise  distribu¬ 
tion)  in  many  situations.34  If  conditions  of  reference  34  are  not  satis¬ 
fied,  robustness  of  the  procedure  may  be  achieved  using  robustizing 
algorithms  developed  in  references  29  and  32.  The  actual  detection 
method  is  based  on  the  use  of  F-statistics  along  radial  vectors  from 
inside  of  the  object  enclosed  by  the  edge.  To  find  the  location  of  a 
region  inside  an  object,  several  simple  procedures  to  accomplish  this 
task  have  been  generated.  It  was  shown  that  the  radial  F-statistic 
along  a  vector  beginning  inside  an  object  varies  linearly  between  two 
constant  values  in  passing  between  two  grey  levels.  The  midpoint  of 
the  linear  portion  represents  a  true  edge  point.  The  pertinent 
F-statistic  is  of  the  form 
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where  V?  =  sample  variance  of  the  extracted  sample,  V3  sample  vari¬ 
ance  of  xhe  training  sample  and  the  notation  of  reference  35  has  been 
used.  The  ideas  outlined  above  have  been  thoroughly  tested  on  various 
corrupted  images  yielding  excellent  results.  The  final  version  of  the 
paper  based  on  these  ideas  is  in  preparation. 
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5.  DoD  AND  OTHER  INTERACTIONS 

With  respect  to  the  contributions  of  Professor  Papoulis,  we  should 
like  to  single  out  three  contributions  which  resulted  from  JSEP- 
supported  investigations  that  have  recently  attracted  general  interest, 
and  for  two  of  which  personalized  names  were  attached. 

(a)  The  Papoulis-Gerchberg  method  of  spectral  estimation 

This  is  an  iterative  algorithm  that  yields  the  power  spectrum  of  a 
single  x(t)  in  terms  of  only  a  segment  of  x(t).  The  method  is  based 
on  the  assumption  that  it  reconstructs  the  unknown  portion  of  the 
signal  by  utilizing  the  available  information  in  the  time  and  frequency 
axis.  The  method  is  particularly  effective  if  the  unknown  spectrum 
consists  of  lines  (hidden  periodicities). 

(b)  The  Papoulis  window 

This  is  an  optimum  window  developed  for  high  resolution  spectral 
estimates  based  on  the  earlier  Blackman-Tukey  techniques. 

(c)  The  generalized  sampling  expansion 

We  developed  a  method  for  reconstructing  a  bandlimited  signal  x(t) 
in  terms  of  the  samples  of  N  functionals  y.  (t)  of  x(t)  sampled  at  1/N 
the  Nyquist  rate.  Recently,  this  result  was  applied  by  several  authors 
to  a  variety  of  problems  in  signal  processing. 

The  remaining  three  items  refer  to  the  contributions  of  Professor 
Kurz. 

(d)  The  robust  detection  techniques  developed  by  Professor  Kurz 
have  been  of  interest  to  the  Navy  in  solving  certain  types  of  under¬ 
water  sounding  problems.  His  former  students  were  involved  with  this 
problem  at  Bell  Laboratories  under  a  contract  with  the  US  Navy.  In 
recent  years,  one  of  his  students  (Dr.  Roger  Dwyer,  Underwater 
Sounding  Laboratory,  US  Navy,  New  London,  CT)  in  collaboration  with 
him  worked  on  applying  statistical  partition  techniques  to  underwater 
sounding . 

(e)  The  work  on  statistical  masking  techniques  and  image  recon¬ 
struction  algorithms  is  of  interest  to  the  IBM  Corp.  They  have  pub¬ 
lished  three  internal  research  reports  based  on  Professor  Kurz'  work  in 
this  area,  and  they  hope  to  include  these  approaches  in  their  future 
software  packages. 
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(f)  At  the  recent  IFIP  Conference,  Professor  Kurz  was  ap¬ 
proached  after  his  presentation  by  people  from  several  agencies  of  the 
U.S.  Government.  The  Bureau  of  The  Census  and  the  Department  of 
Energy,  in  particular,  asked  him  to  give  presentations  to  their  staffs 
on  his  robustizing  approaches  to  statistical  inference. 
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C.  ROBUST  PROCEDURES  IN  ESTIMATION  AND  DETECTION  PROB- 
LEMS 


Professor  L.  Kurz 

Unit  IE1-3 

1.  OBJECTIVE(S) 

Considerable  interest  generated  in  robust  estimation  and  detection 
problems  motivated  the  development  in  recent  years  by  the  senior 
investigator  in  cooperation  with  his  doctoral  students  of  several  classes 
of  rank  and  non-rank  detection  procedures  as  well  as  robust  estimation 
procedures  which  use  nonparametric  preprocessors  to  achieve  a  high 
degree  of  robustness  in  the  presence  of  unknown  noise  and/or  severe 
extraneous  contamination  of  the  data.  The  main  objective  of  future 
research  is  to  concentrate  the  effort  on  classes  of  detectors  based  on 
mixed  statistics  which  preserve  excellent  qualities  of  rank  and  non-rank 
procedures;  also  included  in  the  effort  will  be  the  imposition  of  side 
constraints,  statistical  dependence,  sequential  modes  of  operation  and 
array  processing  of  the  data.  A  parallel  effort  will  include  extensions 
of  the  robustized  estimation  techniques  to  problems  involving  stochastic 
approximation  minimum  variance  least  squares  estimators  (SAMVLS)  and 
Kalman-type  estimation  and  filtering  problems.  The  interaction  of  the 
robust  detection  and  estimation  techniques  in  adaptive  systems  will  be 
studied.  Particular  attention  will  be  given  to  the  application  of  the 
procedures  to  problems  in  radar,  sonar  and  communications. 

2.  APPROACH 

The  statistical  properties  of  partition  detectors  and  their  variants 
are  well  understood.  A  concise  exposition  of  the  present  state  of  the 
art  may  be  found  in  reference  1.  The  mathematical  framework  intro¬ 
duced  therein  permits  a  logical  evolution  of  many  robust  detectors  with 
an  almost  inexhaustible  source  for  developing  such  detectors  useful  in 
various  applications.  In  addition,  two  approaches  introduced  in  a 
recent  paper  of  ours2  for  including  side  constraints  in  the  design  of 
robust  detectors  represent  a  major  step  in  our  ability  to  evolve  many 
flexible  procedures  for  robust  detection  subject  to  various  physical 
constraints.  Inclusion  of  statistical  dependence  and  sequential  modes  of 
operation2'3'4'10  adds  another  important  dimension  to  the  utility  of 
partition  detectors.  The  interaction  of  the  approaches  outlined  above 
represents  an  almost  inexhaustible  source  for  solving  numerous  inter¬ 
esting  and  useful  problems  involving  detection  of  weak  signals  in  an 
unknown  environment. 

In  a  recent  paper,5  a  new  class  of  non-rank  robust  detectors, 
which  are  related  to  the  generalized  quantile  detector1  was  considered. 
Unlike  the  non-rank  detectors  based  on  quantile  statistics  with  fixed 
scoring  vectors,  the  new  class  of  detectors  uses  a  random  scoring 
vector.  This  new  class  of  detectors  is  particularly  useful  in  a  noise 
environment  that  is  relatively  rapidly- varying. 
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The  approach  to  robust  recursive  estimation  problems  initiated  in 
references  6  and  7  with  its  extensions  presented  in  references  8  and  9 
represents  another  source  of  numerous  theoretical  and  practical  prob¬ 
lems  of  interest  in  radar,  sonar  and  communications.  Inclusion  of 
statistical  dependence  and  other  constraints  in  this  approach  will  yield  a 
tool  of  almost  limitless  potential  for  solving  practical  estimation  problems 
when  the  parameters  to  be  estimated  are  severely  corrupted  by  noise  or 
the  noise  is  unknown  or  varying.  Various  additional  forms  of  robust- 
ized  estimators,  including  p.d.f.,  c.d.f.,  and  score  estimators,  permit 
solution  of  numerous  practical  problems.  For  all  approaches  outlined 
above  stress  will  be  placed  on  procedures  which  are  applicable  in  real¬ 
time  processing  and  are  simple  to  implement. 

3.  PROGRESS 

A.  The  results  on  dependent  sampling  in  sequential  partition 
detectors  reported  in  reference  10  were  extended  to  a  more  general 
model  of  Markov-dependent  noise.11  This  paper  involves  proving 
certain  basic  theorems  which  point  to  the  fact  that  considerable  advan¬ 
tage  in  detection  efficiency  of  partition  detectors  is  achieved  if  one  uses 
the  Markov-dependent  structure  of  noise.  The  results  open  the  door  to 
a  whole  new  class  of  robust  detectors  which  utilize  the  structure  of 
dependence  of  the  underlying  noise  to  improve  the  efficiency  of  detec¬ 
tion. 


B.  The  advent  of  modern  high-speed  sampling  techniques  result 
in  additional  stress  on  the  problem  of  signal  processing  in  dependent 
noise  mentioned  in  the  previous  paragraph.  Since  in  many  systems  of 
practical  interest  the  knowledge  of  the  noise  statistics  is  either  inexact 
or  unspecified,  some  adaptive  techniques  which  are  frequently  ineffec¬ 
tive  or  too  complicated  to  implement  were  suggested  in  the  past.  How¬ 
ever,  if  one  applies  the  methodology  of  system  identification,  the  ap¬ 
proach  is  essentially  equivalent  to  the  operation  of  a  prewhitening 
filter.  In  particular,  in  a  recent  paper12  a  procedure  for  robust  iden¬ 
tification  of  the  autoregressive  (AR)  model  was  introduced.  The  robust- 
izing  process  is  in  the  form  of  a  robustized  Robbins-Monro  stochastic 
approximation  (RMSA)  algorithm  and  is  based  on  the  m- interval  poly¬ 
nomial  approximation  method.  The  resulting  algorithm  represents  a 
recursive  robustized  version  of  the  well-known  maximum  entropy  method 
(MEM)  for  spectral  estimation  introduced  by  Burg13  or  of  the  popular 
Widrow  least-mean-square  (LMS)  adaptive  filter.14  Furthermore,  the 
robustized  algorithm  leads  naturally  to  a  robustized  Akaike's  information 
criterion  (AIC).1S  The  simplicity  of  implementation  and  flexibility  make 
the  application  of  the  new  robust  identification  algorithm  particularly 
attractive  in  practical  applications.  The  flexibility  and  robustness  of 
the  new  procedure  were  confirmed  by  extensive  Monte  Carlo  simulations. 

C.  The  standard  Kalman  filter  will  usually  diverge  when  the 
dynamic  noise  is  not  Gaussian.  A  robust  Kalman  filter  based  on  m- 
interval  polynomial  approximation  (MIPA)  method  for  nongaussian  noise 
was  formulated.16  Two  situations  were  studied  in  detail:  (a)  the 
states  ax  aaussi'  and  the  observation  noise  is  nongaussian,  (b)  the 
states  are  .»  ngt  jsian  and  the  observation  noise  is  gaussian.  It  was 
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shown  that  the  MI  PA  Kalman  filter  maintains  efficient  performance  over 
a  broad  class  of  noise  distributions  and  is  both  computationally  attrac¬ 
tive  and  easy  to  implement.  Compared  to  the  min-max  Kalman  filters 
suggested  by  other  authors,  the  new  class  of  filters  is  more  efficient 
and  robust.  The  theoretical  results  were  verified  by  extensive  Monte 
Carlo  simulations. 

D.  A  particularly  efficient  and  flexible  new  approach  to  robust 
detection  was  introduced  in  reference  17.  In  the  paper,  the  Huber- 
Tukey  model  of  mixture  distributions  is  replaced  by  a  more  versatile 
one  which  leads  to  solutions  which  are  more  appealing  from  a  practical 
point  of  view  than  the  robust  detectors  using  the  Huber-Tukey  noise 
model  in  conjunction  with  the  min-max  criterion  of  performance.  The 
new  detector  consists  of  two  parts  —  parametric  and  nonparametric  -- 
which  are  switched  in  or  out  depending  on  one  parameter  of  the  mixture 
distribution  model.  As  part  of  the  detector,  an  efficient  running  esti¬ 
mator  of  this  unknown  parameter  is  included.  When  compared  to  the 
min-max  detector  proposed  by  other  researchers,  the  new  detector  is 
more  efficient  for  deterministic  and  stochastic  signals.  This  divergence 
in  efficiencies  is  particularly  apparent  under  severe  and/or  variable 
noise  conditions. 

E.  Suppression  of  impulsive  noise  and  timing  jitter  to  robustize 
data  transmission  systems  requires  modification  of  signal  selection  and 
equalization  processes.18'19  A  further  contribution  to  this  robustizing 
approach  to  equalization  in  data  transmission  systems  was  presented  in 
two  new  papers.20'21  In  particular,  the  correlatively-coded  signals 
tolerant  to  sampling  offsets  (timing  jitter)  described  in  reference  19 
were  recursively  equalized  at  the  receiver.20  It  was  shown  that  for 
data  transmission  systems  disturbed  by  noise,  controlled  intersymbol 
interference  and  sampling  offsets,  the  recursive  equalizer20  yields 
improvement  when  compared  to  the  nonrecursive  equalizers. 19  In 
addition,  the  results  of  reference  18  were  generalized  to  PAM  systems 
using  partial  response  techniques.21  The  results  presented  in  the 
paper  permit  a  meaningful  study  of  practical  systems  because  even  if 
the  nonlinearity  is  not  added  to  the  system  to  suppress  impulsive  noise, 
all  real  systems  have  a  finite  dynamic  range  of  linear  operation  and  act 
as  soft  limiters  with  a  high  saturation  level.  In  particular,  the  effect 
of  nonlinearities  on  intersymbol  interference  and  signal  shaping  was 
considered.  The  generalization  of  the  equalizer  design  in  the  frequency 
domain  was  presented  resulting  in  suppression  of  degradation  of  system 
performance  because  of  the  nonlinearity  in  the  system.  Several  partial 
response  signals  which  are  of  interest  in  practical  systems  were  con¬ 
sidered.  Finally,  comparisons  between  a  standard  PAM  system  and  a 
partial  response  system  in  the  presence  of  nonlinearities  were  given. 

F.  The  importance  of  data  links  discussed  in  the  preceding  sec¬ 
tion  in  computer-communication  networks  led  naturally  to  the  preliminary 
study  of  transition  probabilities  in  such  networks.22  To  facilitate  the 
generation  of  the  pertinent  transition  probabilities,  an  approximate 
efficient  technique  was  evolved  in  the  paper.  Using  time-dependent 
coefficients  and  state-dependent  orthogonal  functions  in  a  generalized 
Fourier  series,  approximate  solutions  to  the  pertinent  Fokker-Planck 
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equation  were  generated  yielding  both  transient  and  steady-state  solu¬ 
tions  to  the  problem.  For  a  given  mean-square  error,  the  efficiency  of 
the  approximation  is  achieved  by  an  appropriate  choice  of  orthogonal 
functions  which  require  the  minimum  number  of  terms  in  the  series 
expansion.  The  small  number  of  coefficients  in  the  series  are  deter¬ 
mined  as  solutions  of  a  system  of  ordinary  differential  equations  result¬ 
ing  in  a  high  computational  efficiency  of  the  method.  It  should  be 
noted  that  the  procedure  outlined  in  the  paper  is  applicable  in  obtaining 
solutions  to  other  than  Fokker- Planck  partial  differential  equations. 

G.  The  theory  of  partition  detectors  was  extended  to  so-called 
double-threshold  detectors  which  are  particularly  useful  in  radar  appli¬ 
cations.23  Quantization  loss  of  the  commonly-used  double- threshold 
detector  is  minimized  by  using  a  multi-level  threshold  for  the  first 
threshold.  The  overall  performance  is  improved  by  optimal  selection  of 
partitions  (thresholds)  and  scores.  Expressions  for  false  alarm  rates 
and  detection  probabilities  and  upper  bounds  based  on  the  generalized 
Chebysher  inequality  are  presented.  In  radar  applications  the  new 
detector  performs  better  than  detectors  suggested  by  other  research¬ 
ers.  The  new  detector  has  potential  for  application  to  communications 
through  dispersive  media. 
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